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ABSTRACT 


The  study  of  impurities  and  defect  levels  in 

semiconductors  is  of  major  scientific  and 

technological  interest.  In  many  cases,  the 

phenoinona  associated,  particularly  with  deep  levels, 

are  not  well  understood  even  in  this  day.  In  the 

programs  supported  under  this  contract,  a  number  of 

major  contributions  were  made  in  this  broad  area  of 

research.  We  developed  a  new  form  of  spectroscopy 

which  allows  one  to  explore  some  of  the  excited 

states  of  various  centers.  We  were  the  first  to 

identify  and  study  in  detail  a  change  in  the 

structure  of  a  defect  due  to  a  radiation  with  a 

laser,  the  so-called  non- rad i a t i ve  induced  defect 

reaction  process.  We  studied  the  properties  of  some 

very  interesting  deep  levels  in  silicon  and,  in 

particular,  attempted  to  develop  an  understanding  of 

the  origin  of  these  deep  levels  which  are  thought  to 

be  associated  with  clusters  of  one  of  the  classic 

shallow  acceptors  and  iron.  We  developed  a  new 

t' 

theory  of  the  undulation  spectra  that  occurs  in  them 
and  donor  acceptor  pair  luminescence.  Finally,  we 
developed  a  theory  of  the  vibrational  properties  of 


defects  in  sem i conduct ors  and  applied  it  to  the 
study  of  oxygen  and  gallium  phosphide,  one  of  the 
better  known  systems  and  the  one  studied  in  our 
experiments  on  defect  reactions. 


1.  INTRODUCTION 


One  of  the  key  theoretical  and  experimental  problems 
in  modern  semiconductor  physics  is  the 
identification  and  description  of  the  properties  of 
deep  levels.  Deep  levels  are  typically  associated 
with  chemical  or  structural  defects  which  may  be 
inadvertently  introduced  into  a  perfect  crystal 
during  growth  of  processing.  Many  of  the  unsolved 
problems  in  semiconductor  technology  revolve  around 
our  ability  to  identify  and  control  these  dopants. 

In  the  program  we  made  a  number  of  major 
contributions  in  a  generic  sense  to  this  subject. 


II.  MAJOR  RESULTS 


We  developed  a  new  form  of  spectroscopy  called 
selective  excitation  spectroscopy  in  which  one  uses 
a  dye  laser  system  to  selectively  pump  the 
semiconductor  below  band  gap  and  then  detect  the 
radiation  emitted  by  the  semiconductor.  This 


spectroscopic  technique  allows  one  to  selectively 
pump  one  or  more  of  the  defect  levels  in  the 
semiconductor  and,  hence,  obtain  a  very  sensitive 
determi nation  of  their  properties.  This  new 
technique  is  described  in  detail  in  publications 
number  1  and  2.  (The  references  for  these 
publications  are  contained  at  the  end  of  this 
report,  and  an  appendix  to  the  report  includes  a 
copy  of  the  manuscript.) 

One  of  the  major  problems  in  electro-optical  devices 
is  the  degradation  of  the  light  output  during  the 
operation  of  the  device  due  to  changes  in  the  defect 
structure  in  the  semiconductor  which  have  been 
induced  by  the  presence  of  electrons  and  holes. 

While  this  process  has  been  speculated  as  being  the 
origin  of  a  number  of  degradation  mechanisms  for 
lasers  and  light  emitting  diodes,  no  detailed  basic 
scientific  study  has  been  carried  out  in  the  past  to 
show  that  such  a  process  can  actually  occur  for  a 
well-known  center  and  an  attempt  to  identify  its 
mechanism.  In  a  series  of  papers,  see  publications 
numbers  3  and  4,  we  carried  out  the  first  study  of 
such  a  defect  reaction.  We  chose  to  study  the 
well-known  center  and  GaP  associated  with  (Zn,0) 
pairs.  This  center  has  been  extensively  studied  in 


the  past  and,  hence,  is  a  likely  candidate  for 


analysis  for  this  kind  of  phenomona.  By  pumping  the 
luminescence  from  this  system  and  studying  the 
changes  and  the  line  shape  we  were  able  to  show  that 
defect  reactions  did  in  fact  occur. 

Nearest-neighbor  (Zn,0)  pairs  were  found  to  be 
broken  up  by  the  pumping  of  the  above-band-gap 
radiation.  To  our  knowledge,  this  is  the  first  such 
clear  indication  that  defect  reactions  could  in  fact 
occur . 

One  of  the  more  interesting  topics  in  modern 
semiconductor  defectolgy  is  the  question  of  whether 
defects  and  clusters  together  can  be  broken  up  in 
various  processes,  such  as  the  one  described  above. 
One  of  the  centers  that  has  attracted  a  great  deal 
of  interest  recently  is  one  originally  observed  at 
Caltech,  which  occurs  in  indium  doped  silicon 
samples  and  is  thought  to  be  due  to  the  presence  of 
indium  complexed  with  iron.  As  part  of  this 
program,  we  carried  out  a  study  of  the  role  of  iron 
in  the  formation  of  these  clusters  and  silicon  doped 
with  indium  and  silicon  doped  with  thalium.  We 
found  that  the  introduction  of  iron  led  to  the 
formation  of  the  cluster  defect  that  was  responsible 
for  the  luminescence,  and,  hence,  gave  the  first 
clear  indication  that  these  centers  actually 
involved  iron.  The  results  of  this  study  are 


reported  in  publication  number  5.  Two  other 
projects  included  the  first  theoretical  explanation 
of  undulation  spectra  that  did  not  involve  in  detail 
calculations  of  the  crystal  structure.  We  showed 
that  periodicities  in  the  luminescence  spectra  would 
result  from  a  simple  mathematical  argument  that  is 
applicable  any  time  you  have  a  periodic  structure. 
This  is  reported  in  publication  number  6.  We  also 
studied  the  properties  of  the  photol umi nescense  from 
Si-rich  SiGe  alloys  to  try  to  track  the  deep  levels 
occurring  in  this  system.  The  result  of  this  study 
is  published  in  publication  number  7. 

Finally,  we  developed  a  new  vibrational  theory  of 
oxygen  in  GaP.  Green's  function  techniques  were 
used  to  calculate  the  local  mode  properties  of  this 
important  defect  center.  After  careful  calculation, 
we  were  able  to  identify  many  of  the  modes  that 
appeared  in  the  luminescence  experiments  that  we 
carried  out.  The  details  of  these  modes  assisted  us 
in  our  understanding  of  the  defect  reaction  for  zinc 
oxygen  pairs  in  GaP  which  are  the  results  of  highly 
heating  some  of  the  local  modes  of  this  center.  The 
results  of  this  study  are  reported  in  publications 
number  8  and  9. 

In  summary,  this  program  was  a  very  successful  one 
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Selective  excitation  luminescence  in  bulk-grown  GaAs 

A.  T.  Hunter  and  T.  C.  McGill 

California  Institute  of  Technology,  Pasadena,  California  91125 
(Received  6  July  1981;  accepted  for  publication  2  October  1981) 

We  have  measured  the  excited  state  levels  of  two  different  shallow  acceptors  in  bulk-grown  GaAs, 
using  selective  excitation  luminescence.  The  IS- IS  energy  differences  were  measured  to  be  2 1.5 
and  18.5  meV,  respectively.  By  comparing  these  values  to  those  measured  by  two-hole  transition 
luminescence  in  high  quality  epitaxial  GaAs  [Ashen  etal.,  J.  Phys.  Chem.  Solids  36, 1041  (1975)], 
the  acceptors  were  identified  as  Zn  and  C.  The  measured  \S-2P  energy  differences  also  support 
the  identification.  These  studies  demonstrate  that  selective  excitation  luminescence  can  be  used  to 
identify  shallow  acceptors  in  bulk-grown  semi-insulating  GaAs,  and  hence  can  be  used  as  a 
diagnostic  tool  for  bulk-grown  samples. 

PACS  numbers:  78.55.Ds,  71.55.Fr,  78.60.  —  b 

tor  as  Zn  for  the  first  sample,  and  C  for  the  second,  as  will  be 
shown  below. 

The  experiment  is  essentially  that  used  for  ZnSe  in  Refs. 
4  and  5.  A  tunable  laser  is  used  to  create  an  electron  on  a 
donor  and  a  hole  on  an  adjacent  acceptor,  with  the  hole  in  an 
excited  state.  In  principle,  pairs  should  be  created  with  do¬ 
nor  electrons  in  excited  states  as  well,  but  lines  due  to  this 
process  are  not  observed.  The  first  excited  state  of  the  donor 
is  only  ~  1.5  meV  from  the  band  edge,  and  these  weakly 
bound  states  may  be  thermally  ionized  with  high  probability. 
Also,  the  very  extended  wavefunction  of  the  excited  donor 
may  be  substantially  broadened  by  interaction  with  other 
defects.  This  could  lead  to  a  line  breadth  which  is  too  wide  to 
resolve  against  background  luminescence. The  energy  re¬ 
quired  to  create  the  pair  with  an  excited  hole  is  6wlMir  =  Et 
—  E  J  —  Ed  +  e2/eR  +  J  *[R  ),  where  £s  is  the  band  gap, 

E  J  the  energy  of  the  excited  hole  with  respect  to  the  valence 
band,  ED  the  energy  of  the  electron  with  respect  to  the  con¬ 
duction  band,  e2/eR  is  the  Coulomb  attraction  of  the  ionized 
centers  at  distance  R  from  one  another,  and  J  *(R  )  represents 
the  interaction  of  the  donor  and  excited  acceptor  wave  func¬ 
tions.  This  is  illustrated  in  Fig.  1,  where  the  energy  of  the 
occupied  pair  is  plotted  against  the  inverse  of  the  pair  separa¬ 
tion  with  J  *(R )  taken  to  be  zero.  The  laser  will  create  pairs 
with  the  hole  in  the  first  excited  state  at  pair  separation  R ,, 
and  holes  in  the  second  excited  state  at  separation  R:.  The 
holes  very  quickly  relax  to  their  respective  ground  state,  then 
recombine  with  the  electrons,  emitting  light  at 

=  E,  —  Ea  —  Ed  +  e~/eR  +  J(R),  where  £,  is  now  the 
energy  of  the  ground  state  of  the  hole  with  respect  to  the 
valence  band,  and  /  (R )  is  the  interaction  between  the  donor 
and  the  ground  state  acceptor  wave  function.  The  energy 
difference  between  the  laser  and  emitted  light  consists  of  two 
terms  in  this  experiment:  AE  =  (£,  -  £$)  4-  [J *(R  ) 

—  J(R  )].Both/(R  )and/*(R  ) goto zeroforsufficiently sepa¬ 
rated  pairs,  and  in  this  limitd£  is  just  the  difference  between 
the  ground  and  excited  state  energy  of  the  acceptor  wave 
fvinctions. 

The  equipment  and  experimental  procedure  used  are  as 
follows.  A  Coherent  CR-3000K  ion  laser  was  used  to  pump  a 
Coherent  model  590  dye  laser,  with  a  carbo-cyanine  dye. 
The  dye  laser  was  tuned  with  a  3-plate  birefringent  filter 
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The  identification  of  shallow  acceptors  in  GaAs  by  the 
use  of  straightforward  photoluminescence  is  fairly  well  ad¬ 
vanced  for  high  quality  epitaxial  samples. 1  However,  there 
are  additional  problems  in  bulk-grown  GaAs  such  as  high 
defect  concentrations  and  low  luminescence  intensity  that 
limit  the  usefulness  of  these  methods.  Selective  excitation 
luminescence,2'5  in  which  a  tunable,  below  band-gap  laser  is 
used  to  selectively  excite  the  shallow  impurities  of  interest, 
seems  to  overcome  the  problems  in  identifying  these  impuri¬ 
ties  in  bulk-grown  GaAs.  This  technique  has  been  used  suc¬ 
cessfully  in  bulk-grown  InP,23  and  also  has  been  used  in 
ZnTe  (Ref.  2)  and  ZnSe.° 

In  this  letter,  we  report  for  the  first  time  results  for 
GaAs.  We  were  able  to  conclusively  identify  the  principal 
acceptor  in  two  bulk-grown  samples,  neither  of  which  were 
intentionally  doped  with  shallow  acceptors.  Both  samples 
studied  were  grown  by  the  liquid  encapsulated  Czochralski 
method.6  Sample  1  was  chrome  doped,  and  is  semi-insulat- 
tng.&  Sample  2  was  not  intentionally  doped,  and  is  p  type 
with  .V,  -  ,VD  —  tO16  cm-3.6  We  have  identified  the  accep- 


O  i/t»j  >/o, 


I /PAIR  SEPARATION 

FIG.  I .  Schematic  of  the  selective  excitation  luminescence  process.  ‘  Ener¬ 
gy  required  to  create  a  donor-acceptor  pair  is  plotted  against  the  inverse  to 
the  pair  separation  The  dye  laser  creaies  pairs  wiih  I  lie-  hole  ill  tin  I'lt’ilrd 
state.  The  hole  quickly  relaxes  to  the  IA  ground  state,  then  radiatively  re¬ 
combines  with  the  electron.  The  difference  between  the  energy  of  the  laser 
and  that  of  (he  emitted  light  provides  (he  ground-stute-excited  state  energy 
difference. 


placed  in  the  dye-laser  cavity.  This  tuning  element  allowed 
linewidths  of  between  0. 1  and  0.2  meV,  which  allowed  ener¬ 
gy  resolution  of  better  than  0.1  meV  in  most  cases,  because 
the  line  shapes  were  symmetric  Gaussians,  the  peaks  of 
which  could  be  determined  with  high  accuracy.  The  wave¬ 
length  of  light  used  for  the  experiments  was  such  that  the 
ground-state  emission  was  near  the  peak  of  the  donor-accep¬ 
tor  band.  The  laser  light  was  then  focused  onto  a  sample, 
which  was  immersed  in  liquid  helium  pumped  below  the 
lambda  point  in  a  Janis  cryogenic  dewar.  The  emitted  light 
was  analyzed  by  scanning  a  Spex  1404  double  grating  spec¬ 
trometer  and  detected  with  a  GaAs  photomultiplier  tube. 
The  gratings  used  were  replica  gratings,  but  ghost  and  back¬ 
ground  intensity  proved  to  be  insignificant  at  the  wavelength 
differences  of  interest  in  this  experiment.  Luminescence 
background  was  not  insignificant,  but  could  be  easily  distin¬ 
guished  from  the  scattering  processes  we  wished  to  observe. 
This  was  done  by  taking  additional  spectra  with  the  pump 
laser  at  slightly  different  wavelengths.  Those  features  de¬ 
scribed  above  shifted  in  energy  with  the  laser,  while  back¬ 
ground  luminescence  did  not. 

Figure  2  shows  luminescence  spectra  for  the  two  bulk- 
grown  GaAs  samples  used  in  this  study.  In  the  upper  spec¬ 
trum,  the  signal-to-noise  ratio  is  not  large  enough  to  resolve 
the  bound-exciton  line.  In  the  bottom  spectrum,  it  can  just 
barely  be  resolved  on  an  expanded  scale  at  1512  meV.  The 
two-hole  transition78  of  the  bound-exciton  line  (which  leaves 
the  remaining  hole  in  an  excited  state)  cannot  be  resolved  at 
all.  In  principle,  the  donor-acceptor  pair  luminescence 
bands  centered  at  1490  meV  can  be  used  to  identify  the  ac¬ 
ceptor  involved  in  the  transition.1  However,  this  line  has  a 
rather  high  intrinsic  width,  as  does  the  free-electron  to 
bound-hold  band  seen  at  higher  temperatures.  It  is  difficult 
to  differentiate  between  acceptors  with  ionization  energies 
that  differ  by  only  a  few  milli-electron  volts,  especially  when 
the  signal-to-noise  ratio  is  as  poor  as  it  is  in  the  upper  spec¬ 
trum.  Even  with  sufficiently  good  signal-to-noise,  line  shifts 
due  to  high  impurity  concentrations,  which  are  quite  possi¬ 
ble  in  bulk-grown  material,  make  identification  by  the  posi¬ 
tion  of  the  donor-acceptor  band  questionable.1 


ENERGY  (meV) 

FIG.  2.  Photoluminescence  spectre  for  two  samples  of  bulk-frown  GaAs. 
with  above  band-gap  excitation.  Intensity  of  emitted  light  is  plotted  on  the 
vertical  scale,  against  photon  energy  on  the  horizontal  scale. 


Selective  excitation  luminescence  data,  however,  pro¬ 
vides  an  unambiguous  identification  of  the  acceptors.  The 
top  three  spectra  in  Fig.  3  are  from  sample  1 ,  each  taken  with 
the  dye  laser  at  different  photon  energy.  The  spectra  show 
intensity  of  luminescence  plotted  against  energy  difference 
between  the  laser  and  luminescence.  This  means  the  zero  of 
energy  in  each  case  is  the  laser  energy  indicated  above  the 
spectrum  at  the  right.  Four  peaks  are  apparent  which  remain 
at  an  almost  constant  energy  from  the  laser.  The  most  in¬ 
tense  peak  is  located  ~  21.5  meV  below  the  laser  line,  and  is 
seen  in  all  three  spectra.  The  line  does  shift  slightly  as  the 
energy  of  the  pump  laser  is  changed,  which  was  determined 
by  carefully  measuring  the  energy  differences  for  the  spectra 
shown  here,  as  well  as  for  other  spectra.  This  shift  is  due  to 
the  interaction  between  the  donor  and  acceptor  wave  func¬ 
tions  described  earlier,  and  causes  the  energy  difference  to 
decrease  as  the  pair  separation  is  decreased.  The  shift  is  less 
than  0. 1  meV  (which  is  the  uncertainty  of  the  energy  separa¬ 
tion  measurement  for  this  sample)  for  pump  laser  energies  of 
1509.7  meV  and  smaller.  This  pump  energy  creates  pairs 
with  approximately  200-A  pair  separation.  At  a  pump  ener¬ 
gy  of  1512.6  meV,  corresponding  to  ~  130-A  pair  separa¬ 
tions,  the  energy  shift  is  0.5  meV.  In  addition  to  this  fairly 
intense  peak,  three  other  less  intense  peaks  can  be  resolved, 
though  not  in  all  the  spectra.  They  occur  at  19.2,  21.5,  and 
25.0  meV.  Interpreting  the  data  as  was  done  in  Refs.  2  and  3, 
the  principal  line  is  assigned  to  the  15-25  energy  difference, 
and  the  less  intense  lines  are  assigned  to  the  IS-2P 
differences. 

These  results  can  now  be  compared  to  theoretical  calcu¬ 
lations,9  and  to  measurements  made  on  epitaxial  GaAs  by 
other  methods.  These  earlier  results,  and  the  results  of  this 
study,  are  tabulated  in  Table  I.  The  two-hole  shift  measured 
for  Zn-doped  epitaxial  samples  is  21.8  meV'-8  and  is  very 
close  to  the  2 1 . 5-me V  15-25  difference  measured  for  sample 
1.  The  position  of  the  less  intense  peaks  line  up  well  with  the 


FIG.  3.  Selective  excitation  lu¬ 
minescence  spectre  for  two 
samples  of  bulk-grown  GaAs. 
Intensity  of  emitted  light  is 
plotted  against  the  difference 
in  photon  energy  of  emitted 
light  and  that  of  the  laser.  Pho¬ 
ton  energy  of  the  laser  for  each 
spectrum  is  given  over  the 
spectrum  at  the  right.  The  top 
three  spectre  are  for  sample  1. 
the  bottom  four  for  sample  2. 
The  labels  for  the  acceptor  lev¬ 
els  are  those  used  in  Ref.  9. 
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TABLE  I.  Ground  stale  to  excited  state  splittings  for  the  shallow  acceptors  C  and  Zn  in  GaAs.  The  labels  for  the  acceptor  levels  are  those  used  in  Ref.  9. 


l-Si/j-iPi/j 

(meV) 

15„j-25„, 

(meV) 

\S%/2—2Pi vj(/* J 
(meV) 

(meV) 

Zn 

Theory* 

19.3 

23.1 

23.5 

25.4 

Photoconductivity  * 

19.4 

... 

23.2 

25.0 

"Two-hole”  shift* 

... 

21.8 

... 

... 

This  work 

19.2 

21.5 

23.1 

25.0 

C 

Theory* 

14.6 

18.4 

18.8 

20.7 

Photoconductivity  ‘ 

15.2 

... 

19.3 

21.3 

"Two-hole"-shift* 

— 

18.5 

... 

... 

This  work 

15.1 

18.5 

... 

... 

‘Deduced  from  theoretical  values  quoted  in  Ref.  9  by  replacing  the  calculated  15  position  (25.7  meV)  with  that  measured  in  Ref.  1 : 26.0  meV  for  C  and  30.7 
meV  for  Zn. 

"Measured  on  epitaxial  GaAs  from  Ref.  10. 

‘Measured  on  epitaxial  GaAs  from  Ref.  1. 

15-2/*  differences  measured  using  photoconductivity10  on  Therefore,  both  the  breadth  of  the  lines  and  the  increased 

epitaxial  GaAs.  Also  included  in  the  table,  for  reference,  are  intensity  of  the  15-2/*  lines  in  sample  1  could  be  explained  by 

the  energy  differences  deduced  from  the  effective  mass  cal-  a  much  higher  ionized  impurity  concentration  in  this  sample 

culation  of  Baldereschi  and  Lipari,9  by  replacing  the  calcu-  than  in  sample  2.  1 

lated  15  position  with  that  measured  by  Ashen  et  al. 1  The  In  conclusion,  we  have  measured  the  excited-state  lev- 

matchup  between  our  results  and  both  theory  and  earlier  els  of  shallow  acceptors  in  bulk-grown  GaAs,  using  selective 
results  indicate  that  the  acceptor  in  this  sample  is  Zn.  It  also  excitation  luminescence.  While  these  levels  have  been  mea- 
demonstrates  that  selective  excitation  luminescence  pro-  sured  in  GaAs  before,  it  was  done  using  different  methods  in 

vides  an  unambiguous  method  of  determining  the  identity  of  high  quality,  epitaxial  layers.  This  is  the  first  time  the  levels 
shallow  acceptors  in  bulk-grown  semi-insulating  GaAs.  have  been  measured  in  bulk-grown  GaAs,  where  the  other 

We  have  also  used  this  method  to  identify  the  acceptor  methods  are  not  as  useful.  In  doing  so,  selective  excitation 

in  the  second  sample  as  carbon.  The  bottom  half  of  Fig.  3  luminescence  has  been  shown  to  be  useful  in  identifying  the 

shows  four  spectra  taken  on  this  sample,  each  at  a  different  shallow  acceptors  in  bulk-grown  GaAs. 

laser  energy.  Only  two  peaks  are  visible  on  these  spectTa,  but  We  would  like  to  acknowledge  the  support  from  the 

they  are  sufficient  to  identify  the  acceptor.  The  principal  line  Office  of  Naval  Research  under  Contract  No.  N00014-81- 
in  this  case  is  18.5  meV  below  the  laser  energy  (and  is  some-  K-0305.  The  authors  also  gratefully  acknowledge  H.  Ki- 
what  narrower  than  the  principal  line  in  sample  1)  The  18.5-  mura,  O.  J.  Marsh,  and  H.  V.  Winston  of  Hughes  Research 

meV  difference  is  again  the  energy  difference  extrapolated  to  Laboratories  for  providing  the  samples  and  the  results  of 

the  distant  pair  limit,  which  was  determined  by  making  the  Hall-effect  measurements  on  the  samples, 
measurement  for  a  wide  range  of  pump  energies.  The  18.5- 
meV  value  is  the  two-hole  shift  measured  for  carbon-doped 
epitaxial  samples,11  and  identifies  the  dopant  in  our  sample 
as  carbon.  The  second  line,  at  15.2  meV,  confirms  the  identi¬ 
fication,  because  this  is  the  \S~2Pin  shift  measured  for  car¬ 
bon  using  photoconductivity. 10  These  results  are  also  tabu¬ 
lated  in  Table  I.  ‘D.  J.  Ashen,  P.  J.  Dean.  D.  T  J.  Hurle,  J  B.  Mullin,  A.  M.  White,  and  P. 

Two  points  concerning  line  breadth  and  height  could  ,D  Greene,  J.  Phys.  Chem.  Solids  36,  i04i  i  I97SI. 

use  further  discussion.  The  increased  linewidth  in  sample  l  ’*’££**'  °  1  Robbins’  *nd  S  G  Blshop' Sol,d  s,ate  Commun- 32- 379 

could  be  indicative  of  much  higher  ionized  impurity  concen-  >P.  j.  Dean.  D.  j.  Robb, ns.  and  s.  G.  Bishop,  J.  Phys.  c  12.  5567 1 19791. 

tration.  The  electric  fields  of  nearly  charged  impurities  can  ‘H.  Tews  and  H.  Venghaus,  Solid  State  Commun.  30. 219  ( 19793. 

inhomogeneously  broaden  the  impurity  levels  through  a  sec*  ’H- Tews- H-  v«nghaus,  and  P.  J  Dean.  Phys.  Rev.  b  19, 5178  (1979). 

ond-order  Stark  shift. ' '  In  addition  to  shifting  the  impurity  ^  Winstonipnvatecommun.cat.onL 

energy  levels,  the  Stark  effect  also  mixes  5  and  P  wavefunc-  >A.  M.  white,  p.  1.  Dean,  D.  y  Ashen,  J.  B  Mullin,  M.  Webb,  B.  Day,  and 

tions.  Since  the  hole  must  undergo  what  is  apparently  a  par-  P.  D.  Greene,  1.  Phys.  c  6,  l  243 1 1973). 

ity  forbidden  relaxation  process2  to  reach  the  ground  state  Baldereschi  and  N  o  Lipari,  Phys.  Rev  b  9, 1525  (I974|. 

from  a  /Mike  excited  state,  this  mixing  would  enhance  the  ^8|K,rkman' R  A  s,radlin8' and  p  J  u"  Ch“»8- J  Phys.  c  11. 419 

relaxation  from  the  IP  excited  states  to  the  15  ground  state.  "d.  m.  Larsen,  Phys.  Rev.  b  8, 535  (19731. 
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We  observe  for  the  first  time  two  excited  states  of  the  78-meV  acceptor  in  liquid  encapsulated 
Czochralski  GaAs  with  excitation  energies  of  62.9  and  66.9  meV  above  the  ground  state, 
respectively.  These  levels  were  observed  in  two  different  samples  using  selective  excitation 
luminescence  and  electronic  Raman  scattering.  We  conclude  that  these  levels  could  correspond  to 
transitions  from  the  lr  ground  state  to  the  split  1j'2j'  excited  state  of  a  double  acceptor. 

PACS  numbers:  71.55.Fr,  78.50.Ge,  78.30.Gt,  78.55.Ds 


Recently  there  has  been  considerable  interest  in  using 
semi-insulating  GaAs  as  a  substrate  for  device  fabrication. 
Presently  such  material  is  grown  undoped  by  the  liquid  en¬ 
capsulated  Czochralski  (LEC)  method.  It  has  been  observed 
that  if  the  melt  stoichiometry  becomes  Ga-rich  during  the 
growth  process,  the  material  is  p  type. 1  Such  material  also 
possesses  an  acceptor  level  at  78  me  V  above  the  valence  band 
which  is  not  present  in  the  semi-insulating  material.* 3  It  has 
been  suggested  that  the  center  is  a  double  acceptor,  with  first 
and  second  ionization  energies  of  78  and  200  meV,  respec¬ 
tively.4-*  Although  some  p-like  excited  states  have  been  ob¬ 
served,3  they  do  not  give  a  strong  indication  of  the  valency  of 
an  acceptor  since  p-like  excited  states  are  nearly  identical  for 
simple  single  and  double  acceptors,  as  is  the  case  in  germani¬ 
um.7  The  r-like  excited  states,  however,  give  a  better  indica¬ 
tion  of  the  valency.  In  this  letter  we  report  the  observation  of 
two  previously  unknown  excited  states  of  the  78-meV  accep¬ 
tor  using  two  techniques  sensitive  to  the  detection  of  s-like 
excited  states:  selective  excitation  luminescence  and  elec¬ 
tronic  Raman  scattering.  Our  results  from  both  experiments 
indicate  that  this  level  is  due  to  a  double  acceptor. 

Three  different  LEC  samples  were  used  for  this  study. 
Samples  1  and  2  were  both  grown  p  type  from  As-deficient 
melts.  Temperature-dependent  Hall  effect  measurements  on 
sample  1  indicated  l.lx  10,s  cm-3  donor,  5.3 X  10”  cm-3 
carbon,  and  8.1  x  10”  cm-3  78-meV  acceptor  concentra¬ 
tions.  Sample  2,  taken  from  a  crystal  similar  to  sample  1 ,  was 
neutron  transmutation  doped  and  annealed  at  580  *C  for  1  h. 
Neutron  transmutation  doping  added  3.8 X  10”  cm-3  Se 
and  Ge  donors,  plus  additional  damage  related  donors.  Tem¬ 
perature-dependent  Hall  effect  measurements  on  an  unirra¬ 
diated  sample  adjacent  to  sample  2  given  the  same  anneal 
gave  concentrations  of  2.6x  10”  cm-3  donors,  7.4x  10” 
cm-3  carbon,  and  3.6 X  10”  cm- 3  78-meV  acceptors.  Mea¬ 
surements  on  the  unirradiated  samples  showed  the  slope  for 
carbon  below  4.2  K,  indicating  that  the  78-meV  level  was 
fully  occupied."  Temperature-dependent  Hall  measure¬ 
ments  on  a  sample  transmutation  doped  and  annealed  in  a 
manner  identical  to  sample  2  showed  an  0.07 -eV  slope  how¬ 
ever,  indicating  that  the  78-meV  acceptor  was  partially  com¬ 
pensated  in  the  transmutation  doped  and  annealed  samples. 


The  78-meV  acceptor  was  also  identified  by  band-acceptor 
photoluminescence  in  samples  1  and  2.  Sample  3  was  chro¬ 
mium-doped,  semi-insulating  GaAs,  and  photolumines¬ 
cence  experiments  did  not  detect  the  presence  of  a  level  at  78 
meV  in  this  sample.  The  setup  is  nearly  identical  for  both 
experiments.  The  below  band-gap  excitation  was  provided 
by  a  Coherent  model  590  tunable  dye  laser  containing  a  car- 
bocyanine  dye,  using  a  Coherent  CR3000K  krypton  ion  la¬ 
ser  to  optically  pump  the  dye.  Tuning  the  dye  laser  was  ac¬ 
complished  via  a  three-plate  quartz  birefringent  filter.  The 
sample  was  mounted  in  a  Janis  cryogenic  dewar  and  im¬ 
mersed  in  liquid  helium  pumped  below  the  lamda  point.  The 
luminescence  from  the  sample  was  analyzed  by  a  Spex  1404 
scanning  spectrometer,  using  a  cooled  SI  photomultiplier  as 
a  detector. 

The  excited  states  of  the  78-meV  acceptor  were  mea¬ 
sured  using  two  different  techniques:  selective  excitation  lu¬ 
minescence  (SEL)  and  electronic  Raman  scattering  (ERS). 
The  SEL  experiment  is  essentially  that  which  was  described 
in  Ref.  9.  A  below  band-gap  photon  is  used  to  create  a  donor- 
acceptor  pair  by  exciting  an  electron  from  an  ionized  accep¬ 
tor  to  an  ionized  donor,  leaving  the  resulting  hole  on  the 
acceptor  in  an  excited  state.  The  excited  hole  rapidly  relaxes 
to  the  acceptor  ground  state,  and  eventually  recombines 
with  the  electron  on  the  donor,  emitting  a  photon.  If  we 
ignore  any  non-Coulombic  interaction  between  the  donor 
and  acceptor,  the  energy  difference  between  the  incident  and 
emitted  photons  is  the  difference  between  the  ground  state 
and  excited  state  binding  energies.  Changing  the  laser  energy 
alters  only  the  separation  of  the  donor  and  acceptor,  en¬ 
abling  us  to  identify  the  excited  states  since  they  remain  at 
fixed  energy  separation  from  the  laser.  Since  this  technique 
requires  both  the  donor  and  acceptor  to  be  ionized,  it  is  use¬ 
ful  for  studying  compensated  material.  Occupied  acceptors 
are  best  studied  using  ERS.  In  this  experiment  a  below  band- 
gap  photon  can  create  a  virtual  electron-hole  pair.  The  elec- 
t  ron-holc  pair  can  then  recombine,  raising  a  hole  bound  to  an 
acceptor  to  an  excited  state  and  emitting  a  photon  with  the 
remaining  energy.  The  energy  difference  between  the  inci¬ 
dent  and  emitted  photons  is  the  difference  between  the 
ground  state  and  excited  state  binding  energies. 
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FIG.  t  .  Selective  excitation  luminescence  spectra  from  sample  2.  Intensity 
of  emitted  light  is  plotted  against  the  difference  between  the  emitted  photon 
energy  and  that  of  the  laser.  Three  different  laser  energies  were  used  to  dis¬ 
tinguish  between  SEL  features  and  background  processes.  The  laser  photon 
energy  for  each  spectrum  is  at  the  right. 

The  SEL  results  for  sample  2  for  three  different  exciting 
laser  frequencies  are  shown  in  Fig.  1 ,  and  the  ERS  results  for 
sample  1  for  two  different  exciting  laser  frequencies  are 
shown  in  Fig.  2.  Two  transitions  are  evident  in  the  SEL  data, 
one  at  62.5  meV  (line  A)  and  one  at  66.5  meV  (line  B).  The 
broad  background  in  the  spectra  is  due  to  donor-acceptor 


FIG.2  .  Electronic  Raman  scattering  spectra  from  sample  I .  Intensity  of 
emitted  light  is  plotted  against  the  difference  between  the  emitted  photon 
energy  and  that  of  the  laser.  Two  different  laser  energies  were  used  to  distin¬ 
guish  between  ERS  features  and  background  processes.  The  laser  photon 
energy  for  each  spectrum  is  at  the  right. 


TABLE  I.  Observed  line  positions  and  interpretations. 


A 

B 

c 

ImeVr 

AE 

lmeV|K 

SEL 

62.5  meV 

66.5  meV 

4.0 

l-'KS 

62.*>  meV 

66  >*  nirV 

72. *< 

40 

Single 

acceptor*1 

6.V4  meV*1 

Wi.ft  meV 

3  2 

Double 

acceptor* 

l^—lr  '2s1 
l’5,75) 

Ir— lx  '2i' 

1  'S.'S) 

2.6 

‘  We  believe  that  line  C  is  due  to  Raman  scattering  from  2  LO  (f )  phonons, 
and  that  it  is  not  related  to  the  78-meV  level. 

*  Energy  separation  of  lines  A  and  B. 

‘Based  on  work  done  by  Baldereschi  and  Lipari  (Refs.  12  and  13). 

*  This  estimate  is  obtained  by  including  the  observed  central  cell  shift  for  the 
78-meV  acceptor  (see  Ref.  14). 

*  Based  on  the  helium-atom-like  approximation  discussed  in  the  text. 

luminescence.  Three  peaks  are  visible  in  the  ERS  data:  line  A 
at  62.9  meV,  line  B  at  66.9  meV,  and  line  C  at  72.9  meV. 
Lines  A  and  B  appear  shifted  about  0.4  meV  toward  greater 
energy  loss  in  the  ERS  experiment  than  in  SEL.  This  may  be 
expected  since  any  non-Coulombic  interaction  between  the 
donor  and  acceptor  wave  functions  in  SEL  would  tend  to 
shift  the  measured  energy  from  the  actual  binding  energy. 
Therefore,  the  energies  obtained  in  ERS  should  be  more  ac¬ 
curate  measurements  of  the  excited  state  energies.  Line  C 
was  not  seen  in  the  SEL  experiment.  The  observed  line  posi¬ 
tions  are  summarized  in  Table  I. 

The  two  lines  A  and  B  show  up  only  in  material  con¬ 
taining  the  78-meV  acceptor.  Figure  3  compares  ERS  spec¬ 
tra  from  sample  1  with  spectra  from  sample  3,  which  does 
not  contain  the  78-meV  level.  There  is  no  indication  of  lines 
A  or  B  in  sample  3.  This  suggests  that  they  are  both  due  to 
the  78-meV  level.  Although  they  lie  in  the  energy  region  for 
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ENERGY  LOSS  <rr.eV) 

FIG.  3.  Electronic  Raman  scattering  spectra  from  samples  1  and  3.  Intensi¬ 
ty  of  emitted  light  is  plotted  against  the  difference  between  the  emitted  pho¬ 
ton  energy  and  that  of  the  laser  Narrow  lines  at  36. 5  and  33  .6  meV  are  the 
LOi/'l  and  TChfl  phonons,  respectively. 
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two-phonon  Raman  scattering,  we  can  be  sure  that  they  are 
electronic  in  nature  since  the  excitation  intensity  in  SEL  is 
too  low  to  observe  even  single-phonon  Raman  scattering. 
There  is  evidence  of  a  line  at  73  meV  in  sample  3,  which  may 
correspond  to  line  C  seen  in  sample  1 .  Since  wc  do  not  ob¬ 
serve  line  C  in  the  SEL  experiment,  and  it  appears  to  be 
present  in  sample  3,  we  believe  that  it  is  due  to  2  LO  ( r ) 
phonons  (room-temperature  Raman  scattering  measure¬ 
ments  give  a  value  of  72.9  meV  for  this  transition10)  and  not 
related  to  the  78-meV  acceptor. 

We  believe  lines  A  and  B  are  due  to  the  split  lr  '2s' 
(using  atomic  single  particle  notation)  state  of  a  double  ac¬ 
ceptor.  The  r-like  excited  states  of  a  double  acceptor  will  be 
split  into  two  levels  by  the  exchange  interaction.  A  simple 
model  of  a  double  acceptor  in  GaAs  can  be  obtained  by  ana¬ 
logy  to  the  single  acceptor  effective  mass  approximation.  We 
consider  the  holes  to  bey  =  3/2  particles,  with  mass  equal  to 
the  heavy  hole  effective  mass,  and  that  the  electric  fields  are 
screened  by  the  static  dielectric  constant.  For  a  heliumlike 
acceptor  in  this  approximation,  we  obtain  two  Is  '2s1  levels: 
a  tenfold  degenerate  [}S,yS )  level  (with  symmetric  spins),  and 
a  sixfold  degenerate  level  ( 'S.’S )  (with  antisymmetric  spins). 
We  can  scale  the  energy  splitting  in  helium  to  GaAs  in  this 
model  by  multiplying  the  energy  in  helium  by  the  ratio  of  the 
Bohr  radii  divided  by  the  dielectric  constant.  The  experi¬ 
mentally  observed  magnitude  of  the  splitting  in  helium  is 
0.79  eV,"  which  scales  to  2.6  meV  in  GaAs.  Including  the 
effects  of  the  crystal  field  and  the  spin-orbit  interaction  with 
the  valence  band  will  not  change  the  splitting  to  first  order, 
since  they  will  only  affect  the  wave  functions  used,  not  the 
potential.  The  observed  splitting  between  lines  A  and  B  is  4.0 
meV,  which  is  in  good  agreement  with  the  estimate  (see  Ta¬ 
ble  I). 

Another  possible  interpretation  of  lines  A  and  B  is  the 
following:  line  A  could  be  the  lSJ/2— »25J/2  transition  and 
line  B  could  be  the  lSin-*2PJ/2  transition  of  a  simple  accep¬ 
tor.  Theoretical  estimates  based  on  work  done  by  Balderes- 
chi  and  Lipari,'2-13  using  the  observed  value  of  78  meV  for 
the  15]/2  binding  energy,  yield  values  for  these  transitions  of 
63.4  meV  (including  the  central  cell  shift1*)  and  66.6  meV, 
respectively.  We  obtain  good  agreement  between  lines  A  and 
B  and  these  estimates.  However,  there  is  evidence  that  line  B 
is  not  due  to  a  /7-like  final  state,  but  rather  an  r-like  final  state. 
Previous  bEL  measurements  on  carbon  and  zinc  in  GaAs 
reveal  r-like  states  with  approximately  two  to  five  times  the 
intensity  of /7-like  states.’  Carbon  has  also  been  observed  by 
ERS  in  GaAs,”  and  the p-like  transitions  were  observed  to 
be  at  least  ten  times  less  intense  than  ther-like  transitions.  In 
the  experiments  conducted  here,  the  intensities  of  the  A  and 
B  lines  are  nearly  identical  (see  Figs.  1  and  2).  Also,  neither 
experiment  shows  evidence  of  the p-like  excited  states  seen  in 
IR  absorption.1  (In  both  SEL  and  ERS,  all  p-like  states  are 
seen  with  roughly  the  same  intensity.)  Finally,  IR  absorption 
measurements  heavily  favor  p-like  excited  states  over  5-like 
excited  states  since  the  latter  are  forbidden  to  first  order. 


Thus,  one  would  not  expect  to  seer-like  states  in  IR  absorp¬ 
tion.  If  line  B  were  actually  p-like,  it  should  have  been  ob¬ 
served  in  the  IR  absorption  experiment.3  This  leads  to  the 
conclusion  that  lines  A  and  B  are  both  5-like. 

Wc  also  have  evidence  of  a  line  at  4.8  meV  in  F.RS  data 
front  sample  I.  This  line  was  not  seen  in  any  other  samples, 
and  might  be  related  to  the  78-meV  acceptor.  This  line  is  not 
donor  related  since  the  donors  are  ionized  in  this  sample. 
This  line  could  correspond  to  a  transition  between  two  levels 
of  a  split  lr2  ground  state.  Such  splittings  have  been  observed 
for  mercury  and  zinc  in  germanium.716  [The  ground  state  of 
a  double  acceptor  in  the  effective  mass  approximation  is  six¬ 
fold  degenerate  ( lS,sS ),  and  can  be  split  by  interaction  with 
the  valence  band.17] 

In  conclusion,  we  have  observed  two  excited  states  of 
the  78-meV  acceptor  in  GaAs,  at  62.9  and  66.9  meV.  We 
believe  that  these  lines  correspond  to  transitions  from  the 
ground  state  to  5-like  excited  states  of  this  acceptor.  It  is 
likely  that  the  excited  states  involved  are  the  lr  '2s'  fS^S) 
and  Is  '2j‘  ('5, ’5)  states  of  a  double  acceptor.  Our  results 
show  that  SEL  and  ERS  are  good  techniques  for  determin¬ 
ing  the  valency  of  acceptors  since  they  are  sensitive  to  the  5- 
like  excited  states. 
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0305.  The  authors  also  gratefully  acknowledge  M.  H.  Young 
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used  in  this  study. 
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Defect  Reactions  in  GaP:(Zn,0) 
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We  have  observed  photoinduced  reactions  between  pairs  of  zinc  and  oxygen  impurities 
in  gallium  phosphide.  From  photoluminescence  studies  we  find  that  the  nearest-neighbor 
(Zn,  O)  pairs  are  dissociating,  after  which  they  re-form  as  further  separated  pairs.  The 
activation  energy  for  the  dissociation  is  found  to  be  0.60*0.07  eV  for  the  photoinduced 
reaction,  and  2.6*  0.6  eV  for  the  purely  thermal  reaction.  We  tentatively  identify  the 
photoinduced  reaction  as  being  due  to  excitation  of  local  phonon  modes  by  nonradiative 
electron-hole  recombination. 


PACS  numbers:  66.30.Jt,  63.20. Pw,  78.55.Ds 
Defect  complexes  in  semiconductors  can  play  a 
very  important  role  in  determining  the  electrical 
and  optical  characteristics  of  the  material.  One 
simple  complex  which  commonly  occurs  is  a  pair 
of  donor  and  acceptor  impurities.  The  energy  of 
the  luminescence  emitted  from  donor -acceptor 
(DA)  pairs  depends  on  the  separation  of  the  im¬ 
purities,1  and  so  a  luminescence  spectrum  pro¬ 
vides  information  about  the  number  of  pairs  of 
each  possible  separation. 

-  In  this  work  we  use  photoluminescence  to  moni¬ 
tor  the  relative  positions  of  Znc«  and  Op  impuri¬ 
ties  in  GaP.  When  the  GaP  is  subjected  to  laser 
excitation  at  high  temperatures,  we  observe  reac¬ 
tions  in  which  the  impurities  diffuse  through  the 
lattice  to  occupy  new  substitutional  sites.  This 
is  the  first  time  DA  luminescence  has  been  used 
as  a  microscopic  probe  of  defect  reactions. 

From  low-temperature  studies  we  identify  these 
reactions  as  being  the  dissociation  of  nearest- 
neighbor  ( Zn,  O)  pairs  and  the  formation  of  fur  - 
ther  separated  pairs.  Such  a  dynamic  process  of 
reactions  among  lattice  defects  of  known  initial 
and  final  chemical  type  has  not  previously  been 
observed  in  a  semiconductor.  The  photolumi¬ 
nescence  spectra  reported  here  provide  valuable 
information  about  those  phonon  modes  which  are 
likely  to  be  responsible  for  the  (Zn, O)  dissocia¬ 
tion.  Most  other  studies  of  defect  reactions  in 
semiconductors  have  not  included  any  observa¬ 
tions  of  the  vibrational  modes  of  the  defect. 

The  samples  used  in  this  work  were  liquid- 
phase -epitaxially  grown  pn  junctions  on  n-type 
GaP  (single  crystal).  The  growth  technique  is 
described  by  Saul,  Armstrong,  and  Hackett.2 
The  exposed  50- nm -thick  layer  of  the  junction 
is  p  type,  doped  with  Zn  and  O.  The  (substitu¬ 
tional)  O  concentration  is  about  1010  cm'1  and 
the  Zn  concentration  is  about  3  xlO17  cm"3 
(based  on  the  crystal-growth  conditions).  Other 
impurities  include  N  and  S.  The  photolumines¬ 
cence  spectra  were  obtained  using  above-band- 


gap  Ar*  laser  excitation  and  a  liquid-N2 -cooled 
S-l  type  photomultiplier. 

Figure  1  shows  several  high-  and  low-tempera¬ 
ture  photoluminescence  spectra.  In  Fig.  1(a)  the 


FIG.  1.  Photoluminescence  spectra  at  various  tem¬ 
peratures  ( T )  and  laser  power  densities  (f);  (a)  time- 
resolved  spectra  showing  defect  reactions  in  progress 

ami  fl>)  Iho  imliul  nnd  final  status  of  lliu  reaction.  Thu 
labels  m  denote  luminescence  from  an  mth-nearest- 
neighbor  (Zn.O)  pair.  NP  refers  to  a  no-phonon- 
assisted  transition.  Theoretical  predictions  for  the 
line  positions  are  shown  on  the  abscissa  (also  in  Fig. 
2). 
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time  between  successive  scans  is  about  30  min. 

We  see  that  the  spectra  are  changing  with  time. 
This  is  a  direct  observation  of  a  reaction  in 
progress.  We  refer  to  this  reaction  as  an  “agin£’ 
process.  The  initial  and  final  spectra  shown  are 
similar  to  those  seen  by  Dapkus  and  Henry.1 
These  authors  showed  that  the  aging  reaction  is 
due  to  photoinduced  destruction  of  nearest-neigh¬ 
bor  (Zn,  O)  pairs,  and  they  speculated  that  these 
pairs  were  in  fact  dissociating  and  forming  fur¬ 
ther  separated  pairs.  Here,  we  prove  that  this 
speculation  is  correct. 

To  identify  the  luminescence  bands  we  go  to  the 
low -temperature  spectra  shown  in  Fig.  1(b). 

From  the  work  of  previous  authors4'7  we  identify 
the  band  centered  at  1.82  eV  as  being  due  to 
nearest-neighbor  pairs  (Zn,  0)„«lf  and  the  band 
centered  at  1.35  eV  as  due  to  far-separated  pairs 
(ZntO)w_..  The  breadth  of  these  bands  is  due  to 
phonon-assisted  radiative  transitions.  The  (Zn,  S) 
DA  luminescence  at  about  2.05  eV  is  constant 
in  intensity  between  the  unaged  and  aged  material, 
which  indicates  that  the  low -temperature  minor¬ 
ity  carrier  concentration  is  also  constant.  Thus, 
the  reduced  intensity  of  the  m  =  1  band  in  the  aged 
material  must  be  due  to  a  reduction  of  the  (Zn.O), 
concentration.3  Similarly,  the  (Zn, con¬ 
centration  has  increased.  Centered  at  1.59  eV 
is  a  luminescence  band  which  is  seen  in  the  aged 
sample,  but  barely  visible  in  the  unaged  material. 
The  breadth  of  this  band  is  similar  to  those  of 

( Zn,  O)  i  and  (Zn.O)„, _ and  the  resolved  phonon 

structure  is  characteristic  of  a  complex  involving 


0.  Furthermore,  if  we  compare  the  no-phonon 
transition  location  with  theoretical  predictions,1’* 
we  see  that  the  position  of  this  new  band  is  as  ex¬ 
pected  for  ( Zn,  0)2.  We  thus  identify  the  new 
luminescence  as  being  due  to  second-nearest- 
neighbor  (Zn,  0)a  pairs.  The  reduction  in  m  - 1 
pair  concentration  has  been  accompanied  by  an 
increase  in  concentration  of  m  >  1  pairs.  From 
the  observed  luminescence  intensities,  and  cal¬ 
culations  (following  Ref.  9)  of  the  initial  distribu¬ 
tion  of  pair  separations,  we  find  that  the  increase 
in  concentration  of  m  >  1  pairs  is  roughly  equal 
to  the  decrease  in  concentration  of  m  =  1  pairs. 

We  also  note  that  in  Fig.  1(a)  the  m  >  1  band  is 
increasing  at  the  same  rate  as  the  m  =  1  band  is 
decreasing.  Thus,  we  conclude  that  the  nearest- 
neighbor  pairs  are  dissociating  and  forming  fur¬ 
ther  separated  pairs. 

In  Fig.  2  the  discrete  DA  lines  for  second  near¬ 
est  neighbor  and  further  separated  pairs  are  re¬ 
solved.  All  except  the  m  =  2  and  m  -  3  lines  have 
previously  been  reported,7  although  we  reidentify 
a  few  of  the  lines  as  phonon-assisted  transitions. 
Replicas  of  the  m  =  2  line  are  seen  involving 
phonons  of  energies  9.2,  21.7,  27.2,  and  46.8 
meV.  The  46.8-meV  mode  is  probably  very  simi¬ 
lar  to  a  bulk  optical  (O)  mode  and  we  have  labeled 
it  as  such.  The  lower  energy  modes  are  local 
modes  (resonances)  of  the  (Zn,  O)  complex. 

These  modes  form  a  broad  band  at  -20  meV  from 
the  no-phonon  band  in  the  m  —  *»  DA  luminescence 
[Fig,  1(b) J.  Similar  localized  phonons  of  energy 
~25  meV  are  seen  for  the  isolated  neutral  O  don- 


PHOrON  ENERGY  (eV) 

HG.  l.ou -tcmivraiiiiv  •  Imioi—  ici.eplor  pair  luminescence.  The  lines  are  labeled  by  (m)  for  mth-nearest- 
neighbor  pairs  with  i  superscript  indicating  the  type  ot  phonon  involved  in  the  transition  (O’optical,  LOC  'local, 
NP  -  nu-phonuni .  Those  laliels  with  no  superscripts  refer  to  no-phonon  transitions. 
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FIG.  3.  The  dissociation  rates  of  nearest-neighbor 
(Zn,  O)  pairs.  The  rates  labeled  "recombination  en¬ 
hanced"  are  for  photoinduced  dissociation  of  the  pairs, 
in  the  limit  of  high  incident  light  intensity.  The  rates 
labeled  "thermal"  are  for  pair  dissociation  by  purely 
thermal  means. 


or.10  For  a  defect  with  small  mass,  a  localized 
low-energy  phonon  mode  implies  a  reduced  spring 
constant.  From  the  phonons  observed  here,  we 
estimate  that  the  O-Ga  spring  constant  is  about 
30%' of  the  P-Ga  spring  constant.  This  certainly 
implies  a  reduction  in  the  O-Ga  bond  strength 
compared  to  that  of  P-Ga.  A  weak  O-Ga  bond 
seems  reasonable  considering  the  small  size 
and  large  electronegativity  of  O  compared  to  P. 

The  mechanism  responsible  for  the  (Zn, O), 
dissociation  can  be  identified  by  measuring  the 
activation  energies  for  the  photoinduced  and  pure¬ 
ly  thermal  reactions. ll,“  We  have  performed 
these  measurements  and  a  complete  account  of 
them  will  be  given  elsewhere.  The  results  are 
shown  in  Fig.  3.  We  find  activation  energies  of 
0,60  *  0.07  eV  for  the  photoinduced  dissociation 
and  2.6  *0.6  eV  for  thermal  dissociation  of  the 
pairs.  The  difference  between  these  activation 
energies  is  roughly  equal  to  the  recombination 
energy  ('1.9  eV).  Thus,  we  identify  the  (Zn.O), 
dissociation  to  be  a  "recombination -enhanced" 
reaction  in  which  the  energy  available  from  elec¬ 
tron-hole  recombination  is  transferred  to  local¬ 
ized  vibrational  modes  which  leads  to  a  reaction. 
The  strong  phonon  coupling  which  is  required  for 
recombination-enhanced  reactions  is  clearly  evi¬ 
dent  in  the  photoluminescence  spectrum,  Fig. 

1(b),  The  broad  (Zn, 0)„.,  band  is  primarily 
composed  of  a  series  of  6.0-meV  phonon  repli- 
cas/,s  and  multiple  occupation  of  this  6.0-meV 


mode  probably  provides  the  “critical  reaction 
path"  for  the  dissociation  to  proceed. 

In  conclusion,  we  have  observed  photoinduced 
dissociation  of  nearest -neighbor  (Zn,  O)  pairs  in 
GaP,  and  the  formation  of  further  separated 
pairs.  As  a  result  of  the  strong  electron-phonon 
coupling  in  the  (Zn, 0)„«,  system,  it  is  possible 
that  the  dissociation  proceeds  by  the  oxygen  (or 
zinc)  ion  jumping  from  a  substitutional  to  an 
interstitial  site  and  migrating  away.  The  ion  is 
then  captured  by  a  vacancy  to  form  a  further 
separated  pair  with  the  nearest  zinc  (or  oxygen) 
ion.  From  the  local  phonon-mode  energies  it 
seems  that  the  oxygen  ion  is  weakly  bound  in  the 
GaP  lattice,  and  so  we  suggest  that  it  may  be  the 
oxygen  (rather  than  the  zinc)  ion  which  makes  the 
jump  in  the  dissociation.  The  GaP:(Zn,0)  sys¬ 
tem  has  provided  us  a  unique  opportunity  for 
studying  molecular  defect  reactions  occurring 
in  an  ordered  lattice. 
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Photoinduced  reactions  in  which  pairs  of  ZnG,  and  Op  ions  in  GaP  dissociate  and  sub¬ 
sequently  re-form  are  studied.  The  decay  of  nearest-neighbor  (Zn.O)  pair  luminescence  is 
measured  versus  laser  power  and  temperature.  The  dissociation  rate  of  the  pairs  is  de¬ 
duced  with  the  use  of  a  model  for  the  electron-hole  recombination  kinetics  which  includes 
the  effects  of  saturation  and  spatially  dependent  carrier  generation.  The  activation  energy 
for  the  photoinduced  dissociation  is  found  to  be  0.60  +0.07  eV.  Dissociation  of  the 
nearest-neighbor  pairs  by  purely  thermal  mechanisms  is  also  observed,  with  an  activation 
energy  of  2. 6+0. 6  eV.  From  these  activation  energies  and  the  observed  pre-exponential 
factors,  the  photoinduced  dissociation  is  determined  to  be  a  recombination-enhanced  defect 
reaction.  The  dissociated  nearest-neighbor  pairs  tend  to  re-form  as  pairs  with  small 
inter-ion  separation. 

I.  INTRODUCTION  observed  in  a  system  which  is  as  well  characterized 

as  GaP:(Zn.O).  Thus,  this  system  provides  a  good 
The  properties  of  gallium  phosphide  doped  with  opportunity  for  studying  the  strong  electron- 
zinc  and  oxygen  have  been  extensively  studied,  pri-  phonon  interaction  which  is  often  present  at 

manly  because  of  the  use  of  this  material  in  light-  “deep-level”  defects  in  semiconductors.  The 

emitting  diodes.  These  diodes  are  observed  to  de-  understanding  of  this  interaction  is  important  since 

grade  under  forward-bias  ope.-ation.1-5  In  a  previ-  it  significantly  affects  the  performance  of  semicon- 

ous  paper,4  we  have  shown  that  this  degradation  is  ductor  devices.5,6 

due  to  the  dissociation  of  nearest-neighbor  Photoluminescence  (optical  emission)  data  will 

(Zn.O),,.!  pairs.  When  the  GaP  is  subjected  to  be  presented  throughout  this  paper,  and  some  ex¬ 
laser  excitation,  we  observe  luminescence  from  radi-  perimental  details  are  summarized  in  Sec.  II.  In 

ative  electron-hole  recombination  at  the  (Zn,0)„  Sec.  Ill,  equations  are  written  describing  the 

pairs,  where  m  refers  to  an  mth-nearest-neighbor  electron-hole  recombination  kinetics.  We  explicitly 

pair.  At  elevated  temperatures,  the  nearest-  include  the  effects  of  saturation  of  the  (Zn.O)m 

neighbor  pair  luminescence  intensity  decreases  with  electron  occupation,  and  of  spatially  varying 

time  due  to  the  dissociation  of  the  pairs.  After  the  electron-hole  generation.  Several  limits  of  these 

nearest-neighbor  pair  dissociates,  it  re-forms  as  a  equations  that  are  applicable  to  the  experiments 

further  separated  pair.  In  this  paper,  we  will  ex-  described  here  will  be  discussed.  Simple  relations 

amine  the  kinetics  of  this  reaction.  The  word  “ki-  between  the  observed  photoluminescence  intensities 

netics”  here  has  a  twofold  meaning:  First,  we  will  and  the  defect  concentrations  are  derived.  Using 

discuss  the  electron-hole  recombination  kinetics  in  these  relations,  we  present  in  Sec.  IV  of  this  work 

order  to  relate  the  observed  photoluminescence  in-  an  analysis  of  the  decay  rates  of  the  (Zn,0)m  _  ( 

tensities  with  the  impurity  concentrations.  Second,  luminescence  intensities.  The  observed  decay  rates 

we  will  study  some  aspects  of  the  kinetics  of  the  are  shown  to  be  proportional  to  the  electron  occu- 

atomic  diffusion  process  in  which  the  (Zn.O)m  pation  of  the  pairs.  At  high  laser  powers  the  de- 

pair  dissociates  and  later  re-forms  as  a  further  cay  rates  saturate  and  equal  the  dissociation  rate  of 

separated  pair  (Zn.O)m  >i.  the  pairs.  From  the  temperature  dependence  of  the 

The  major  result  of  this  work  is  that  the  dissociation  rates,  activation  energies  for  the  reac- 

(Zn.O)m_|  dissociation  is  a  “recombination-  tions  are  deduced.  We  show  that  the  (Zn,0)m., 

enhanced"  reaction  in  which  energy  available  from  pairs  can  dissociate  by  a  recombinatioh-enhanced 

electron-hole  recombination  is  transferred  to  local-  process,  or  by  purely  thermal  mechanisms.7  In 
ized  vibrational  modes  which  leads  to  the  dissocia-  Sec.  V,  data  are  presented  on  the  distribution  of 

tion.  This  type  of  reaction  has  not  previously  been  m  >  1  pairs  that  are  formed  by  the  dissociation  of 
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the  m  =  1  pairs.  A  brief  discussion  is  given  of  the  formed  in  a  furnace  purged  with  He.  Following 

process  by  which  the  Zn  and  O  impurities  diffuse  each  anneal  the  sample  was  quenched  in  water  and 

through  the  lattice  to  form  the  (Zn,0)„  >i  pairs.  etched  for  3  min  in  3:1:1  H;S04:H202:H;0  at 

The  observed  degradation  rates  are  shown  to  be  60  *C. 

dependent  on  the  details  of  this  diffusion  process. 

The  conclusions  which  can  be  drawn  from  this 

work  are  presented  in  Sec.  VI.  III.  ELECTRON-HOLE  RECOMBINATION 

KINETICS 


II.  EXPERIMENTAL 


The  samples  used  in  this  work  were  liquid- 
phase-epitaxially  grown  p~n  junctions  on  n-type 
GaP  (single  crystal).  The  growth  technique  is 
described  by  Saul,  Armstrong,  and  Hackett.8  The 
exposed  50-jxm  thick  layer  of  the  junction  is  p 
type,  doped  with  Zn  and  O.  The  (substitutional)  O 
concentration  is  about  1016  cm-1  and  the  Zn  con¬ 
centration  is  about  3x  1017  cm-3  (based  on 
crystal-growth  conditions).  Other  impurities  in¬ 
clude  N  and  S.  The  photoluminescence  spectra 
were  obtained  using  a  liquid-N2-cooled  S-l-type 
photomultiplier  and  a  SPEX  1269  spectrometer 
with  a  1200  lines/mm  grating  blazed  at  1  fim.  For 
the  low-temperature  work  ( T  <  100  K)  reported 
here  the  sample  was  mounted  in  a  variable- 
temperature  cryostat  and  illuminated  with  the 
4545-A  line  of  an  Ar+  laser.  The  sample  was 
moved  vertically  to  permit  observation  of  different 
regions  of  the  crystal  and  the  change  in  photon 
collection  efficiency  between  these  runs  was  very 
small.  For  the  high-temperature  work  (300 
K  <  T  <  700  K)  the  sample  was  mounted  on  a 
heating  block  and  illuminated  with  the  4880-A  line 
of  an  Ar+  laser,  focused  to  a  beam  diameter  of  55 
/xm  full  width  at  half  maximum  (FWHM).  The 
maximum  sample  heating  AT  can  be  estimated 
from 


AT  = 


P 

lirka 


(1) 


where  P  is  the  absorbed  laser  power,  k  is  the  ther¬ 
mal  conductivity  [k  =0.4  W/cm  K  at  473  K  for 
GaP  (Ref.  9)],  and  a  is  the  beam  width  [FWHM 
equal  to  2crv'(21n2)].  This  formula  is  derived  by 
solving  the  heat-diffusion  equation  in  an  infinite 
slab  of  thickness  /  with  one  side  at  a  constant  tem¬ 
perature  and  the  other  side  exposed  to  a  surface 
heat  source  with  a  Gaussian  profile  in  the  limit 
a  <<  /.  For  a  typical  power  level  of  P  =  200  mW 
we  find  A T  =  35  K.  For  the  results  reported  here 
the  sample  was  in  air,  although  similar  results  were 
obtained  when  the  sample  was  in  a  He  environ¬ 
ment.  Annealing  studies  (T  >  1000  K)  were  per- 


The  kinetics  of  electron-hole  recombination  in 
GaP:(Zn.O)  have  been  extensively  studied.10-16  It 
is  not  the  purpose  of  this  work  to  obtain  a  detailed 
understanding  of  the  recombination  kinetics. 
Rather,  we  wish  to  arrive  at  approximate  values 
for  some  of  the  kinetic  parameters.  From  these 
values  we  will  determine  the  regime  of  carrier  con¬ 
centrations  in  which  our  experiments  lie,  and  this 
will  enable  us  to  make  a  reasonable  interpretation 
of  the  data  to  be  presented  in  Secs.  IV  and  V. 

We  use  a  simple  but  accurate  model  which  is 
similar  to  that  used  for  nearest-neighbor  pairs  by 
Henry,  Bachrach,  and  Schumaker.16  Each  (Zn,0)„ 
center  is  considered  to  be  a  single-electron  trap, 
with  an  associated  electron  release  rate  Rm,  capture 
rate  Cm,  and  a  total  recombination  rate  l/rm.  For 
m  —  oo,  i.e.,  isolated  oxygen,  a  second  electron  lev¬ 
el  has  been  observed.5  Here  we  will  be  considering 
centers  with  m  <  50  for  which  this  level  has  not 
been  observed  and  so  we  will  ignore  it.  For  the 
present  discussion,  the  hole  occupation  of  the 
(Zn,0)m  centers  will  be  implicitly  included  in  the 
electron  recombination  rate  l/rm.  Later,  we  will 
consider  the  dependence  of  I/rm  on  the  various 
hole  concentrations.  The  concentration  nm  of 
(Zn,0)m  centers  which  are  occupied  by  electrons  is 
given  by 

dnm  nm 

-T~  = - Rmr'm+Cm(Nm-nm)n0  , 

«f  rm 

m  =  1,2,3 .  (2) 


where  Sm  is  the  total  concentration  of  (Zn,0)m 
centers  and  n0  is  the  free-electron  concentration 
In  steady  state,  the  concentration  of  occupied 
(Zn.O)m  centers  is  found  by  equating  Eq.  (2)  to 
zero  to  yield 


»o/nm 

l+«o/"rm 


(3a) 


where 


(3b) 
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Here,  fm  represents  the  fractional  electron  occupa¬ 
tion  of  the  centers.  The  intensity  of  luminescence 
emitted  from  an  mth-nearest-neighbor  pair  is  given 
by 

=  nm  =  ^ mf  m  •  (4) 

where  l/r^*  is  the  radiative  decay  rate.  Later  in 
this  work  we  will  discuss  processes  in  which  the 
impurity  concentrations  .Vm  vary  with  time.  How¬ 
ever,  these  time  variations  are  very  slow  compared 
with  those  expressed  in  Eq.  (2);  therefore  Eq.  (3)  is 
still  valid  at  each  point  in  time. 

Let  us  consider  the  minority-carrier  (electron) 
concentration  n0.  The  variables  of  interest  in  these 
experiments  are  the  laser  power  P,  the  temperature 
T,  and  the  impurity  concentrations  IV).  The  im¬ 
purities  labeled  IV,  include  the  (Zn,0)„  centers  as 
well  as  various  intrinsic  impurities  that  could  be 
undergoing  photoinduced  reactions.3,17  In  steady 
state  n0  is  in  general  a  function  of  T,  P,  and  IV,. 
Neglecting  spatial  diffusion  of  the  carriers  we 
write 

n0  =  Gr0,  (5) 

where  G  is  the  electron-hole  generation  rate  and  r0 
is  defined  by  this  equation  to  be  the  steady-state 
minority-carrier  lifetime.  We  will  assume  that  r0 
has  no  spatial  dependence.  This  neglects  effects 
such  as  surface  recombination.  The  spatial  depen¬ 
dence  of  G  is  derived  by  considering  the  incident 
laser  beam  that  is  perpendicular  to  the  surface  ix,y 
plane)  of  the  crystal.  The  photon  density  has  a 
Gaussian  dependence  on  p  —  [(x2  +yl)]wl  and 
varies  exponentially  with  z  (the  direction  into  the 
crystal).  Thus,  the  generation  rate  can  be  ex¬ 
pressed 


G  = 


aP 

hvjlira2 


exp 


(6) 


where  P  is  the  absorbed  laser  power,  1  /a  is  the  ab¬ 
sorption  depth,  a  is  the  beam  radius,  and  hv,  is 
the  incident  photon  energy.  The  minority-carrier 
lifetime  r0  is  determined  by  all  of  the  various 
recombination  channels  in  the  material.  In  our  ex¬ 
periments  we  observe  quantities  which  are  indepen¬ 
dent  of  r0,  and  so  we  do  not  need  to  know  its  pre¬ 
cise  value.  The  minority-carrier  concentration  can 
be  estimated  by  observing  the  luminescence  intensi¬ 
ty  as  a  function  of  laser  power."  For  the  nearest- 
neighbor  pairs  at  T=  500  K  we  observe  a  signifi¬ 
cant  saturation  effects  (f  \  =0.5)  to  begin  at 


P=  100  mW  (FWHM  =  55  (im).  At  this  tempera¬ 
ture  l/r|=  10h  s_1  (Ref.  16)  <  </?|  =  10®  s_1  (Ref. 
14).  The  maximum  minority-carrier  concentration 
is  then  given  by  n0  =  nM  =  /? j  /C,  =  5x  1015  cm-3. 
This  also  equals  the  concentration  of  majority  car¬ 
riers  (holes)  introduced  by  the  laser.  From  Eqs.  (5) 
and  (6)  this  value  of  n0  corresponds  to  a  lifetime  of 
r0=  1  ns,  which  seems  quite  reasonable. 

Much  of  the  complicated  electron-hole  recom¬ 
bination  kinetics  has  been  hidden  in  the  rates  1  /r0, 

1  /rm ,  and  1  /t^.  The  rate  of  1  /rm  is  the  sum  of 
nonradiative  and  radiative  rates.  There  are  two 
nonradiative  decay  channels:  Auger  decay14  and 
multiphonon-emission  (MPE)  recombination.3 
There  are  also  two  radiative  decay  channels:  The 
bound  electron  can  recombine  with  a  hole  bound 
onto  the  same  (Zn.O)m  center,  or  with  a  hole 
bound  onto  a  distant  impurity,  namely,  a  Zn  accep¬ 
tor.10,18  The  rates  Wrm  and  1  /r^  are  therefore 
functions  of  the  temperature  and  various  hole  con¬ 
centrations.  The  concentration  of  holes  introduced 
by  the  laser  is  small  compared  with  the  free- hole 
concentration  of  p  =  1017  cm-3  for  T >  300  K. 
Furthermore,  at  these  temperatures  the  shallowly 
bound  holes  are  in  equilibrium  with  the  free 
holes.13,16  Thus,  \/rm  and  l/rJJJ11  will  be  functions 
only  of  temperature  and  will  have  very  little  spatial 
or  time  dependence. 

Now  let  us  specialize  to  the  experiments 
described  here.  For  ease  of  discussion  we  separate 
these  into  two  cases.  First  we  consider  the  high- 
temperature  results  which  will  be  presented  in  Sec. 
IV.  In  these  experiments  the  sample  is  exposed  to 
the  laser  at  some  temperature  300  K  <  T  <  700  K 
and  some  laser  power  P.  The  intensity  of 
(Zn,0)w_t  luminescence  is  observed  to  decrease 
with  time,  presumably  due  to  some  changing  im¬ 
purity  concentration  IV, .  The  luminescence  intensi¬ 
ty  is  given  by 

/'(r)=i  I  dVN^p,z,t)f,(p,z,t)  ,  (7) 

M 

where  the  spatial  dependence  of  f  ^  is  given  by  Eqs. 
(3a),  (5)  and  (6).  At  low  power  levels  the  electron 
occupation  / ,  is  proportional  to  the  minority- 
carrier  concentration  n0.  In  general,  n0  changes 
with  (V,  and  results  in  changes  in  the  (Zn,0)m  „  \ 
luminescence  intensity.  We  have  observed  these  ef¬ 
fects.  At  low  powers  a  change  in  n0  will  cause  the 
entire  luminescence  spectrum  to  change  uniformly. 
However,  at  high  powers  the  luminescence  band 
saturates;  /,  — *  1  for  most  of  the  centers  contribut- 


6332 


R.  M.  FEENSTRA  AND  T.  C.  McGILL 


25 


ing  to  the  band.  In  that  case  the  luminescence  in¬ 
tensity  is  independent  of  n0.  The  observed  de¬ 
creasing  intensity  is  due  to  a  decreasing  concentra¬ 
tion  of  (Zn,0)t  pairs,  and  the  intensity  is  directly 
proportional  to  the  pair  concentration.  The  data 
presented  in  Sec.  IV  is  in  this  high-power  regime. 

The  second  case  to  discuss  are  the  low- 
temperature  results  presented  in  Sec.  V.  Here  we 
observe  luminescence  from  different  regions  of  a 
sample  which  have  different  pair  concentrations 
Nm .  Care  was  taken  to  keep  the  temperature,  laser 
power,  and  photon  collection  efficiency  constant 
between  the  different  regions.  The  luminescence 
intensity  from  various  shallow  impurities  is  ob¬ 
served  to  be  identical  in  the  regions,  so  that  the 
electron  and  hole  concentrations  are  probably  con¬ 
stant.  This  implies  that  fm  and  1  /r^  are  con¬ 
stant.  Therefore,  from  Eq.  (4),  the  ratio  of 
luminescence  intensities  equals  the  ratio  of  pair 
concentrations  in  the  different  regions. 

IV.  REACTION  RATES 

In  Fig.  1  we  show  several  high-temperature  pho¬ 
toluminescence  spectra.  As  discussed  in  our  previ¬ 
ous  publication4  these  spectra  consist  of  a  band 
centered  at  1.7  eV  due  to  nearest-neighbor 
(Zn.O)„  pairs,  and  a  band  centered  at  1.5  eV 
due  to  further  separated  pairs  (Zn,Oi„  > ,.  There  is 
also  luminescence  at  2.1  eV  due  to  various  shallow 


FIG.  1.  High-temperature  photoluminescence  spectra 
from  (Zn.O)  donor-acceptor  pairs.  The  labels  m  denote 
luminescence  from  an  mth-nearest-neighbor  pair. 

Several  successive  scans  are  shown.  The  luminescence 
intensities  are  changing  with  time,  indicating  defect  re¬ 
actions  in  progress. 


impurities  and  possibly  free-exciton  or  band-to- 
band  recombination.  In  Fig.  1  the  time  between 
successive  scans  is  about  30  min.  We  see  that  the 
intensity  of  the  (Zn,0)m  luminescence  is  decreas¬ 
ing  and  the  intensity  of  the  (Zn,0)m  >  j  lumines¬ 
cence  is  increasing.  For  the  high  laser  power  used 
in  Fig.  1  the  m  =  1  and  m  >  1  bands  are  saturated 
(at  least  initially,  see  below),  and  so  their  changing 
intensities  are  due  to  changing  pair  concentrations. 
Clearly  the  nearest-neighbor  pairs  are  dissociating 
at  some  rate,  and  further-separated  pairs  are  form¬ 
ing  at  the  same  rate.  We  refer  to  this  reaction  as 
an  “aging”  process.  This  reaction  was  first  studied 
in  luminescence  by  Dapkus  and  Henry2  and  was 
identified  in  our  previous  paper.4 

The  (Zn,0)m  _ ,  dissociation  reaction  is  not  a 
thermal  one,  since  the  temperature  in  Fig.  1  is  only 
470  K  (sample  heating  due  to  the  laser  is  about  20 
K),  and  the  (Zn.O)m  . ,  pairs  are  thermally  stable 
up  to  300  K.  Rather,  the  reaction  is  photoin- 
duced.2  The  incident  light  creates  electron-hole 
pairs  in  the  material,  and  the  electrons  and  holes 
may  be  captured  and  may  recombine  at  the  (Zn.O) 
centers.  The  energy  to  dissociate  the  (Zn.O)m  » j 
pairs  comes  from  some  electron  (or  hole)  capture 
or  recombination  event  at  the  defect.  In  our  previ¬ 
ous  work4  we  identified  this  event  to  be  electron- 
hole  recombination.  Here,  we  will  give  a  detailed 


FIG.  2.  Decay  of  the  nearest-neighbor  (Zn,0)m 
luminescence  at  various  temperatures  T.  The  data  were 
fitted  to  a  theoretical  form  corresponding  to  the 
/(0,0)— *0  curve  shown  in  Fig.  3. 
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account  of  how  we  arrived  at  this  identification. 

The  decay  (intensity  versus  time)  of  the 
(Zn.O)m  „  |  luminescence  is  nonexponential.  Figure 
2  shows  several  decay  curves.  Typically,  the 
luminescence  intensity  initially  changes  relatively 
fast,  then  it  slows  down  to  a  long-lived  tail.  We 
attribute  this  nonexponential  behavior  to  a  nonuni¬ 
form  minority-carrier  (electron)  density  in  the  ma¬ 
terial.  This  nonuniformity  is  due  to  the  Gaussian 
profile  of  the  incident  laser  beam  and  the  exponen¬ 
tial  dropoff  of  its  intensity  into  the  sample.  A 
nonuniform  electron  density  results  in  a  nonuni¬ 
form  electron  occupation  of  the  (Zn,0)j  centers. 
Those  centers  with  high  occupation  (f  1 )  disso¬ 
ciate  at  some  rate,  and  those  centers  with  low  oc¬ 
cupation  (/,  — *0)  dissociate  at  a  much  slower  rate. 
The  observed  decay  curves  are  integrals  over  the 
sample  of  all  the  individual  exponential  decays. 

The  result  is  a  nonexponential  decay.  For  high 
laser  powers  the  observed  luminescence  is  initially 
due  almost  entirely  to  the  saturated  centers,  tn 
that  case  the  luminescence  intensity  is  proportional 
to  the  pair  concentration  and  independent  of  n0,  as 
discussed  in  Sec.  III.  At  later  times  we  observe 
those  centers  whose  occupation  is  lower.  Then,  the 
luminescence  intensity  may  depend  on  n0.  We  can 
monitor  n0  by  observing  the  intensity  of  lumines¬ 
cence  from  the  shallow  impurity  levels.  As  shown 
in  Fig.  1,  the  variation  in  nQ  occurs  over  a  much 
shorter  time  period  than  the  (Zn,0)m„|  decay  time. 
We  observe  that  this  is  true  for  P  >  100  m W  at  all 
of  the  temperatures  discussed  here.  Thus,  the  final 
part  of  the  decay  curves  are  not  affected  by  n0  be¬ 
cause  it  is  a  constant.  Therefore,  we  need  not  fur¬ 
ther  concern  ourselves  with  the  time  dependence  of 
«o  or/i- 

The  degradation  observed  by  Henry  and 
Dapkus3  is  also  reported  to  be  nonexponential.  In 
their  case  the  degradation  was  accomplished  using 
forward  biased  p-n  junctions  in  which  the 
minority-carrier  concentration  varies  exponentially 
with  distance  from  the  junction.  This  produces 
nonexponential  degradation  characteristics  similar 


to  those  reported  here.  Another  example  of  nonex¬ 
ponential  behavior  is  seen  in  carrier  capture  or 
emission  rates  as  observed  in  deep-level  transient 
spectroscopy  (DLTS). 19-20  The  carrier  capture  rate 
is  proportional  to  the  free-carrier  density.  This 
density  is  spatially  varying  at  the  edge  of  the  de¬ 
pletion  region.  The  resulting  capture  curves  are 
nonexponential  with  a  long-lived  tail  on  them.21 
However,  with  a  large  enough  reverse  bias,  the  size 
of  the  depletion  region  is  large  compared  with  the 
edge  region,  and  these  effects  can  be  minimized. 

In  our  experiments  the  laser  power  is  never  large 
enough  to  produce  an  essentially  uniform  excita¬ 
tion  region,  and  so  the  nonexponential  effects  are 
always  significant. 

Let  us  now  derive  a  functional  form  for  the  ob¬ 
served  decay  curves.  From  Eqs.  (2)— (6)  we  find 
that  the  electron  occupation  of  the  (Zn.O),  centers 
is  given  by 

2 

/(0,0)exp  — ^ — az 
la2 

la1 

(8) 

where  we  have  dropped  the  subscript  1.  We  as¬ 
sume  that  the  dissociation  rate  of  a  given  center  is 
proportional  to  its  occupation.  This  assumption  is 
justified  below.  Thus,  the  concentration  of  (Zn.O), 
centers  at  (p,z)  varies  with  time  aborting  to 

N (p,z,t)=N (0)exp[  —f(p,z)rt]  ,  (9) 

where  r  is  the  reaction  rate  and  N(0)  is  the  spatial¬ 
ly  uniform  impurity  concentration  at  t  =0.  We  are 
now  ready  to  derive  the  main  theoretical  result  of 
this  paper:  the  time  dependence  of  the  (Zn.O)t 
luminescence  intensity  in  a  model  where  the  pairs 
are  dissociating  at  a  rate  proportional  to  their  spa¬ 
tially  dependent  electron  occupation.  From  Eqs. 
(7)— (9)  the  total  intensity  of  (Zn.O)|  luminescence 
is  evaluated  to  be 


+.rr  J  0  U 


—f(0,0)urt 

1 +/(0,0)(u  -1) 

1  +/(0,0)(u  —  1 ) 


Figure  3  shows  several  examples  of  Eq.  (10)  for 
various  values  of  /" (0,0).  For  y (0,0) — ►  1  the  decay 
is  purely  exponential.  This  corresponds  to  an  in¬ 
finitely  large  excitation  region  with  uniform  occu¬ 


pation  of  all  the  centers.  For  /(0,0)-*0  the  decay 
follows  an  exponential  integral  form.  In  our  exper¬ 
iments  /(0,0)  <0.9.  From  Fig.  3  we  see  that  Eq. 
(10)  is  quite  well  approximated  by  its  /(0,0)— *0 
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limit, 

I(t)ss~[y+\a(rQt)+E\(rQt)]+b  ,  (11) 

rol 

where  a  is  a  constant,  y  is  Euler’s  constant, 
r0=/(0,0)r  is  the  observed  decay  rate,  and  £|  is 
an  exponential  integral  of  the  first  kind.  To  this 
form  we  have  added  another  constant  b  to  account 
for  the  background  intensity  due  to  (Zn.O)m  > 
Equation  (11)  provides  a  good  fit  to  all  of  the  ob¬ 
served  decay  curves.  This  is  the  real  justification 
for  using  it. 

The  data  shown  in  Fig.  2  were  fitted  to  the 
functional  form  Eq.  (11).  The  decay  curves  shown 
in  the  figure  have  had  a  background  intensity  b 
(determined  from  the  fit)  subtracted  from  them. 
The  decay  rates  determined  from  the  fits  are  plot¬ 
ted  in  Fig.  4  versus  the  incident  laser  power.  At 
low  powers  the  rates  increase  linearly  with  power; 
at  high  powers  the  rates  approach  some  limit:  the 
saturated  decay  rate.  The  decay  rates  saturate 
when  the  electron  occupation  saturates.  This  justi¬ 
fies  the  above  assumption  that  the  dissociation  rate 
is  proportional  to  the  electron  occupation.  The  sa¬ 
turation  power  level  increases  with  temperature  as 
expected  from  Eq.  (3b).  The  fact  that  the  decay 
rates  do  indeed  saturate  provides  evidence  that 
sample  heating  and  large  hole  concentrations  due 
to  the  laser  excitation  are  not  significantly  affect- 


TIME 


FIG.  3.  Theoretical  curves  for  the  decay  of  the 
(Zn,0)„_i  luminescence.  /( 0,0)  is  the  maximum  elec¬ 
tron  occupation  of  the  centers.  /(0.0)— »1  for  very  large 
laser  powers  and  /(0,0)— *0  for  small  laser  powers.  For 
the  experiments  described  here,  /iQ.O)  <0.9. 
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LASER  POWER  (mW) 

FIG.  4.  Decay  rates  of  the  (Zn,0)m  _ ,  luminescence 
vs  the  laser  power.  For  high  powers  the  decay  rates  sat¬ 
urate  due  to  saturation  of  the  electron  occupation  of  the 
centers.  The  solid  lines  are  drawn  as  a  guide  to  the  eye. 

ing  our  results. 

The  saturated  decay  rates  are  plotted  versus  tem¬ 
perature  in  Fig.  5.  We  find  an  activation  energy  of 
0.60+0.07  eV  with  a  preexponential  factor  of  ap¬ 
proximately  3x  I04  s~ Similar  activation  ener¬ 
gies  have  been  observed  in  other  GaP:(Zn.O)  degra¬ 
dation  studies.1,1  For  the  data  in  Fig.  5  we  have 

T(K) 


i/ktiev") 


FIG  5.  Dissociation  rates  of  (Zn,0)«,_i  pairs.  The 
rates  labeled  “recombination-enhanced"  are  for  photoin- 
duced  dissociation  of  the  pairs,  in  the  limit  of  high  laser 
power.  The  rates  labeled  "thermal"  are  for  pair  dissoci¬ 
ation  by  purely  thermal  means. 
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not  included  any  corrections  to  the  temperature 
due  to  sample  heating  since  these  corrections  pro¬ 
duce  a  negligible  change  to  the  activation  energy 
within  the  quoted  uncertainties.  In  addition  to 
measuring  the  activation  energy  for  the  photoin- 
duced  reaction,  we  have  also  measured  the  barrier 
for  the  nearest-neighbor  dissociation  to  prpceed  by 
purely  thermal  means.  The  procedure  consisted  of 
observing  the  absolute  intensity  of  luminescence  (at 
1 .77  eV)  from  a  sample,  putting  it  in  a  furnace  for 
some  amount  of  time,  removing,  quenching,  and 
etching  the  sample,  and  finally  observing  again  the 
luminescence  intensity.  The  annealing  tempera¬ 
tures  were  greater  than  600  *C  (the  initial  state  of 
the  crystal1)  and  so  the  luminescence  intensity  de¬ 
creased  with  annealing  time  due  to  dissociation  of 
the  nearest-neighbor  pairs.7  The  resulting  exponen¬ 
tial  decay  rates  are  shown  in  Fig.  5.  The  activa¬ 
tion  energy  is  2.6+0.6  eV  with  a  preexponential 
factor  of  approximately  10‘°  s_l.  This  activation 
energy  represents  the  barrier  for  dissociation  of 
(Zn.O)i  without  a  bound  electron  on  the  center  (let 
us  call  this  the  ground  state  of  the  complex). 

Now  let  us  consider  what  mechanisms  may  be 
responsible  for  the  photoinduced  dissociation  reac¬ 
tion.  From  the  literature6,22  we  identify  two 
mechanisms  which  have  been  observed  in  other 
systems.  First,  the  (Zn.O)i  dissociation  could  be  a 
"recombination-enhanced”  defect  reaction  (RED!.) 
in  which  energy  available  from  electron-hole 
recombination  is  transferred  to  localized  vibration¬ 
al  modes  which  leads  to  a  reaction.  This  is  the 
mechanism  we  believe  to  cause  the  reactions  ob¬ 
served  here.  The  activation  energy  of  0.60  eV  we 
identify  roughly  as  the  difference  between  the  ther¬ 
mal  barrier  (2.6+0.6  eV)  and  the  recombination  en¬ 
ergy  (  —  1.9  eV).  Since  all  of  the  recombination  en¬ 
ergy  is  used  in  surmounting  the  barrier,  the  reac¬ 
tion  is  the  result  of  nonradiative  recombination  by 
multiphonon  emission  (MPE).  Although  MPE 
recombination  has  not  been  reported  for  (Zn.O)|  it 
has  been  reported1  for  isolated  O  and  it  seems  like¬ 
ly  that  it  also  occurs  for  (Zn,0)t.  The  fraction  of 
MPE  recombination  events  which  lead  to  a  pair 
dissociation  is  estimated  to  be  10~4  (roughly,  this 
is  the  value  of  ij  as  defined  in  Ref.  6). 

The  second  mechanism  that  must  be  considered 
is  a  “charge-state”  reaction  in  which  the  (Zn,0)t 
dissociation  would  occur  thermally  from  the  excit¬ 
ed  state  of  the  complex  in  which  the  bound  elec¬ 
tron  is  present.  At  high  laser  powers  the  pairs  are 
almost  always  in  this  state.  The  observed  0.60  eV 
would  be  identified  as  the  thermal  reaction  barrier 


for  this  excited  state.  This  is  rather  different  than 
the  ground-state  barrier  of  2.6  eV.  A  nearest- 
neighbor  Coulomb  attraction  of  0.55  eV  cannot 
completely  account  for  this  difference.  Further¬ 
more,  the  observed  3x  104  s-1  is  about  9  orders  of 
magnitude  less  than  a  phonon  frequency.  This  im¬ 
plies  a  multijump  process  in  which  the  (Zn.O)  pair 
dissociates  and  recombines  109  times  before  per¬ 
manently  dissociating,  or  in  which  a  third  particle 
is  needed  for  the  reaction  and  must  diffuse  to  the 
site.  Although  we  cannot  prove  that  these  process¬ 
es  are  not  occurring,  we  consider  them  to  be  un¬ 
likely  and  we  prefer  the  simpler  explanation  of  the 
recombination-enhanced  reaction. 

V.  KINETICS  OF  ATOMIC  DIFFUSION 

In  this  section  we  will  consider  some  aspects  of 
the  (Zn.O)  pair  dissociation  and  subsequent  re¬ 
formation.  The  strong  electron-phonon  coupling 
required  for  a  REDR  is  clearly  evident  in  the 
(Zn.O)|  photoluminescence  spectra.4  This  strong 
coupling  suggests  that  the  reaction  may  proceed 
without  the  presence  of  any  additional  defects. 
Thus,  it  seems  quite  likely  that  the  oxygen  (or  zinc) 
simply  jumps  from  its  substitutional  site  to  an  in¬ 
terstitial  site  and  migrates  away.  The  mobile  ion  is 
subsequently  captured  by  a  vacancy  to  form  a 
(ZntO)m  pair  with  the  nearest  zinc  (or  oxygen)  ion. 
The  observed  activation  energy  of  0.60  eV  will  in 
general  contain  contributions  from  each  of  the  pro¬ 
cesses:  dissociation,  diffusion,  and  re-formation. 
From  Fig.  5  it  is  clear  that  the  observed  reaction  is 
a  first-order  kinetic  process.  Also,  as  shown  in 
Fig.  1  the  (Zn,0)m  >  i  pairs  form  at  exactly  the 
same  rate  as  the  (Zn.O)m  _ ,  pairs  dissociate.  Thus, 
the  rate  limiting  process  is  the  dissociation,  and  the 
observed  rates  are  indeed  the  rates  for  dissociation. 
We  emphasize  that  it  is  the  dissociation  of  the 
pairs  which  we  believe  to  be  a  photoinduced  reac¬ 
tion.  The  subsequent  diffusion  and  re-formation 
can  probably  occur  thermally  at  the  temperatures 
considered  here. 

We  have  studied  the  pair  re-formation  process 
by  measuring  the  number  of  pairs  of  a  given  shell 
number  m  which  are  produced  by  the  aging  pro¬ 
cess.  Figure  6  shows  low-temperature  spectra  in 
an  unaged  and  an  aged  sample.  The  discrete  lines 
in  the  spectra  are  (ZntO)m  donor-acceptor  (DA) 
luminescence.  The  broad  background  is  due  to  ra¬ 
diative  processes  involving  mainly  recombination 
of  electrons  bound  to  m  =  1  and  2  pairs  with  holes 
bound  to  distant  acceptors.  The  pairs  responsible 
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FIG.  6.  Low-temperature  donor-acceptor  luminescence  from  an  aged  (solid  line)  and  an  unaged  (dashed  line)  sample. 
The  lines  are  labeled  by  m  for  snth-nearest-neighbor  pairs  with  a  superscript  indicating  the  type  of  phonon  involved  in 
the  transition  (O  is  optical,  LOC  is  local,  NP  is  no-phonon).  Those  labels  with  no  superscripts  refer  to  no-phonon  tran¬ 
sitions.  Note  the  intensity  and  energy  scale  change  at  1.57  eV.  The  data  on  the  left  had  a  flat  background  subtracted 
from  it. 


for  the  discrete  DA  lines  are  labeled  in  accordance 
with  previous  identifications.4,21  As  discussed  in 
Sec.  UI,  the  ratio  of  line  intensities  is  equal  to  the 
ratio  of  pair  concentrations  in  the  samples.  In  the 
aged  material,  the  line  intensity  is  the  sum  of  that 
from  the  new  centers  produced  by  the  aging  and 
that  from  the  centers  which  existed  prior  to  the 
aging.  Thus,  the  concentration  of  centers  produced 
by  the  aging  is  given  by 

Nm=Nm  —  1  ,  (12) 

where  Im  and  I'„  are  the  luminescence  intensities 
in  the  unaged  and  aged  sample,  respectively,  and 
Nm  is  the  pair  concentration  in  the  unaged  materi¬ 
al.  Nm  can  be  computed24  by  considering  equilibri¬ 
um  between  oppositely  charged  ions  at  the  growth 
temperature  of  the  crystal.  For  small  m  this  is  a 
very  crude  model,  but  it  has  been  shown  to  give 
reasonable  results.7 

The  concentration  of  new  m  >  1  centers  pro¬ 
duced  by  the  aging  process  is  dependent  on  the 
sample  conditions  prior  to  aging.  Figure  7  shows 
typical  results,  based  on  the  integrated  DA  line  in¬ 
tensities  of  Figure  6.  We  see  that  the  dissociated 
nearest-neighbor  pairs  tend  to  re-form  as  pairs 
with  low  m.  The  distribution  of  new  centers  con¬ 
tains  some  information  about  the  pair  re-formation 
process.  At  present  we  do  not  have  a  theory  which 
allows  us  to  uniquely  identify  the  re-formation 
mechanism.  In  general,  it  appears  that  the  mobile 
ion  and/or  the  vacancy  tend  to  stay  in  the  vicinity 
of  a  dissociated  center  (i.e.,  an  isolated  Zn  or  O). 

In  any  case,  we  always  observe  that  for  1  <  m  <  10 
the  distribution  of  new  centers  increases  sharply 
with  decreasing  m.  In  fact,  we  feel  that  this  distri¬ 


bution  can  be  extended  to  include  m  =  1 .  This  ex¬ 
tension  is  reasonable  since  we  cannot  think  of  any 
physical  process  which  would  allow  the  formation 
of  m  >  1  pairs  while  prohibiting  the  formation  of 
m  =  1  pairs.  Thus,  we  conclude  that  a  significant 
number  of  the  dissociated  m  =  1  pairs  re-form  as 


FIG.  7.  The  distribution  of  new  pairs  formed  by  the 
aging  process.  N'm  is  the  concentration  of  the  mth- 
nearest-neighbor  pairs,  normalized  to  the  total  concen¬ 
tration  of  puirs  (estimated)  in  the  sample.  zm  is  the 
number  of  equivalent  mth-nearest-neighbor  sites. 
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other  m  =  1  pairs.  These  pain  are  subsequently 
dissociated  and  the  entire  aging  process  stops  when 
all  of  the  dissociated  m  =  1  pain  re-form  as  m  >  1 
pain. 

Clearly  the  observed  aging  rates  will  be  propor¬ 
tional  to  the  probability  of  not  re-forming  as  a 
m  =  1  pair.  In  order  to  minimize  the  device  degra¬ 
dation  rate,  the  probability  for  m  =  1  re-formation 
should  be  maximized.  This  probability  is  depen¬ 
dent  on  the  sample  conditions.  For  example, 

Bergh1  found  that  the  presence  of  Cu  in 
GaP:(Zn,0)  produces  an  increase  in  the  observed 
degradation  rates.  We  would  interpret  this  as  be¬ 
ing  due  to  the  formation  of  complexes  between  the 
dissociated,  mobile  ion(s)  and  the  Cu  impurity, 
thereby  eliminating  the  chance  of  producing  anoth¬ 
er  m  =  1  pair.  The  net  dissociation  rate  can  be  af¬ 
fected  by  many  other  similar  processes.  The  rates 
of  these  processes  must  be  controlled  in  order  to 
produce  devices  with  the  longest  lifetime. 

VI.  CONCLUSIONS 

In  this  paper  we  have  given  a  detailed  descrip¬ 
tion  of  photoinduced  reactions  between  nearest- 
neighbor  zinc  and  oxygen  impurities  in  gallium 
phosphide.  These  reactions  were  first  observed  by 
Henry  and  Dapkus.2,3  In  a  previous  work4  we 
demonstrated  that  the  (Zn.O),  pairs  were  dissociat¬ 
ing  and  subsequently  re-forming  as  further  separat¬ 
ed  pain.  Here,  we  report  measurements  of  the  dis¬ 
sociation  rates.  The  observed  luminescence  decay 


curves  are  analyzed  by  a  model  of  the  electron-hole 
recombination  kinetics  which  includes  the  effects 
of  saturation  and  a  spatially  dependent  generation 
rate.  From  the  resulting  dissociation  rates,  in  the 
limit  of  high  laser  power,  the  activation  energy  for 
the  photoinduced  dissociation  is  found  to  be 
0.60+0.07  eV.  We  also  report  observations  of  the 
purely  thermal  dissociation  of  the  pairs,  with  an 
activation  energy  of  2.6+0.6  eV.  From  these  ac¬ 
tivation  energies  along  with  the  observed  preex¬ 
ponential  factors,  we  conclude  that  the  dissociation 
is  a  recombination-enhanced  defect  reaction.  A 
nonradiative  electron-hole  recombination  event 
puts  the  (Zn.O),  center  into  a  highly  excited  vi- 
bronic  state,  resulting  in  the  dissociation.  The 
atomic  motion  in  the  dissociation  and  re-formation 
is  not  well  understood.  However,  we  propose  that 
in  the  dissociation  the  oxygen  (or  zinc)  ion  jumps 
to  an  interstitial  site.  This  interstitial  ion  then  dif¬ 
fuses  thermally  some  short  distance  and  is  cap¬ 
tured  by  a  vacancy  to  form  a  pair  with  the  nearest 
zinc  (or  oxygen)  ion. 
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Role  of  Fe  in  new  luminescence  lines  in  Si:Tl  and  Si:In 
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We  report  on  the  experimental  results  which  indicate  that  complexes  of  Fe  with  T1  are 
likely  to  be  the  origin  of  the  recently  observed  A,B,C  lines.  We  have  also  confirmed  that 
diffusion  of  Fe  produces  substantial  increase  in  similar  lines,  labeled  P,Q,R,  observed  in 
Si:In. 


The  luminescence  spectra  of  Si:In  have  been 
studied  in  much  detail.1-4  In  particular  long-lived 
lines  were  observed  in  the  work  of  Lightowlers  and 
Vouk2  and  by  Mitchard  et  a/.,5  who  labeled  these 
lines  P,Q,R.  Weber  et  al.  suggested  that  these 
lines  were  a  result  of  an  Fe-In  complex  with  a 
<100)  axial  symmetry.6  More  recently  a  similar 
set  of  luminescence  lines  were  observed  by  Thewalt 
et  al.  in  Si:Tl.7  In  this  Report  we  present  evidence 
that  the  luminescence  observed  by  Thewalt  et  al.  in 
Si:Tl  is  created  by  the  introduction  of  Fe  into  Tl- 
doped  Si  samples.  We  confirm  the  observations  of 
Weber  ei  al.  that  the  long-lived  lines  observed  in 
Si: In  are  also  enhanced  upon  the  introduction  of 
Fe  into  those  samples. 

Iron  is  introduced  into  the  samples  as  follows. 

A  layer  of  iron  is  evaporated  onto  the  surface  of 
the  sample  after  it  has  been  cleaned  using  a  chemi¬ 
cal  etch  (3HNOj:1CH3COOH:0.4HF).  The  iron  is 
then  diffused  into  the  sample  by  annealling  it  at 
1 100*C  for  1  h  and  then  quickly  quenching  the 
sample  in  distilled  water.  The  surface  is  then 
cleaned  again  using  the  same  chemical  etch.  We 
estimate  that  this  procedure  results  in  an  iron  con¬ 
centration  of  1015  cm-3.8  The  photoluminescence 
spectra  were  obtained  using  a  Coherent  Model  CR 
3000  K.  krypton-ion  laser  for  above-band-gap  exci¬ 
tation.  The  laser  was  operated  in  cw  mode  at  a 
peak  output  power  of  4  W.  The  luminescence  was 
then  collected  from  the  edge  of  the  sample  and 
directed  through  the  entrance  slits  of  a  Spex  Model 
1404  spectrometer  and  measured  with  a  Ge  detec¬ 
tor. 

In  Fig.  1,  we  present  the  spectrum  of  Si:Tl  be¬ 
fore  (upper  panel)  and  after  (lower  panel)  Fe  dif¬ 
fusion.  The  spectra  were  taken  at  4.2  K  and  are 
labeled  with  the  notation  used  by  Thewalt  et  al. 
The  spectrum  in  the  top  half  of  Fig.  1  is  plotted 
on  a  different  arbitrary  intensity  scale  than  that  in 
the  lower  panel.  The  large  peak  is  the  thalium 
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bound-exciton  line  and  the  smaller  peaks  at  lower 
energies  are  probably  due  to  donor-acceptor 
luminescence  resulting  from  the  unintentional  pres¬ 
ence  of  P  in  the  samples.  This  same  spectrum  is 
presented  in  the  bottom  half  of  the  figure  along 
with  the  spectrum  from  the  sample  diffused  with 
iron.  To  indicate  the  intensities  of  the  various 
lines  relative  to  the  bound-exciton  line  for  Si:Tl, 
the  two  spectra  in  this  panel  have  been  scaled  so 
that  the  T1  bound-exciton  line  is  of  the  same  inten¬ 
sity  in  both.’  These  spectra  show  that  after  the 
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FIG.  I.  (a)  Phololuminescence  spectrum  of  Si:Tl  (a) 
before  Fe  diffusion,  and  (b)  after  Fe  diffusion  at  4.2  K. 
taken  using  a  Ge  detector.  The  line  labeled  BE  is  the  T1 
bound-exciton  line. 
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diffusion  of  Fe  the  A,C,D  lines  in  Si:Tl  are  in¬ 
creased  in  intensity  from  not  being  observable  in 
the  spectra  before  diffusion  to  dominating  the 
spectrum  after  diffusion.  When  some  care  was 
taken  not  to  contaminate  the  samples  and  when  we 
subjected  them  to  the  same  cleaning  and  annealing 
procedure  without  iron  this  luminescence  did  ap¬ 
pear.  The  luminescence  was,  however,  more  than 
an  order  of  magnitude  less  intense  relative  to  the 
T1  bound-exciton  line  than  in  samples  into  which 
iron  was  diffused.  We  attribute  the  appearance  of 
the  luminescence  to  the  presence  of  iron  in  the  as- 
grown  material  or  to  slight  contamination  during 
the  processing.  This  result  is  consistent  with  the 
work  of  Mitchard  et  al?  in  Si:In  in  which  P,Q,R 
luminescence  was  observed  without  intentionally 
introducing  iron  into  their  samples. 

In  Fig.  2  we  present  the  spectrum  of  Si:In  before 
and  after  iron  diffusion  and  following  the  labeling 
scheme  of  Mitchard  et  al. 5  Again  the  spectra  were 
taken  at  4.2  K  and  are  presented  so  that  the  In 
bound-exciton  line  has  the  same  intensity  in  both. 
The  spectra  show  that  after  the  diffusion  of  Fe  the 
P,Q,R  lines  in  Si:In  are  increased  in  intensity  from 
being  barely  observable  to  dominating  the  spectra. 
We  have  also  observed  in  the  case  of  Si:In  that  the 
luminescence  can  be  diminished  by  annealing  the 
samples  at  high  temperatures  (1100*0  and  allow¬ 
ing  them  to  cool  slowly.  This  can  be  attributed  to 
the  iron  precipitating  out  of  the  samples.' 

In  summary,  we  have  shown  that  the  introduc¬ 
tion  of  Fe  in  Si:Tl  produces  the  luminescence  re¬ 
ported  by  Thewalt  et  al.,  implying  that  this  lumin¬ 
escence  is  due  to  an  Fe-Tl  complex.  We  have  also 
made  observations  which  confirm  the  conclusions 
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FIG.  2.  Photoluminescence  spectra  of  Si:In  before 
and  after  Fe  diffusion  (as  labeled)  at  4.2  K  taken  using  a 
Ge  detector.  The  line  labeled  BE  is  the  In  bound- 
exciton  line. 


of  Weber  et  al.b  that  the  P,Q,R  luminescence  ob¬ 
served  in  Si:In  is  due  to  an  Fe-In  complex. 
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Periodicity  in  the  undulation  spectra  of  GaP:N 

R.  M.  Feenstra  and  T.  C.  McGill 
California  I nsimue  of  Technology.  Pasadena.  California  91125 
(Received  22  March  1982) 

The  approximately  periodic  undulations  observed  in  luminescence  spectra  of  GaP:N  are 
shown  to  arise  from  the  long-range  order  of  the  crystal  lattice.  The  period  of  the  oscillations  is 
equal  to  the  separation  of  the  lattice  planes. 


The  luminescence  spectra  of  GaP:N  codoped  with 
acceptors  displays  approximately  periodic  undula¬ 
tions.  '  A  variety  of  explanations  have  been  proposed 
for  this  undulatory  behavior. J-s  Finally,  it  was  shown 
by  Street  and  Wiesner3  that  the  luminescence  is  due 
to  a  series  of  zero-phonon  lines  of  excitons  bound  to 
N-acceptor  pairs.  The  undulatory  behavior  of  the 
spectra  results  from  the  properties  of  Nir),  the  dis¬ 
tribution  of  acceptor  sues  located  at  distance  r  from 
the  donor  site.  Street  and  Wiesner  observed  that 
fluctuations  in  Nir)  existed,  and  they  speculated 
that  these  fluctuations  were  purely  random  in  origin. 
In  this  paper,  we  show  that  Nir)  is  composed  of  a 
series  of  functions  which  are  approximately  periodic 
in  r.  The  period  of  oscillations  is  given  by  2ir/|S|. 
where  G  is  a  reciprocal-lattice  vector  This  period 
equals  the  perpendicular  separation  of  those  lattice 
planes  labeled  by  G.  This  leads  us  to  the  following 
geometrical  interpretation  for  the  periodic  behavior: 
Consider  a  sphere  of  radius  r  centered  at  the  donor 
site.  As  r  increases  the  sphere  intersects  new  lattice 
planes.  The  total  number  of  points  Hr)  enclosed  by 
the  sphere  will  show  a  sharp  increase  each  lime  a  new 
plane  is  intersected,  and  N (r)-dL  ir  )/dr  will  show  a 
corresponding  peak.  Thus,  Nir)  will  show  oscillatory 
behavior  with  period  given  by  the  lattice  plane  spac¬ 
ing.  In  particular,  the  low-frequency  oscillations  in 
N  ir)  result  from  lattice  planes  of  high  symmetry 
1(100),  (110),  etc  ] 

It  is  straightforward  to  numerically  compute  Nir), 
where  Nir)  dr  is  the  number  of  acceptor  sites  located 
a  distance  between  r  and  r  +  dr  from  the  donor  site. 

At  zero  temperature,  Nir)  will  be  a  series  of  5  func¬ 
tions,  one  for  each  donor-to-acceptor  shell.  Experi¬ 
mentally,  these  discrete  lines  are  broadened.  To 
simulate  this  line  broadening  we  replace  each  S  func¬ 
tion  by  a  Gaussian  of  width  c r.  The  averaged  num¬ 
ber  of  sites  Nir)  is  then  computed  by  summing  each 
Gaussian  weighted  by  the  number  of  sites  in  each 
shell.  Figure  1  shows  1 V (r)  for  a  type-ll  spectrum  in 
a  zinc-blende  lattice  (donor  and  acceptor  ions  on  dif¬ 
ferent  sublattices).  For  a  broadening  of  <7-0  3  A 
the  results  are  basically  identical  with  those  of  Street 
and  Wiesner.’  The  approximately  periodic  behavior 
of  Nir)  is  clearly  seen.  For  <r-0.8  A  the  period  of 
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the  oscillations  as  measured  from  Fig.  1  is  2.73  A 
-a o/2  where  a0  is  the  lattice  constant.  This  period 
will  be  theoretically  justified  below. 

Let  us  now  derive  an  analytic  form  for  Nir). 
Consider  a  lattice  composed  of  all  possible  acceptor 
sites.  For  simplicity,  assume  that  this  lattice  has  only 
one  atom  per  unit  cell.  The  acceptor  sites  are  then 
specified  by 

T  —  /»|3*i  +/i2a2  +  njaj  ,  (1) 

where  a  t,  a2,  and  a  3  are  the  primitive  translation 
vectors  and  m,  nj,  and  /i3  run  over  all  integer  values. 
Let  Vc  denote  the  volume  of  the  unit  cell,  and  let  the 
reciprocal-lattice  vectors  be  labeled  by  G.  Let  the 
donor  position  be  specified  by  R*  The  total  number  • 
of  lattice  points  enclosed  by  a  sphere  of  radius  r,  cen¬ 
tered  at  1?*,  can  be  expressed  as 

L(r.R,)-2e(r-|T-R,|)  ,  (2) 

T 


FIG.  1.  Broadened  distribution  Nir)  of  acceptor  sites  lo¬ 
cated  a  distance  r  from  donor  site  for  (N.  acceptor)  pairs  in 
GaP.  Upper  curve:  <r  -0.3  A  broadening  with  the  peaks  la¬ 
beled  as  in  Ref  5.  Lower  curve:  <r— 0.8  A  broadening 
showing  the  lowest  frequency  oscillations  with  period  2.73 
A  —  a0/2. 
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where_9  is  a  unit  step  function.  The  function 
L(r. R„)  is  periodic  in  and  so  it  can  be  expanded 
in  a  Fourier  series, 

L(r,  Rrf)  —  ^ciGlexpOG  Rrf)  .  (3) 

G 

Following  Kendall  and  Rankin,6  the  coefficients 
c(G)  are  given  by 

c(G)--£-  f  i.(r.R,)exp(-iGRaWR, 

V e  *  ceil 

-  —  ^  J"ct|i0(/--|f-Rrf|)exp(-/G-Rrf)  dKj 

-jrirx.*' p<jg  K)  tfK 


-kwMt5W) 


where  j\  is  a  spherical  Bessel  function  of  the  first 
kind.  The  distribution  of  acceptor  sites  located  a  dis¬ 
tance  r  from  the  donor  can  now  be  evaluated  to  be 


N(r) 


dL  ( r,  Rrf) 
dr 


(4) 


It  is  clear  from  Eq.  (4)  that  N(r)  is  approximately 
periodic  in  /-with  frequencies  given  by  151/2 rr. 
Depending  on  the  value  of  R*  some  frequencies  may 
not  be  observed.  In  particular,  if  cos(G  -  R*)  -0 
then  the  frequency  151/27 t  is  forbidden.  The  above 
discussion  can  easily  be  generalized  to  include  more 
than  one  type  of  acceptor  site  (i.e.,  more  than  one 
atom  per  unit  cell  in  the  lattice).  N{r)  then  contains 
a  summation  of  the  form  Eq.  (4)  for  each  inequiv¬ 
alent  type  of  acceptor  site.  Inequivalent  donor  sites 
can  be  handled  in  a  similar  manner.  Finally,  we  note 
that  the  role  of  acceptors  and  donors  in  the  above 
discussion  can  be  interchanged  and  the  results  are  the 
same. 

These  results  can  be  applied  to  the  problem  of 


GaP:(N,  acceptor).  The  acceptor  ions  are  located  on 
a  face-centered-cubic  lattice  with  Vc  -  /4  where 

do —  5.45  A  is  the  lattice  constant.  The  frequencies 
of  oscillation  are  found  to  be  -/J/ao,  2 /a0, 

2%/2/ao .  The  donor  is  located  at  Rrf 

-a0(  7.7.7)  so  that  the  frequency  t/3/oo  is  forbid¬ 
den.  Therefore,  the  lowest  frequency  will  be  2/a0  as 
observed  in  the  numerical  calculation  shown  in  Fig.  1 
for  <7—0.8  A.  For  a  broadening  of  <7—  0.3  A  some 
of  the  higher-frequency  components  are  observed. 

Now  that  we  have  a  mathematical  reason  for  the 
observed  undulations,  it  is  interesting  to  seek  some 
simple  geometrical  interpretation  of  the  results. 

From  the  form  of  Eq.  (4)  this  interpretation  is  obvi¬ 
ous.  The  reciprocal-lattice  vectors  label  lattice  plane 
in  real  space.  The  period  of  oscillation  2tr/|Gl  is 
simply  equal  to  the  perpendicular  distance  between 
lattice  planes.  Consider  the  number  of  lattice  points 
contained  in  a  sphere  of  radius  r.  As  r  increases,  the 
sphere  will  intersect  new  lattice  planes.  Each  time 
such  an  intersection  occurs  there  will  be  a  sharp  in¬ 
crease  in  L(r)  and  a  corresponding  peak  in  N(r). 
Thus,  the  frequency  |Gl/2ir  is  associated  with  those 
lattice  planes  labeled  by  G.  As  a  corollary,  we  note 
that  those  acceptor  ions  which  contribute  to  the 
IO|/2ir  fluctuations  in  A f(r)  will  be  concentrated 
around  thejlirections  from  the  donor  which  are 
parallel  to  G. 

In  conclusion  we  have  shown  that  the  apparent 
periodicity  of  GaP:(N,  acceptor)  luminescence  spectra 
is  the  result  of  long-range  order  in  the  crystal  lattice. 
In  principal  the  crystal  symmetry  could  be  deduced 
from  the  spectra,  but  in  practice  this  may  be  difficult 
since  it  would  require  detailed  knowledge  of  the  rela¬ 
tionship  between  luminescence  energy  and  the 
donor-acceptor  separation. 
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Photoluminescence  of  Si-rich  Si-Ge  alloys 

G.  S.  Mitchurd  and  T.  C.  McGill 
California  Institute  of  Technology.  Pasadena.  California  91125 
(Received  27  July  1981) 

A  detailed  study  of  the  photoluminescence  spectrum  of  Si'.,  Ge,  is  presented  for 
x  =s0. 1.  Undoped  and  In-doped  alloy  samples  are  investigated.  Spectral  features  are 
identified,  and  the  effects  of  the  disordered  nature  of  the  alloy  on  the  luminescence  spec¬ 
trum  are  discussed.  In  particular,  mechanisms  for  the  observed  broadening  of  bound- 
exciton  luminescence  are  discussed,  and  the  shift  in  alloy  luminescence  energy  is  obtained. 

At  low  temperatures  a  luminescence  feature  is  observed  which  may  result  from  recombi¬ 
nation  of  excitons  bound  to  local  fluctuations  in  alloy  composition. 


I.  INTRODUCTION  mechanisms  for  the  observed  bound-exciton 

luminescence  broadening  are  discussed.  Comparis- 
The  properties  of  alloy  semiconductors  have  on  of  the  luminescence  energies  from  Si  and  the  al- 

been  of  general  interest  for  some  time.  Such  alloys  loy  samples  results  in  values  for  the  band-gap  shift 
provide  a  convenient  system  with  which  to  experi-  and  free-exciton  binding  energy  for  the  values  of 

mentally  and  theoretically  study  the  effects  of  the  composition  parameter  x  which  we  were  able 

disorder,  which  can  be  varied  with  alloy  composi-  to  study.  Finally,  some  evidence  is  presented 

tion.  The  band  gap  can  also  be  varied  with  alloy  which  suggests  that  we  observe  at  low  tempera- 

composition,  and  for  this  reason  considerable  effort  tures  luminescence  which  is  the  result  of  recombi- 

has  been  directed  towards  the  development  of  in  nation  of  excitons  bound  to  local  fluctuations  in  al- 

trinsic  and  extrinsic  alloy  optoelectronic  devices.  loy  composition. 

The  study  of  the  luminescence  properties  of  alloy  The  remainder  of  this  paper  is  organized  in  the 

semiconductors  can  provide  useful  and  relatively  following  way.  In  Sec.  II,  we  discuss  the  Sii_xGex 

easily  interpreted  information  regarding  the  pro-  alloy  samples  studied,  and  the  general  features  of 

perties  of  the  alloys  and  the  consequences  of  their  the  experimental  technique.  In  Sec.  Ill,  the  results 

disordered  nature.  Most  attention  has  been  direct-  of  measurements  on  the  undoped  Sii_xGex  sam- 

ed  towards  the  III-V  ternary  alloys,  where  lumines-  pies  are  presented  and  discussed.  In  Sec.  IV,  we 

cence  processes  have  been  studied  extensively.1-3  consider  the  results  obtained  for  one  In-doped 

Luminescence  of  II-VI  ternary  alloys  has  also  been  Sii_xGex  sample.  In  Sec.  V,  we  discuss  the  gen- 

studied,  most  notably  Hg!_xCdxTe  (Ref.  4)  where  eral  mechanisms  for  luminescence  line  broadening 

ir  detector  applications  are  particularly  important.  in  the  alloy,  compare  the  alloy  luminescence  ener- 

Another  alloy  material  of  interest  is  the  binary  gies  with  those  observed  in  Si,  and  consider  the 

alloy  Si|_xGex.  The  luminescence  properties  of  evidence  for  excitons  bound  to  composition  fluc- 

Si(_xGex  are  not  particularly  well  known,  al-  tuations.  Finally,  in  Sec.  VI,  the  results  are  sum- 

though  free-  and  bound-exciton,5  donor-acceptor,6  marized. 

and  electron-hole  droplet7  recombination  has  been 
observed  in  Ge-rich  alloys.  In  this  paper  we  report 

the  first  detailed  measurements  of  luminescence  II.  EXPERIMENTAL  CONSIDERATIONS 

from  Si-nch  alloys,  in  particular  alloys  for  which 

x  =0. 1.  Intrinsic  and  impurity-related  lumines-  A.  Si-Ge  alloy  samples 

cence  is  discussed,  as  well  as  certain  properties  of 

the  luminescence  which  arise  as  a  result  of  the  In  this  paper,  we  present  measurements  of  the 

compositional  disorder  of  the  alloy.  In  particular,  photoluminescence  properties  of  three  Si,_xGex 

we  have  been  able  to  identify  the  luminescence  samples,  two  of  which  were  not  intentionally  doped 

features  which  result  from  free-exciton  and  and  the  third  of  which  was  doped  with  In.  The 

bound-exciton  recombination  in  the  alloy.  Also,  crystals  were  grown  by  the  Czochralski  technique. 
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Impurity  concentrations  and  alloy  compositions 
were  determined  at  Hughes  Research  Laboratories, 
where  the  crystals  were  grown.  Impurity  concen¬ 
trations  were  established  on  the  basis  of  Hall-effect 
measurements.  These  measurements  show  that 
residual  concentrations  of  B  and  P  impurities  are 
present  in  all  the  samples.  Alloy  compositions 
were  obtained  from  a  variety  of  techniques;  results 
of  the  electron  microprobe,  density,  and  x-ray  dif¬ 
fraction  measurements  were  in  excellent  agreement. 
The  samples  studied,  their  impurity  concentrations, 
and  their  alloy  compositions  are  given  in  Table  I. 


B.  Experimental  technique 

The  Si-Ge  alloy  samples  described  above  were 
mechanically  etched,  and  then  mounted  in  a  Janis 
variable-temperature  cryostat.  The  luminescence 
resulted  from  above  band-gap  optical  excitation  by 
a  Spectra- Physics  Model  166  Ar4"  laser,  which 
could  be  operated  in  pulsed  or  cw  mode.  The 
luminescence  was  collected  from  the  edge  of  the 
sample,  wavelength-analyzed  with  a  Spex  Model 
1269  spectrometer,  and  detected  with  an  RCA 
7i02  S-l  photomultiplier  cooled  to  liquid-nitrogen 
temperature.  The  photomultiplier  output  was  pro¬ 
cessed  with  standard  photon-counting  electronics, 
acquired  on  a  Nuclear  Data  ND-60  multichannel 
analyzer,  and  finally  passed  to  a  PDP  1 1/34  com¬ 
puter  for  analysis  and  storage. 

III.  EXPERIMENTAL  RESULTS  AND 
DISCUSSION  FOR  UNDOPED  Si,_,Ge, 

In  this  section,  the  results  of  the  application  of 
the  photoluminescence  technique  to  the  undoped 
samples  are  considered.  Two  undcped  Sii_,Ge, 
samples  were  available,  C077  (x  =0  11)  and  C021-3 
(x  =0.067).  As  indicated  in  Table  I,  these  samples 


have  low  residaal  concentrations  of  B  and  P  im¬ 
purities,  in  spite  of  not  being  intentionally  doped. 

A.  Typical  photoluminescence  spectrum 
and  identification  of  phonon  replicas 

A  typical  low-temperature  photoluminescence 
spectrum  of  sample  C077  is  shown  in  Fig.  1.  Two 
broad  features  are  visible  at  the  high-energy  end  of 
the  spectrum  and  are  labeled  FE  (free  exciton)  and 
BEp  (bound  exciton).  Replicas  of  these  lines,  la¬ 
beled  FE(TO)  and  BEP(TO),  are  observed  approxi¬ 
mately  58  meV  lower  in  energy.  Since  this  is  the 
transverse-optical  (TO)  phonon  energy  in  Si,  the 
low-energy  lines  are  interpreted  as  being  due  to  TO 
phonon  replicas  of  the  higher-energy  lines. 

B.  Identification  of  free-exciton  luminescence 

In  this  section  we  consider  the  effect  of  increas¬ 
ing  the  sample  temperature  on  the  photolumines¬ 
cence  spectrum  of  sample  C077.  In  Fig.  2  we  see 
that  as  the  temperature  is  increased  the  line  labeled 
BEP  thermalizes  with  respect  to  the  line  labeled 
FE.  Above  10  K.  (Fig.  3),  the  line  labeled  BEp  is 
no  longer  visible,  and  the  line  labeled  FE  assumes 
a  shape  characteristic  of  free-exciton  (FE)  recom¬ 
bination  in  Si.  In  intrinsic  Si,  FE  recombination 
luminescence  has  been  fitted  very  accurately  with  a 
line  shape,1 

HE)ccV'E-E0exp  ~~  ,  (1) 

kBT 

where  HE)  is  the  luminescence  intensity  at  photon 
energy  E,  £0  is  the  FE  threshold  energy,  and  T  is 
the  temperature.  In  this  expression  wc  assume  par¬ 
abolic  FE  bands  and  a  Boltzmann  distribution  of 
FE  center-of-mass  kinetic  energies.  The  high- 
temperature  luminescence  from  sample  C077  can 


TABLE  I.  Si|_,Ge,  alloys  studied  using  photoiuminescence.  The  samples  were  grown 
and  characterized  at  Hughes  Research  Laboratories. 


x 


Sample 

(cm-’) 

AT  p 

(cm'1] 

(cm-5) 

Electron 

microprobe 

Density 

X-ray 

diffraction 

C077 

J.ft  <  ID" 

2.x  v  to” 

0.1115 

0. 1 1 3 

0.113 

0)21-1 

2. 1  "*  l()'4 

4  1  s  I014 

0.0677 

0.0683 

C0V3 

5.0  x  10” 

5..1  x  It)14 

2.5  x  to1" 

0  104 
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FIG.  I.  Typical  photoluminescence  spectrum  of  un¬ 
doped  Si  i  _,Ge,  sample  C077.  Refer  to  the  text  for  an 
explanation  of  the  line  assignments. 


FIG.  2.  Temperature-resolved  spectra  of  undoped 
Si|_,Ge,  sample  C077  from  4.2  to  8.9  K.  The  scale 
factors  (e.g.,  X2)  give  the  relative  intensity  magnifica¬ 
tion.  Note  that  thermalization  of  the  line  labeled  BEp  is 
observed.  In  addition,  the  BEp  peak  position  shifts  to 
lower  energy  as  the  temperature  is  increased.  At  high 
temperatures,  the  luminescence  of  the  line  labeled  FE 
assumes  a  shape  characteristic  of  free-exciton  lumines¬ 
cence  in  Si. 


FIG.  3.  Temperature-resolved  spectra  of  undoped 
Sii_xGe.  sample  C077  from  10  to  25  K.  'The  scale  fac¬ 
tors  (e.g.,  X  2)  give  the  relative  intensity  magnification. 
Note  that,  at  these  temperatures,  the  luminescence  of 
the  line  labeled  FE  assumes  a  shape  characteristic  of 
free-exciton  luminescence  in  Si. 


also  be  fitted  very  well  with  the  line  shape 
described  by  Eq.  (1).  An  example  of  this  fit  is 
shown  in  Fig.  4.  In  all  cases,  the  fit  temperature 
obtained  in  this  manner  was  within  1  K  of  the 
measured  bath  temperature.  Also,  the  threshold 
energy  remained  constant  within  0.05  meV.  On 
the  basis  of  this  analysis,  the  line  labeled  FE  is 
identified  as  resulting  from  no-phonon  (NP)  FE 
recombination.  Note  that  this  intrinsic  NP 
luminescence  is  greatly  enhanced  in  the  alloy,  since 
the  Ge  atoms  can  act  as  momentum-conserving 
scattering  centers.5 


C.  Identification  of  bound-exciton  luminescence 

The  identification  of  the  FE  line  obtained  in  the 
preceding  section,  splitting  between  the  FE  thres¬ 
hold  and  the  peak  position  of  the  line  labeled  BEP, 
and  the  thermal  behavior  shown  in  Fig.  2,  suggest 
that  the  line  labeled  BEp  is  due  to  the  NP  recom- 
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FIG.  4.  Frec-exciton  luminescence  from  undoped 
Si,_xGe,  sample  C077  and  least-squares  fit  of  the 
theoretical  line  shape  [Eq.  (1)].  The  sample  temperature 
obtained  as  a  result  of  the  fit  Tr„  is  shown  with  the 
measured  bath  temperature 

bination  of  excitons  bound  by  about  4  meV  to  a 
shallow  level.  The  pump-power  dependence  at  low 
power  shown  in  Fig.  5  supports  this  proposal.  As 
this  figure  demonstrates,  the  FE  to  BEP  lumines¬ 
cence  intensity  ratio  is  independent  of  pump 
power. 

Candidates  for  such  a  level  are  clearly  B  and  P, 
which  have  binding  energies  in  the  neighborhood 
of  4  meV  in  Si  (3.9  and  4.7  meV,  respectively*), 
and  which  are  the  most  common  shallow  impuri¬ 
ties  in  the  undoped  material.  Photoluminescence 
measurements  of  Si  containing  approximately 
equal  background  concentrations  of  B  and  P  show 
that  P  bound-ex citon  (BEp)  luminescence  is  more 
intense  than  B  bound-exciton  (BEb)  luminescence 
in  the  NP  region  by  at  least  an  order  of  magnitude. 
It  seems  reasonable  to  conclude,  therefore,  that  the 
line  labeled  BE  in  the  luminescence  from  sample 
C077  is  primarily  due  to  NP  BEp  recombination. 

The  high-power  pump-power  dependence 
presented  in  Fig.  6  supports  this  conclusion.  In 
undoped  Si,  low-energy  lines  due  to  bound- 
multiexciton  complexes  (BMEC)  appear  as  the 
pump  power  is  increased.  The  splitting  between 
the  BEp  line  and  the  first  BMEC  line  is  2.2  meV 
for  B  and  3.6  meV  for  P.11  A»  iliuwti  in  Fig.  6,  in 
the  luminescence  spectrum  of  sample  C077  the 


FIG.  5.  Pump-power  dependence  of  the  luminescence 
from  undoped  Si|_,Ge,  sample  C077  at  low  pump 
powers.  The  scale  factors  (e.g.,  X  6)  give  the  relative  in¬ 
tensity  magnification.  Note  that  the  BEr  to  FE 
luminescence  intensity  ratio  is  independent  of  pump 
power  at  these  power  levels. 


BEp  line  develops  a  low-energy  shoulder  as  the 
pump  power  is  increased,  which  resolves  into  a 
separate  line  at  high  pump  powers.  The  separation 
between  this  line  and  the  BEP  line  is  about  3.8 
meV,  which  is  consistent  with  the  splitting  between 
the  P  BE  and  first  BMEC  in  Si.  Since  no  line  is 
observed  which  would  correspond  to  B  BMEC 
luminescence,  our  conclusion  regarding  the  origin 
of  the  BEP  line  in  the  alloy  seems  justified. 

D.  Comparison  of  pbotoiutninesccnce  spectra 

The  interpretation  of  the  photoluminescence 
spectrum  of  sample  C077  developed  in  the  previ¬ 
ous  sections  is  consistent  with  our  investigation  of 
the  x  =0.067  sample,  C021-3.  In  Fig.  7,  the 
luminescence  spectra  of  samples  C077  and  C021-3 
are  compared.  We  see  that  the  ratio  of  BEp  to  FE 
luminescence  is  considerably  greater  for  sample 
C02I-3.  This  is  consistent  with  Table  I,  which  in- 
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COMPARISON  OF  LUMINESCENCE 
FROM  SAMPLES  C077  AND  C02I-3 


sanviie  CO  7  7 
i  =3  O.l  i 


'born  *- 

pomp  power  s  iQmvV 


5*  resolution 


sample  C077 

A  0.1  I 

Na  =  3.6  xiO'Jcm3 
N--  2.8  *  '0‘s  cm 


sample  CO 2i- 

3 

1  5 

0.067 

nb 

-  2.3*10'* 

cm'3 

‘  4.1  X  IQ1* 

cm*1 

BEJ  TO) 

ENERGY  (meV) 

FIG-  6.  Pump-power  dependence  of  luminescence 
from  undoped  Si,_,Ge„  sample  C077  at  high  pump 
powers.  The  scale  factors  le.g.,  X7)  give  the  relative  in¬ 
tensity  magnification.  Note  the  new  feature  which  ap¬ 
pears  on  the  low-energy  side  of  the  BEp  line  at  high 
pump  powers.  This  feature  is  discussed  in  the  text. 


dicates  that  impurity  concentrations  in  sample 
C021-3  are  approximately  an  order  of  magnitude 
higher.  In  addition,  we  see  that  the  luminescence 
spectrum  of  sample  C021-3  is  shifted  to  higher  en¬ 
ergy.  This  feature  will  be  discussed  in  detail 
below. 


IV.  EXPERIMENTAL  RESULTS  AND 
DISCUSSION  FOR  ln-DOPED  Si,_,Ge, 


■OTO  MOO  1110  1120  M30  1140 

ENERGY  (mev) 

FIG.  7.  Comparison  of  the  luminescence  spectra  of 
samples  C077  and  C021-3.  The  scale  factors  (e.g., 

X  1.5)  give  the  relative  intensity  magnification.  Note 
that  BEp  luminescence  is  much  more  intense  from  sam¬ 
ple  C021-3  than  from  sample  C077.  Also,  note  that  the 
luminescence  from  sample  C021-3  is  at  higher  energy 
than  the  luminescence  from  sample  C077.  These 
features  are  discussed  in  the  text. 


comes  visible  only  at  high  temperatures.  At  some¬ 
what  lower  energies  a  broad  feature  labeled  L  is 
visible,  which  moves  to  lower  energy  as  the  tem¬ 
perature  is  increased.  The  L  luminescence 
develops  a  low-energy  shoulder  at  high  tempera¬ 
tures,  which  has  been  labeled  FE(TO). 


The  subject  of  the  investigations  reported  in  this 
section  was  the  In-doped  Si,_xGex  sample  C093, 
for  which  x  =0. 10. 


A.  Typical  photoluminescence  spectra 

Typical  luminescence  spectra  obtained  from 
Si|_xGe, :  In  sample  COO  3  are  shown  in  Fig.  8. 
Spcclru  measured  at  vuiuni*  Iciiipcruluicx  arc 

shown,  t  wo  fculuics  uic  (uommciil  m  the  spcctiu. 
The  higher-energy  luminescence,  labeled  FE,  be- 


B.  Identification  of  free-exciton  luminescence 

To  begin  with,  we  note  that  the  luminescence  la¬ 
beled  FE  has  a  line  shape  and  temperature  depen¬ 
dence  characteristic  of  NP  FE  recombination  in  Si, 
as  described  in  detail  in  Sec.  III.  Equation  (1), 
which  describes  the  NP  FE  line  shapes,  can  be  fit¬ 
ted  to  the  FE  luminescence  from  sample  C093  very 
accurately.  The  FE  threshold  energy  obtained 
from  such  fits,  1 1 15.8  mcV,  remained  constant 

wiil  iu  utxiui  0.3  mcV.  However,  the  fit  tempera¬ 
tures  Tn,  were  systematically  higher  than  the  mea- 
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FIG.  8.  Temperature-resolved  spectra  of  Si,_,Ge,  In 
sample  C093  from  4.2  to  30  K.  The  scale  factors  (e.g., 

X  2)  give  the  relative  intensity  magnification.  Refer  to 
the  text  for  a  discussion  of  the  line  assignments  shown 
in  the  figure.  Where  FE  luminescence  is  observed,  the 
sample  temperature  obtained  from  a  least-squares  fit  to 
the  FE  line  shape,  Tr,„  is  shown  with  the  measured  bath 
temperature  T*, h.  The  discrepancy  between  these  re¬ 
sults  is  discussed  in  the  text. 

sured  bath  temperatures  T^,k.  These  temperatures 
are  indicated  in  Ftg.  8.  It  seems,  therefore,  that  ei¬ 
ther  the  FE  luminescence  is  subject  to  some 
broadening  as  a  result  of  the  disordered  nature  of 
the  alloy,  or  the  sample  is  being  locally  heated  by 
the  laser  pump.  Of  these  two  possibilities,  the 
second  seems  most  likely.  First,  one  would  expect 
some  broadening  mechanism  to  affect  the  FE 
threshold  behavior  as  well  as  generally  broaden  the 
line.  That  is,  we  would  expect  to  see  the  FE 
threshold  “smear  out”  and  become  less  abrupt. 
However,  this  behavior  is  not  observed.  As  previ¬ 
ously  mentioned,  the  line  shape  is  described  by  Eq. 
(1)  very  well.  Second,  as  indicated  in  Table  I  sam¬ 
ple  C093  has  impurity  concentrations  which  are 
about  3  orders  of  magnitude  higher  than  those  in 
sample  C077.  As  a  result,  much  higher  pump 
powers  are  needed  to  observe  FE  luminescence  at 
all  in  sample  C093.  At  these  laser  intensities,  sam¬ 


ple  heating  must  be  considered  a  very  likely  cause 
of  the  observed  FE  broudening. 

If  this  explanation  is  accepted,  the  luminescence 
labeled  FE  can  be  interpreted  as  resulting  from 
NP-FE  recombination  in  the  Si|_j,Ge„:  In  alloy. 
Since  the  feature  labeled  FE(TO)  is  separated  from 
the  FE  peak  by  approximately  the  Si  TO  phonon 
energy,  we  interpret  this  shoulder  on  L  as  being 
due  to  TO  phonon-assisted  FE  recombination. 

C.  Identification  of  bound-exciton 
luminescence 

Our  interpretation  of  the  luminescence  labeled  L 
is  necessarily  considerably  more  complicated  than 
that  of  the  FE  luminescence.  There  are  a  number 
of  processes  which  are  undoubtedly  contributing  to 
the  L  line.  Primarily,  BEb,  BEp,  and  BE)n 
luminescence  must  be  considered.  Indium,  of 
course,  is  the  majority  dopant,  and  will  dominate 
the  NP-BE  luminescence.  The  B  concentration  is 
about  an  order  of  magnitude  greater  than  the  P 
concentration.  It  follows  from  our  discussion  of 
BE  luminescence  from  undoped  Si!_xGex  that  B 
and  P  will  probably  contribute  equally  to  the  NP- 
BE  luminescence.  In  addition,  the  L  luminescence 
is  sufficiently  broad  that  transverse-acoustic  (TA) 
phonon  replicas  of  the  FE  and  BE  lines  may  con¬ 
tribute  as  well.  It  is  unlikely  that  NP-BMEC  or 
electron-hole— droplet  (EHD)  luminescence  will 
have  any  significant  influence  on  the  spectra 
presented  here.  The  impurity  concentrations  in 
sample  C093  preclude  the  observation  of  such  ef¬ 
fects  at  the  pump  powers  utilized  in  this  study. 

I.  Time-resolved  spectra  and  isolation 
of  B  bound-exciton  component 

In  an  effort  to  separate  the  BE  components 
which  produce  the  L  luminescence,  the  time  evolu¬ 
tion  of  the  luminescence  from  sample  C093  was 
measured  at  low  pump  power.  The  time-resolved 
spectra  are  shown  in  Fig.  9  for  low-power  laser  ex¬ 
citation.  As  this  figure  demonstrates,  the  L 
luminescence  is  composed  of  a  relatively  sharp 
high-energy  component  with  a  long  lifetime,  and  a 
broad  low-energy  component  with  a  short  lifetime. 
Consideration  of  the  BEb,  BEP  and  BE|„  decay 
times  in  Si  (1055,  272,  and  2.7  nsec,  respectively10) 
leads  to  the  suggestion  that  the  sharp,  long-lifetime 
component  is  primarily  a  result  of  BEb  recombina¬ 
tion.  The  broad,  short-lifetime,  low-energy  should¬ 
er  is  therefore  probably  due  to  BE[n  recombination. 
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FIG.  9.  Time-resolved  spectra  of  Si|_,Ge,In  sample 
COR 3.  The  scale  factors  (e.g.,  X5)  give  the  relative  in¬ 
tensity  magnification.  At  long  times,  a  sharp,  long- 
lifetime  component  of  L  luminescence  is  isolated. 

2.  Temperature-resolved  spectra  and  isolation 
of  In  bound-exciton  component 
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FIG.  10.  Gated,  temperature-resolved  spectra  of 
Sii-,Ge,:In  sample  C093  at  high  pump  power.  The 
scale  factors  (e.g.,  X  3)  give  the  relative  intensity  mag¬ 
nification.  At  high  temperature,  the  BEb  luminescence 
component  thermalizes,  and  the  BE|„  component  is  iso¬ 
lated.  At  high  temperatures,  weak  FE  luminescence  is 
also  observed. 


A  further  separation  of  the  BE  components  we  see  the  BE0  component  thermalize  with  respect 

which  produce  the  L  luminescence  can  be  accom-  to  the  BE|„  component.  Above  10  K,  only  the 

plished  by  measuring  the  temperature  dependence  BE(n  component  remains, 

of  the  photoluminescence  spectrum.  This  measure-  There  is  one  difficulty  with  this  interpretation 

ment  is  shown  in  Fig.  8.  We  see  that  the  high-  which  should  be  mentioned.  As  the  temperature  is 

energy  side  of  the  L  luminescence  thermalizes  as  raised  above  10  K,  the  BE|„  component  peak  posi- 

the  temperature  is  increased.  This  observation  is  non  moves  to  lower  energy,  as  expected  from  the 

consistent  with  the  interpretation  of  L  lumines-  behavior  of  BEP  luminescence  in  undoped 

cence  presented  in  the  previous  section.  At  high  Si|  _xGe, .  This  effect  is  observed  in  Fig.  8  and 

temperatures,  only  the  lower-energy  BE|„  com-  Fig.  10.  However,  the  BEtn  peak  moves  too  far;  at 

ponent  is  observed.  The  higher-energy  BEb  com-  30  K.,  it  is  already  approximately  25  meV  below 

ponent  has  thermalized.  the  FE  threshold,  while  the  BE|„  binding  ener  y  in 

This  effect  is  seen  clearly  in  Fig.  10,  which  Si  is  only  13.7  meV  11  In  spite  of  this,  the  bulk  of 

shows  the  temperature  dependence  of  L  lumines-  the  experimental  evidence  supports  the  interpreta- 

cence  which  is  gated  in  time.  At  low  tempera-  non  that  this  high-temperature  luminescence  is  due 

tures,  the  BE0  peak  is  reasonably  well  isolated.  At  to  BE|„  recombination  in  the  alloy.  It  may  be  that 

the  high  pump  powers  used  for  this  measurement,  contributions  from  FE  and  BE  TA-assisted 

however,  a  low-energy  shoulder  is  clearly  observed.  luminescence  combine  to  produce  the  anomalously 

We  interpret  this  shoulder  as  being  due  to  a  BE|„  low  high-temperature  BE,„  peak  position.  Howev- 

component  which  has  not  been  eliminated  by  the  er.  this  seems  unlikely  when  the  relative  intensities 

time  resolution.  As  the  temperature  is  increased,  of  TA  and  NP  luminescence  in  Si  are  considered. 
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It  is  also  possible  that  the  high-temperature  RE|„ 
binding 'energy  in  the  alloy  is  increased  as  a  result 
of  some  complicated  effect  which  may  involve  the 
unusually  large  radii  of  the  Ge  and  In  atoms  rela¬ 
tive  to  the  Si— Si  bond  length.  Of  course,  this 
mechanism  is  also  unlikely,  practically  by  defini¬ 
tion.  Unfortunately,  all  we  can  state  at  this  point 
is  that  the  anomalous  temperature  dependence  of 
the  BE[a  luminescence  remains  unexplained. 


V.  DISCUSSION 

In  the  previous  sections,  we  have  presented  the 
experimental  photoluminescence  data  for  the  un¬ 
doped  and  In-doped  Si-Ge  alloy  samples,  and  have 
attempted  to  interpret  the  spectra  qualitatively.  In 
this  section,  we  continue  our  qualitative  discussion 
of  the  experimental  results.  In  particular,  we  focus 
on  some  effects  which  result  from  the  composi¬ 
tional  disorder  of  the  alloy. 


A.  Bound-exciton  line  broadening 

To  begin  with,  we  see  that  the  BE  line  is  consid¬ 
erably  broadened  in  the  alloy.  In  fact,  the  BE  line 
in  the  undoped  Si|_xGe,  samples  is  approximately 
ten  times  broader  (—3  meV)  than  its  counterpart 
in  Si  (—0.4  meV).  In  general,  this  broadening  can 
be  accounted  for  on  the  basis  of  the  compositional 
inhomogeneity  of  the  alloy  material.  The  energy 
of  the  luminescence  which  results  from  the  recom¬ 
bination  of  an  exciton  bound  to  a  particular  impur¬ 
ity  will  depend  on  the  local  configuration  of  Si  and 
Ge  atoms.  This  luminescence  energy  can  be 
schematically  written  as 

£pL(k|)=£c(|c})-£FE(|c|)-£BE(|c)) .  (2) 

where  £pt  is  the  photon  energy,  £c  is  the  band- 
gap  energy,  £ is  the  free-exciton  binding  energy, 
and  £be  is  the  bound-exciton  binding  energy,  and 
where  |c)  denotes  the  particular  configuration  for 
the  impurity  being  considered. 

We  can  consider  two  simple  limits  to  the  general 
expression,  Eq.  (2).  First,  we  assume  that  the  local 
configuration  does  not  affect  the  free-exciton  bind¬ 
ing  energy  or  bound-exciton  binding  energy.  Then 
Eq.  (2)  indicates  that  the  BE-line  broadening  is  the 
result  of  variations  >n  band-gap  energy  produced 
by  local  fluctuations  in  the  alloy  composition  x. 
This  limit  was  first  discussed  by  Alferov  et  al 


We  can  roughly  estimate  this  effect  on  the  basis  of 
the  following  argument.  Given  a  volume  Ana  V3. 
there  will  be  Ana  3  <x.  N  / 3  atoms,  where  N  is  the 
density  of  lattice  sites  in  the  alloy,  if  the  Ge  atoms 
are  distributed  randomly  the  typical  fluctuation  in 
the  number  of  Ge  atoms  in  the  volume  will  be 
(4/3 na}xN)ul.  Therefore,  the  typical  fluctuation 
in  the  composition  parameter  x  will  be 


Ax- 


( jnaixN)'/2 
7  rra  iN 


(3) 


If  AE =aAx,  where  AE  is  the  change  in  gap  ener¬ 
gy  produced  by  a  change  in  composition  Ax,  then 
Eq.  (3)  becomes 


A £  —  a 


\na3N 


1/2 


(4) 


N  can  be  estimated  for  the  alloy  from  the  lattice- 
constant  measurements  of  Dismukes  et  al.13  For 
x=0. 1,  we  find  N~ 6x  lO^cm-3.  From  the  ab¬ 
sorption  measurements  of  Braunstein  et  al.,14  we 
can  estimate  a  —  600  meV.  Assuming  that  a  is  ap¬ 
proximately  an  exciton  Bohr  radius  in  Si,  i.e., 
about  40  A,  we  obtain  A E  —  5  meV,  which  is  rea¬ 
sonable  agreement  for  a  crude  calculation. 

In  spite  of  the  agreement,  this  limit  does  not  ex¬ 
plain  the  temperature-dependence  data  presented  in 
Figs.  2  and  8.  The  data  clearly  shows  that  the  BE 
peak  shifts  to  lower  energy  as  the  temperature  in¬ 
creases.  However,  the  band-gap  fluctuation  limit 
does  not  produce  this  behavior.  In  this  limit,  we 
expect  that  as  the  temperature  increases  excitons 
will  thermalize  to  bound-exciton  states  in  the  re¬ 
gions  where  the  band  gap  is  larger,  which  produces 
a  BE  peak  which  moves  to  higher  energy. 

The  second  limit  to  Eq.  (2)  we  will  consider  ac¬ 
counts  for  the  observed  temperature  dependence. 

In  this  limit,  we  assume  that  the  local  configura¬ 
tion  does  not  affect  the  band-gap  energy  or  the  FE 
binding  energy.  We  assume  that  the  gap  energy 
and  the  FE  binding  energy  are  the  bulk  average 
values  for  the  alloy.  Then  the  BE  line  broadening 
is  the  result  of  variations  in  the  BE  binding  energy 
caused  by  variations  in  the  local  configuration  of 
Si  and  Ge  atoms.  It  is  possible,  therefore,  that  the 
BE  line  is  composed  of  several  overlapping 
luminescence  lines  from  BE's  with  slightly  dif¬ 
ferent  binding  energies.  These  binding  energies 
will  increase,  in  principle,  from  the  BE  binding  en¬ 
ergy  in  Ge  to  the  BE  binding  energy  in  Si. 

Qualitatively,  this  limit  produces  the  tempera¬ 
ture  dependence  observed  in  Figs.  2  and  8.  At  low 
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temperatures,  we  see  luminescence  from  all  BE’s. 
As  the  temperature  increases,  though,  we  expect 
that  luminescence  from  the  most  tightly  bound 
BE’s  will  dominate.  As  we  have  previously  men¬ 
tioned,  these  BE’s  will  have  the  most  Si-like  local 
configuration.  In  addition  to  explaining  the 
temperature-dependence  data,  this  picture  of  BE- 
line  broadening  is  consistent  with  time-resolved 
spectra  of  the  alloy  luminescence.  An  analysis  of 
the  relevant  decay  times  shows  that  the  decay  of 
all  components  of  BE  luminescence  is  capture- 
limited  in  the  cases  studied  here.13  This  is  con¬ 
sistent  with  our  observation  that  all  components  of 
the  BE  line  decay  at  the  same  rate  —  the  line 
shape  does  not  change  and  the  peak  position 
remains  fixed. 


B.  Comparison  of  Si  and  Si-Ge 
luminescence  energies 

If  we  assume  that,  at  high  temperatures,  the  BE 
luminescence  is  dominated  by  components  which 
have  the  most  Si-like  local  configuration,  then  a 
procedure  for  directly  comparing  Si  and  Si,_xGex 
luminescence  energies  presents  itself.  We  propose 
that  the  difference  in  recombination  energy  be¬ 
tween  the  high-temperature  BE-peak  position  in 
the  Si-Ge  alloy  and  the  corresponding  BE-peak  po¬ 
sition  in  Si  is  just  due  to  the  shift  in  band  gap  of 
the  alloy. 

Figure  1 1  schematically  illustrates  the  conse¬ 
quences  of  this  proposal  for  the  x  =0. 1 1  Si|_xGe, 
sample,  C077.  The  high-temperature  BEp  lumines¬ 
cence  line  position,  1 108.71  meV,  was  obtained 
from  the  8.0-K  spectrum.  This  was  the  highest 
temperature  at  which  BEp  luminescence  could  be 
reliably  observed,  and  it  was  assumed  that  the 
luminescence  was  dominated  by  recombination  of 
BEp’s  in  essentially  a  Si  environment.  The  P 
luminescence  line  position  in  Si,  1150.11  meV,  was 
obtained  from  direct  measurement  of  lightly  doped 
Si:P.  If  this  change  in  line  position  is  attributed 
solely  to  a  change  in  gap  energy,  these  measure¬ 
ments  imply  a  band-gap  decrease  of  41.4  meV  for 
x  =0. 1 1.  This  value  fluctuated  by  about  0.6  meV 
due  to  large-scale  inhomogeneities  in  the  crystals 
being  examined.  If  a  rough  graphical  extrapolation 
of  the  absorption  data  published  by  Braunstein 
et  al.{*  is  attempted,  we  obtain  an  alloy  band-gap 
shift  of  approximately  60  meV  for  x  =0. 1 1,  which 
is  in  reasonable  agreement  with  the  luminescence 
data. 


SCHEMATIC  COMPARISON  OF 
UNDOPED  Si  AND  UNDOPED  Si,.,Ge, 
I UMINFSCENCE  ENERGIES 
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UNDOPED  Si 

(SAMPLE  C 20204  A) 

Schematic  line-position  comparison  between 

Si  and  Si,_,Ge,  sample  C077.  This  comparison  is 
based  on  the  assumption  that  the  difference  between  the 
energy  of  high-tempcrature  BEp  alloy  luminescence  is 
equal  to  the  difference  in  band-gap  energies  for  Si  and 
Sit_,Ge,  sample  077.  This  fixes  the  alloy  gap  energy. 
The  spectra  are  then  compared  by  lining  up  the  Si  and 
Si|_,Ge,  gap  energies. 


In  addition  to  the  high-temperature  Sij_xGe, 
BEp  line  position  and  the  Si  BEp  line  position.  Fig. 
1 1  also  shows  the  FE  threshold  and  p  positions 
for  Si  and  the  x  =0. 1 1  alloy.  The  FE  threshold 
position  in  the  alloy,  1 1 14.64  meV,  was  obtained 
by  fitting  Eq.  (1)  to  the  high-temperature  FE 
luminescence  spectra  (Fig.  3),  as  illustrated  in  Fig. 
4.  The  NP-FE  threshold  position  in  Si,  1154.68 
meV,  was  obtained  by  fitting  Eq.  (1)  directly  to 
NP-FE  luminescence  from  pure  Si.  The  Si  gap  en¬ 
ergy,  1169.38  meV,  was  determined  from  the  NP- 
FE  threshold  energy  by  adding  the  FE  binding  en¬ 
ergy,  which  is  14.7  meV.  Then  the  alloy  gap  ener¬ 
gy,  1 127.98  meV,  was  obtained  by  subtracting  the 
band-gap  shift,  41.4  meV,  from  the  Si  gap  energy. 

As  indicated  in  Fig.  11,  this  method  implies  that 
the  FE  dissociation  energy  in  the  x  =0. 1 1  alloy  has 
decreased  by  1.36  meV  from  its  value  in  Si,  to 
13.34  meV.  This  compares  favorably  with  a  rough 
estimate  based  on  a  linear  interpolation  between 
the  FE  dissociation  energies  in  Si  and  Ge,  which 
are  14.7  and  4.15  meV,  respectively.  The  linear  in- 


5360 


G.  S.  MITCHARD  AND  T.  C.  McGILL 


25 


terpolation  results  in  an  estimate  of  13.54  meV  for 
the  FE  dissociation  energy  in  the  x  =0. 1 1  alloy. 

A  similar  analysis  has  been  applied  to  the 
x  =0.067  sample,  C021-3.  As  shown  in  Fig.  7, 
compared  to  sample  C077  the  luminescence  from 
sample  C021-3  is  shifted  to  higher  energy  by  15.5 
meV.  This  implies  a  band-gap  decrease  of  25.9 
meV  for  the  x  =0.067  alloy.  This  result  and  the 
result  previously  obtained  for  sample  C077  are 
consistent  with  the  alloy-composition  measure¬ 
ments  presented  in  Table  I  if  a  linear  relationship 
between  the  alloy  composition  and  change  in  band 
gap  is  assumed.  The  composition  ratio  between 
the  two  samples  is  0.609  while  the  ratio  of  the 
measured  band-gap  shifts  is  0.626. 

To  compare  the  Si  and  Si^^Ge,  luminescence 
energies  for  the  In-doped  sample,  C093,  we  note 
that  samples  C077  and  C093  differ  in  composition 
by  only  about  1%.  We  therefore  assume  that  the 
FE  binding  energies  in  the  two  samples  are  the 
same.  This  leads  to  the  comparison  schematically 
illustrated  in  Fig.  12.  The  gap  energy  for  sample 
C093,  1 129.14  meV,  was  determined  by  adding  the 
assumed  FE  binding  energy,  13.34  meV,  to  the  FE 


SCHEMATIC  COMPARISON  OF 
SI :  in  AND  Si,.,Ge„:  In 
LUMINESCENCE  ENERGIES 


Si.-«GVn 


(SAMPLE  C093) 

FE 

(MI5. 8) 

.  ,  .  J  , 

Egoe 

(1129.14) 

1 

MOO 

i  i 

mo  M20 

1  l 

11130 

ENERGY  (meV) 

1140 

i-r 

H50  M60 

t  i 

ill  70 

0 140.981 

^  r 

U150.U)  UI5466) 

t 

(1169.38! 

BE,, 

BEP  FE 

Eqap 

S t  in 


(SAMPLE  Z  10401  • ) 

FIG.  12.  Schematic  line-position  comparison  between 
Si  mui  Sit  ,Gc,:In  sample  C093.  This  comparison  is 
bused  on  the  aMuitiplion  (hut  the  I  T.  Iiimhtig  energy  in 
sample  CW.)  and  sample  C'077  is  identical.  flits  fixes 
the  alloy  gup  energy.  The  spcclru  are  then  compared  by 
lining  up  the  Si  and  Sit  ,Gc,  :ln  gap  energies. 


threshold  energy,  1 1 15.8  meV.  This  results  in  a 
measured  band-gap  shift  of  40.2  meV  for  the 
x  =0. 104  alloy.  The  band-gap  shifts  and  alloy- 
composition  measurements  are  reasonably  con¬ 
sistent,  assuming  a  linear  relationship  between  the 
alloy  composition  and  change  in  band  gap.  The 
composition  ratio  between  samples  C077  and  C093 
is  0.920  while  the  ratio  of  measured  band-gap 
shifts  is  0.971.  Samples  C021-3  and  C093  agree 
exactly  —  the  composition  ratio  and  the  ratio  of 
measured  band-gap  shifts  are  both  0.644. 

C.  Excitons  bound  to  composition  fluctuations 

One  feature  in  the  luminescence  spectrum  of 
sample  C077  remains  to  be  discussed.  At  high 
temperatures  (Fig.  3),  the  FE  line  assumes  a  shape 
characteristic  of  FE  recombination  in  Si  and  is 
well  described  by  Eq.  (1),  as  previously  discussed. 
The  FE  threshold  energy  oHained  by  fitting  Eq. 

(1)  is  independent  of  temperature,  as  expected. 
However,  at  low  temperatures  this  is  no  longer  the 
case.  Spectra  taken  at  various  temperatures  below 
6  K  are  shown  in  Fig.  13.  We  see  that  as  the 
temperature  decreases,  the  line  shifts  to  lower  ener¬ 
gy  and  assumes  a  shape  which  is  no  longer  charac¬ 
teristic  of  FE  recombination. 

These  measurements  of  low-temperature  FE 
luminescence  suggest  that  a  low-energy  shoulder 
appears  on  the  FE  luminescence  at  low  tempera¬ 
tures.  This  shoulder  may  be  due  to  the  recombina¬ 
tion  of  excitons  bound  by  only  about  0. 1  meV  to  a 
very  shallow  level.  This  interpretation  is  supported 
by  the  temperature-resolved  spectra,  which  show 
that  the  shoulder  has  completely  thermalized  by  6 
K.  Above  6  K,  Eq.  (1)  correctly  describes  the  line 
shape,  and  the  threshold  energy  is  independent  of 
temperature.  This  behavior  is  consistent  with  a 
binding  energy  of  tenths  of  meV.  Also,  the 
pump-power  dependence  at  1.6  K  (Fig.  14)  shows 
that  the  FE  to  BEP  intensity  ratio  is  independent 
of  pump  power,  as  expected  for  independent 
exciton-binding  centers  in  the  absence  of  satura¬ 
tion. 

These  features,  and  the  fact  that  there  is  no 
known  defect  or  impurity  in  Si  which  binds  an  ex- 
citon  by  only  0.1  meV,  leads  us  to  suggest  that  the 
new  feature  may  be  due  to  the  recombination  of 
excitons  weakly  bound  to  local  fluctuations  in  the 
alloy  composition  (HE, ).  The  poxxibilily  of  such  a 
stale  was  firxl  pointed  out  by  Uunmovskii  and 
LIVos,"'  who  extended  the  bulk  fluctuation  theory 
of  Alferov  ct  at ,IJ  to  include  excitons.  A  line 
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FIG.  13.  Temperature-resolved  spectra  of  undoped 
Si,_,Ge,  sample  C077  at  low  temperatures.  The  scale 
factors  (e.g.,  X  2)  give  the  relative  intensity  magnifica¬ 
tion.  Note  that  the  free-exciton  line  shifts  and  broadens 
as  the  temperature  is  decreased.  For  reference,  the 
free-exciton  threshold  energy  measured  at  high  tempera¬ 
ture  is  shown. 


FIG.  14.  Pump-power  dependence  of  luminescence 
from  undoped  Sii_zGex  sample  C077  at  1.6  K.  The 
scale  factors  (e.g.,  X  36)  give  the  relative  intensity  mag¬ 
nification.  Note  that  the  BEp-FE  intensity  ratio  is  in¬ 
dependent  of  pump  power. 


which  may  be  due  to  fluctuation  bound  excitons 
has  just  recently  been  reported  in  GaAS|_xP,,17 
and  has  characteristics  which  are  similar  to  those 
of  the  line  reported  here. 

One  of  the  most  important  features  of  this 
model  is  the  prediction  that  BEF  luminescence 
should  have  a  very  long  lifetime.  As  indicated 
during  the  discussion  of  BEP  decay,  the  dominant 
decay  mechanism  for  donor  or  acceptor  BE's  in  Si 
is  the  Auger  process.  This  results  in  observed  de¬ 
cay  rates  which  can  be  orders  of  magnitude  larger 
than  the  radiative  rates  obtained  from  measured 
oscillator  strengths.10'18  For  instance,  the  BEP  ra¬ 
diative  decay  rule  is  estimated  to  Ik  5x  I0j  sec- ', 
whereas  t lie  observed  Auger  decay  rale  is  .1.7 
x  10ft  sec  l.m  When  an  cxciton  is  bound  to  a  lo¬ 
cal  composition  fluctuation,  however,  the  Auger 


process  is  not  possible  and  long-lifetime  lumines¬ 
cence  is  expected.  Indeed,  this  was  the  case  for  the 
luminescence  feature  observed  in  GaAsi_xPx.16 

Figure  15  presents  the  best  low-temperature 
time-resolved  spectra  available  for  the  Si0  pGeo.  i  al¬ 
loy.  Although  the  spectra  are  somewhat  noisy,  we 
see  that  the  BEP  luminescence  does  not  appear  to 
have  a  particularly  long  lifetime  compared  to  the 
BEP  luminescence  line.  This  apparent  discrepancy 
can  be  investigated  by  estimating  the  relevant  de¬ 
cay  times.  These  estimates  show  that  in  this  case 
we  are  in  a  regime  where  the  decays  are  capture 
limited,15  which  explains  why  all  components  in 
the  spectrum  appear  to  decay  with  the  same  rate, 
f  fence,  while  the  low-temperature  time-resolved 
speciru  do  not  show  long  lifetimes  for  the  BEf 
luminescence,  the  data  ure  not  inconsistent  with 
the  model  proposed  above. 
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FIG.  15.  Time-resolved  spectra  of  undoped  Sit_,Ge, 
sample  C077  at  1.6  K.  The  scale  factors  (e.g.,  X4)  give 
the  relative  intensity  magnification.  Note  that  FE 
luminescence  does  not  dominate  the  long-time  spectrum. 

VI.  SUMMARY  AND  CONCLUSION 

In  this  paper,  we  have  presented  the  results  of  a 
study  of  the  photoluminescence  spectra  of  Si-rich 
Sii-jGe*  alloy  materials.  Undoped  and  In-doped 
samples  were  considered  in  the  range  0.067  <;x 


^0.11.  The  luminescence  features  observed  in  the 
spectra  were  identified  as  free-exciton  and  bound- 
exciton  luminescence.  The  free-exciton  lumines¬ 
cence  is  described  very  accurately  by  the  Si  no¬ 
phonon  free-exciton  line  shape,  Eq.  (1).  Asa  re¬ 
sult,  a  no-phonon  free-exciton  threshold  energy  was 
obtained  for  each  sample. 

In  addition  to  this  identification  of  the  lumines¬ 
cence  features,  mechanisms  for  the  observed 
bound-exciton  —  luminescence  broadening  were  dis¬ 
cussed.  On  the  basis  of  this  discussion,  the  Si-Ge 
luminescence  energies  were  compared  with  those 
from  Si,  and  values  for  the  alloy  band-gap  shift 
and  change  in  free-exciton  binding  energy  were  ob¬ 
tained.  Finally,  evidence  was  presented  which  may 
indicate  that  luminescence  from  excitons  bound  to 
local  composition  fluctuations  is  observed  in  the 
low-temperature  luminescence  spectra. 
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Vibrational  modes  of  oxygen  in  GaP  including  nearest-neighbor  interactions 

R.  M.  Feenstra,*  R.  J.  Hauenstein,  and  T.  C.  McGill 
California  Institute  of  Technology,  Pasadena,  California  91125 
(Received  10  March  1983;  revised  manuscript  received  16  August  1983) 

A  theoretical  treatment  of  the  vibrational  modes  of  substitutional  defects  in  zinc-blende  crystals  is 
presented.  The  dynamics  of  the  perfect  lattice  are  described  with  the  use  of  a  two-parameter  model 
including  nearest-neighbor  interactions  only.  The  defect  consists  of  an  impurity  atom  with  variable 
strength  nearest-neighbor  force  constants.  For  the  case  of  oxygen  in  GaP,  the  theoretical  defect  vi¬ 
brational  modes  are  compared  with  experiment.  Agreement  between  theory  and  experiment  is  ob¬ 
tained  for  values  of  the  defect  force  constants  which  are  13%  of  the  bulk  values  for  the  O0  and  O* 
defects,  indicating  that  the  oxygen  impurity  is  quite  weakly  bonded  in  the  lattice. 


I.  INTRODUCTION 

An  oxygen  atom  on  a  P  site  in  GaP  is  a  defect  with 
many  interesting  and  controversial  properties.  The  vibra¬ 
tional  states  of  the  defect  are  observed  in  phonon-assisted 
transitions  appearing  in  a  variety  of  optical  spectra.1-4 
From  the  relatively  low  energy  of  those  phonon  modes 
which  involve  oxygen  motion,2  it  seems  that  the  oxygen 
impurity  is  quite  weakly  bonded  in  the  lattice,5  and  recent 
theoretical  studies6  support  this  conclusion.  In  this  paper, 
we  theoretically  examine  the  vibrational  modes  of  an  oxy¬ 
gen  defect  with  Td  symmetry.  We  use  a  phenomenologi¬ 
cal  model  for  the  lattice  dynamics  with  the  force  con¬ 
stants  between  the  defect  and  neighboring  atoms  being 
treated  as  parameters  in  the  theory.  We  compare  our  re¬ 
sults  with  experiment  for  the  O0  and  0+  defects  (with  and 
without  the  first  donor  electron  present,  respectively).  We 
find  reasonable  agreement  between  theory  and  experiment 
for  defect  force  constants  which  are  about  15%  of  the 
bulk  values,  thus  supporting  the  current  view  that  the  0° 
and  0+  defects  are  weakly  bonded  in  a  Td  symmetric  site. 

In  this  paper  the  dynamics  of  the  perfect  lattice  are 
described  using  a  two-parameter  model  which  includes 
nearest-neighbor  interactions  only.  By  symmetry,  only 
two  parameters  are  required  for  the  most  general  descrip¬ 
tion  of  the  bulk  phonons  in  this  model,7  and  these  param¬ 
eters  are  chosen  to  fit  known  phonon  dispersion  curves.  A 
substitutional  defect  consists  of  an  impurity  atom  of  given 
mass  with  variable-strength  nearest-neighbor  force  con¬ 
stants.  The  properties  of  a  defect  are  thus  determined  by 
only  two  parameters.  Therefore,  we  are  able  to  obtain 
quite  general  results  for  defects  within  our  very  simplified 
model  for  the  bulk  crystal.  More  complicated  models  are 
available  for  the  bulk  phonons,  and,  in  particular,  the  15- 
parameter  deformable-dipole  model  of  Kune*  seems  to 
give  a  good  description  of  many  III-V  compounds.  How¬ 
ever,  within  this  model  a  substitutional  defect  is  also 
described  by  15  (or  more)  parameters,  and  the  results  for 
the  defect  vibrational  modes  may  be  somewhat  dependent 
on  the  choice  of  parameter  values.  Thus,  the  complicated 
models  yield  specific  results  for  a  defect  within  a  realistic 
model  for  the  bulk.  Because  of  this  phenomenological  na¬ 
ture  of  our  calculations  we  feel  that  it  is  necessary  to  ex- 
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amine  the  results  from  both  the  simple  and  complicated 
dynamical  models.  In  this  paper  we  present  results  from 
the  simple  two-parameter  model  and  results  from  the 
deformable-dipole  model  are  planned  to  be  the  topic  of  a 
forthcoming  publication.9 

Previous  workers5,6-10  have  considered  the  vibrational 
states  of  the  GaP:0  defect  using  molecular-type  models  of 
the  defect  and  its  immediate  neighbors.  Here,  we  use  the 
Green's-function  technique,"  which  allows  us  to  compute 
the  vibrational  states  of  the  defect  including  the  response 
of  the  entire  lattice.  The  advantage  of  this  technique  is 
that  it  allows  us  to  couple  the  vibrations  of  the  defect  into 
the  bulk  crystal  and  see  which  types  of  vibrational  modes 
remain  localized  around  the  defect.  The  term  “localized” 
will  be  used  here  to  refer  to  both  strictly  localized  vibra¬ 
tional  modes  and  resonant  vibrational  modes  in  which  the 
vibrational  amplitude  is  relatively  peaked  near  the  defect. 
The  quantity  we  will  use  to  describe  the  vibrational  modes 
of  a  defect  is  the  local  density  of  states  (LDOS).  Strictly 
localized  modes  show  up  as  5  functions  in  the  LDOS  and 
resonant  modes  appear  as  peaks  with  some  nonzero  width. 

In  Sec.  II  we  present  the  details  of  our  computational 
technique.  Our  method  is  quite  similar  to  other  vibration¬ 
al  Green’s-function  computations.7,12  The  results  of  our 
calculations  are  presented  in  Sec.  III.  We  show  the  com¬ 
puted  dispersion  curves  and  density  of  states  (DOS)  for 
the  bulk  phonons  in  GaP.  A  model  is  used  for  the  substi¬ 
tutional  O  defect  in  which  the  nearest-neighbor  O-Ga 
force  constants  are  all  varied  by  the  same  amount.  For 
various  values  of  the  defect  force  constant  the  LDOS  for 
A  |  and  T 2  defect  vibrational  modes  are  presented.  We 
find  basically  one  A ,  and  one  T 2  vibrational  mode  of  the 
defect.  In  Sec.  IV  we  compare  our  theoretical  results  with 
experiment.  For  defect  force  constants  which  are  15%  of 
the  bulk  values,  we  identify  our  computed  T2  mode  with 
the  24.7-meV  resonance  observed  for  the  O0  defect.2  Our 
computed  A\  mode  is  identified  with  the  19-meV  mode 
observed,1,3  for  the  O0  and  0+  defects.  We  argue  that 
these  phonons  are  the  most  localized  models  of  the  defect. 
Other  modes  observed  in  the  optical  spectra  are  either 
bulklike  phonons  or  relatively  delocalized  defect  modes. 
The  conclusions  which  can  be  drawn  from  this  work  are 
presented  in  Sec.  V. 
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II.  COMPUTATIONAL  METHOD 


A.  Bulk  phonons 


The  lattice  dynamics  of  the  perfect  crystal  are  described 
using  a  two-parameter  model  consisting  of  nearest- 
neighbor  interactions  only.  As  shown  by  Grimm  et  al.,1 
in  the  harmonic  approximation,  invariance  of  potential 
energy  with  respect  to  rigid-body  translations,  rotations, 
and  symmetry  operations  of  the  crystal  lattice  result  in 
only  two  parameters  required  to  describe  the  bulk  pho¬ 
nons.  In  this  model  the  matrix  of  second  derivatives  of 
the  potential  energy  with  respect  to  atomic  positions  has 
diagonal  Cartesian  blocks  of  the  form 


<t>(  Ga;Ga )  =  <t>(  P;P)  = 


4 A 
0 
0 


0 

4A 

0 


0 

0 

4 A 


(la) 


For  convenience  we  decompose  the  parameters  A  ar.d  B 
into  bond-stretching  and  bond-bending  components.13 
The  contribution  to  the  potential  energy  of  the  system 
from  a  bond  (spring)  connecting  atoms  i  and  j  is  written 
as 

^  =  7/1  I  (Ui-u;)-eiy  |  2+  y/2  I  ( u,  —  u j)xeij  | 2  ,  (2a) 

with 


en  = 


R  ,-R; 

I  R|  R;  I 


(2b) 


where  u,  is  the  displacement  of  the  ith  atom  from  its 
equilibrium  position  R,.  The  parameters  /,  and  /2  are 
defined  as  bond-stretching  and  bond-bending  force  con¬ 
stants,  respectively,  and  are  related  to  A  and  B  bv 
A  =(/i  +2/2)/3  and  B  =(/,  — /2)/3. 


and  off-diagonal  blocks  typically  of  the  form 


B.  Green’s  functions 


d>(Ga;P)=<t>(P;Ga)=- 


A  B  B 
BAB 
B  B  A 


From  the  eigenvalues  and  eigenvectors  of  the  dynamical 
dkj  matrix  we  compute  the  perfect  crystal  Green’s  function 
G°.  Maradudin 1 1  has  shown  that  the  matrix  elements  of 
G°  can  be  expressed  as 


ag,  2,  1  ^  £(/i,k  )£f(n,k  )exp[/ic-(R(a— Rmfc)] 

Glmab(w  )=—==■  lim  2 - T— - 2 - 


, -  Mill  ^ 

Nv  mamb  c—*o+  - 


a)2  +  ie  —  a>2  - 

n,  k 


(3) 


where  /  and  m  label  the  N  unit  cells,  a  and  b  label  atomic 

types,  and  a  and  &  label  Cartesian  coordinates,  a)2  -  are 

_  * 

the  eigenvalues  and  £J(n,  k)  are  the  eigenvectors  of  the 
(n,  k)th  normal  mode.  The  equilibrium  position  of  an 
atom  is  denoted  by  Rta,  and  ma  is  the  mass  of  the  <zth 
atom. 

To  evaluate  the  LDOS  of  a  crystal  containing  a  defect  it 
is  necessary  to  compute  that  portion  of  the  Green’s  func¬ 
tion  which  is  in  the  space  of  the  defect.  The  defects  we 
consider  here  consist  of  a  substitutional  impurity  atom  in 
a  zinc-blende  crystal  with  bonds  of  variable  strength  con¬ 
necting  the  impurity  atom  to  its  neighboring  atoms.  This 
defect  consists  of  a  five-atom  cluster  with  Td  symmetry. 
The  positions  of  the  five  atoms  are  specified  by  1 5  Carte¬ 
sian  coordinates,  which  form  a  15-dimensional  representa¬ 
tion  of  the  Td  group.  This  reducible  representation  T  can 
be  decomposed  into  the  irreducible  representation  of  the 
Td  group  according  to 

r =A  |  +E  +  T\  +  37*2  .  (4) 


I 

E  and  T 1  coordinates  involve  motion  of  the  outer  atoms 
with  the  inner  atom  stationary.  The  definition  of  all  the 
collective  coordinates  is  given  in  Appendix  A. 

To  produce  a  Green’s  function  suitable  for  computation 
we  apply  two  coordinate  transformations.  The  first 
transformation  is  quite  trivial  and  its  effect  is  to  drop  the 
masses  ma  and  mh  from  Eq.  (3)  and  to  absorb  them  into 
the  defect  perturbation  matrix  5 L.  This  corresponds  to 
expressing  G°  with  respect  to  the  reduced  coordinates 
<?,  =  ^/ntj  u,  where  u,  are  Cartesian  coordinates.  Second¬ 
ly,  we  transform  to  the  collective  coordinates  Q,  defined 
in  Appendix  A.  In  this  coordinate  system  ( s )  and  (r)  label 
the  irreducible  representation,  a  and  p  label  the  occurrence 
of  the  representation,  and  p  and  v  label  the  row  or  column 
of  the  representation.  The  Green’s  function  in  the  space 
of  the  defect  can  now  be  expressed  as 


,  2v  1  Q^(n,k)Q:r%(n,k) 

Hr)op(u l)=—  hm  2  - 2 - : - 7 - 

H  f^o*  r  co2  +  i€- a>*- 

«.  k  n,  k 


(5) 


This  equation  classifies  the  collective  coordinates  of  the 
five-atom  cluster.  The  A 1  coordinate  consists  of  a 
"breathing”  type  of  motion  in  which  the  central  atom  is 
stationary  and  the  outer  atoms  move  radially.  The  motion 
of  the  central  impurity  atom  forms  a  basis  for  a  T2  repre¬ 
sentation  of  the  Td  group.  The  other  two  sets  of  Tz  coor¬ 
dinates  involve  motion  of  the  outer  atoms.  Similarly,  the 


The  summation  on  k  extends  over  the  entire  first  Bril- 
louin  zone.  This  can  be  divided  into  a  sum  over  the  ir¬ 
reducible  wedge  W  of  the  zone  and  a  sum  over  the  group 
elements  R  which  define  the  star  of  k.  Since  the  Q  coor¬ 
dinates  form  bases  for  irreducible  representations  of  the 
Td  group,  they  are  orthogonal  when  summed  over  R ,  and 
we  find  that 
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where  h  -  is  the  number  of  points  in  the  star  of  k  and  d, 
is  the  dimension  of  the  representation  (s)  [compare  with 
Eq.  (41)  of  Ref.  14].  We  see  that  the  Green’s  function  is 
zero  unless  (s)=(r)  and  n=v,  in  which  case  it  is  indepen¬ 
dent  of  it.  Thus  the  Green’s  function  has  been  block- 
diagonalized  by  the  transformation  to  the  collective  coor¬ 
dinates.  The  transformed  matrix  consists  of  lxl,  2x2, 
3x3,  and  9x9  blocks  for  the  A\,  E,  Tlt  and  37"2  repre¬ 
sentations,  respectively.  Furthermore,  these  matrix  ele¬ 
ments  do  not  depend  on  the  row  or  column  of  the  repre¬ 
sentation,  so  that  these  blocks  can  be  further  reduced  to 
1X1,  lxl,  1X1,  and  3x3  symmetric  blocks,  respective¬ 
ly,  for  a  total  of  nine  independent  matrix  elements.  Keep¬ 
ing  the  zero  elements  in  mind  we  can  simplify  the  nota¬ 
tion  by  dropping  the  subscripts  ft,  v,  and  (r).  Further¬ 
more,  the  Green’s  function  can  be  explicitly  divided  into 
its  real  and  imaginary  parts  by  defining  the  partial  DOS 
(“spectral  density”), 

l  A  - 

2  2 Q(s^n,k)Q^p(nS) 

v  <i.  r  €  w  a’  * 

X 8(<a2  —  <u2  -)  ,  (7) 

W,  K 

from  which  Eq.  (6)  can  be  expressed  as14,15 

G?,l<rpicoz)= p  r  ]  (a>*) ,  (8) 

u  CO  — (O 

where  P  signifies  a  principle-value  integral.  Equations  (7) 
and  (8)  are  in  a  form  suitable  for  computation.  The  sum 
over  k  (in  the  first  Brillouin  zone)  was  reduced  to  a  sum 
over  the  irreducible  wedge  (-jj-th  of  the  Brillouin  zone  for 
the  zinc-blende  structure)  and  was  performed  using  the 
Gilat-Raubenheimer  method.16,17  The  principal-value  in¬ 
tegral  was  then  computed  following  Bemholc  and  Pan- 
telides.15  The  results  presented  here  were  computed  using 
5950  points  in  the  irreducible  wedge  corresponding  to 
256000  points  in  the  entire  zone.  'The  computations  were 
all  performed  as  a  function  of  <u2,  using  an  interval  size  of 
4.12  meV2.  This  corresponds  to  an  accuracy  of  about  0.05 
meV  in  the  optical  branches  and  0.5  meV  near  the  bottom 
of  the  acoustic  branch.  The  number  of  points  used  in  the 
k -space  summation  was  consistent  with  this  a)1  interval 
size.  For  plotting  purposes  the  Green’s  functions  were 
converted  to  a  linear  energy  scale  by 

G(aj)da>  =  G((i)2)dcj2  ,  (9) 

where  G  are  the  functions  which  will  be  shown  here. 

The  Green’s  function  for  the  imperfect  crystal  in  the 
space  of  the  defect  is  computed  from11 


G  =(1  — G°6L)-IG°  ,  (10) 

where  G°  and  6 L  are,  respectively,  the  Green’s-function 
and  defect-perturbation  matrices  in  the  15x15  space  of 
the  defect.  The  diagonal  elements  of  —  ImG/jr  give  the 
LDOS  of  the  imperfect  crystal.  For  representations  which 
have  multiple  occurrences  in  the  space  of  the  defect,  the 
LDOS  for  each  occurrence  provides  a  measure  of  the  frac¬ 
tion  of  total  kinetic  energy  in  that  particular  type  of 
motion,  and  this  quantity  is  used  in  the  following  sections 
to  describe  the  composition  of  the  7*2  defect  modes. 

C.  Defect-perturbation  matrix 

The  defect-perturbation  matrix  8L  is  defined  by11 

5L  =  —  o)2&M+  A<t> ,  (11) 

where  AM  and  A4>  are,  respectively,  the  mass-defect  and 
force-constant-defect  matrices.  As  discussed  by  Maradu- 
din,11  A<t>  must  be  invariant  under  rigid-body  translations, 
rotations,  and  the  symmetry  operations  of  the  imperfect 
crystal.  In  the  present  model,  assuming  the  same  value  of 
the  force  constants  change  AA  and  A B  for  each  of  the 
four  bonds  surrounding  the  mass  defect,  the  defect- 
perturbation  matrix  has  the  form  given  by  Talwar,  Vande- 
vyver,  and  Zigone  (Table  2  of  Ref.  12).  The  Ad>  com¬ 
ponent  of  this  matrix  explicitly  satisfies  invariance  with 
respect  to  translations  and  crystal  symmetry  operations. 
To  further  satisfy  rotational  invariance  it  is  necessary  to 
assign  A A  =  A B  (At  —Bt  in  the  notation  of  Ref.  12),  or 
equivalently  the  change  in  bond  bending  must  be  zero, 
A/2  =0.  The  physical  origin  of  such  a  constraint  is  as  fol¬ 
lows:  Bond  bending  clearly  violates  rotational  invariance 
for  an  isolated  XY 4  molecule  in  space  since  a  rotation  of 
the  molecule  produces  a  nonzero  energy  change,  and  the 
invariance  constraints  for  such  an  isolated  molecule  are 
identical  with  those  for  the  “defect  molecule”  considered 
here.  However,  in  considering  very  weakly  bonded  defects 
one  would  like  to  uniformly  reduce  both  the  bond- 
stretching  and  bond-bending  force  constants  to  zero  to 
produce  a  defect  which  experiences  zero  restoring  force 
(i.e.,  A/|  // 1  =  A fi  //2  — *•  —  1  or  equivalently  AA  /A 
=  AB/B—  —  1).  In  order  to  arbitrarily  vary  the  bond¬ 
bending  force  constant  it  is  necessary  to  introduce  an  ad¬ 
ditional  defect  force  constant  which  has  the  form  of  a 
second-nearest-neighbor  bond-bending  interaction  [as  de¬ 
fined  in  Eq.  (2a)]  among  the  four  atoms  surrounding  the 
mass  defect.  Denoting  this  force  constant  by  ju,  rotational 
invariance  is  satisfied  by  the  condition  n=  —  A/2/4.  The 
final  form  of  the  defect  perturbation  matrix  is  given  in 
Appendix  B.  We  note  that  the  constraint  of  rotational  in¬ 
variance  affects  our  results  only  through  bond-bending  in¬ 
teractions  which  are  relatively  small  for  the  GaP  lattice. 
In  a  preliminary  report18  of  our  results  we  neglected  rota¬ 
tional  invariance,  and  a  comparison  between  those  results 
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and  the  results  to  be  presented  here  reveals  that  the  effects 
of  rotational  invariance  condition  are  indeed  negligible. 
Nevertheless,  rotational  invariance  is  a  fundamental  prop¬ 
erty  of  the  system,  and  in  this  paper  we  completely  in¬ 
clude  its  effects. 

III.  RESULTS 


A.  Bulk  phonons 

Shown  in  Fig.  1  are  phonon  dispersion  curves  and  DOS 
computed  from  the  two-parameter  model  described  in  Sec. 
II  A.  The  values  of  the  model  parameters  /,  and  / 2  have 
been  chosen  to  produce  a  reasonable  fit  to  the  actual  pho¬ 
non  energies  of  GaP,  by  matching  the  TA*  and  TOp  ener¬ 
gies.  From  the  dynamical  matrix  given  by  Grimm  et  al.1 
these  energies  are  found  to  be 
,  4  A 

ct>2(TOr)  =  — 7-  (12a) 

M 


and 


«r(TA*)  = 


M 

M 


\_ 

2 


_j_ 

2 


1- 


AM 

ma+mb 


(l  2b) 


where  A  and  B  are  specified  in  Eq.  (1)  and 
M  =mamh/(ma  +mh)  with  a  denoting  Ga  and  b  denoting 
P.  For  GaP  the  actual  energies  are  known  to  be19 
£(TA*)=13.1  meV  and  £(TOr)  =  45.4  meV,  from 
which  the  values  for  the  force  constants  are  determined  to 
be/,  =6.87  eV/A2  and /2  =0.525  eV/A2.  Our  dispersion 
curves  can  be  compared  with  the  shell-model  results  of 
Yamell  et  al.20  and  the  deformable-dipole  model  results  of 
Kune.8  The  major  deficiency  of  the  present  computation 
is  that  the  LO  branch  (near  40  meV)  should  actually  be 
higher  in  energy  due  to  electrostatic  interactions  which 


B.  Perfect-crystal  Green’s  functions 

Shown  in  Fig.  2  are  the  elements  of  the  Green’s- 
function  matrix  for  the  perfect  crystal  as  computed  from 
Eqs.  (6)  and  (8)'.  These  functions  are  computed  using  a  P 
site  as  the  origin.  We  plot  —  ImG°/tr  for  the  At,  E,  T,, 
and  the  diagonal  elements  of  the  T2  representations. 
Shown  on  the  top  of  Fig.  2  is  the  bulk  DOS.  The 
Green’s-function  elements  give  the  partial  DOS  for  each 
representation.  The  appearance  of  the  Green’s-function 
elements  may  be  understood  by  considering  the  type  of 
atomic  motion  they  represent.  For  example,  the  A  , 
Green’s  function  is  a  breathing  motion  which  is  made  up 
entirely  from  longitudinal  phonons.  With  our  P-site  ori¬ 
gin,  the  A 1  Green’s  function  is  very  small  in  the  optical 
branch  because  the  optical  modes  involve  very  little 
motion  of  the  Ga  atoms  due  to  the  mass  difference  be¬ 
tween  Ga  and  P.  The  Ga-site  Green’s  function  elements 
can  be  obtained  from  the  P-site  elements  by  reflecting  the 
G(cj2)  functions  through  a  mirror  plane  located  at  4j^,„/2, 
where  occurs  at  the  top  of  the  optical  branch  (e.g., 
the  A 1  Ga-site  partial  DOS  has  a  small  peak  in  the  acous¬ 
tic  branch  and  a  large  peak  in  the  optical  branch).  By 
summing  the  partial  densities  for  both  a  P-site  origin  and 
a  Ga-site  origin  we  obtain  the  entire  DOS, 

Im(TrG?)-— Im(TrGoi)  =  5Z)(cl>2)  ,  (13a) 

rr  ir 

where 

TrG=[GM,)  +  2G(£)  +  3G<7'I)  +  3G(7'?) 

+  3G(r?)  +  3G(Tc2)]  .  (13b) 

The  factor  of  5  on  the  right-hand  side  of  Eq.  (13a)  arises 
from  the  presence  of  Five  atoms  in  the  defect  cluster.  The 
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FIG.  I.  Phonon  dispersion  curves  and  DOS  NiE)  of  GaP 
computed  using  the  two-parameter  valence-force  model 
described  in  the  text.  Acoustic  (,4)  and  optical  (O)  branches  are 
labeled  according  to  the  polarization  of  the  phonons  ( T  is  trans¬ 
verse  and  L  is  longitudinal).  Some  branches  may  have  mixed 
polarization. 


ENERGY  (meV) 

FIG  2.  Imaginary  part  of  the  Green’s  function  for  the  per¬ 
fect  crystal  using  a  P  site  for  the  origin  in  space.  At  the  top  is 
shown  the  bulk  DOS.  Elements  of  the  Green’s-function  matrix 
shown  here  are  the  partial  DOS  in  each  of  the  ^4 , ,  E,  Tu  T‘ ,  T*. 
and  T\  CC’s. 
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DOS  obtained  from  Eq.  (13)  can  be  compared  with  that 
directly  computed  from 

D({i):)  =  -77  S  -)  .  (14) 

N  “  «•  k 

>1.  k 

In  our  calculations,  the  D(coz)  computed  in  these  two  ways 
agree  to  an  absolute  accuracy  of  0.004  (for  an  energy  scale 
corresponding  to  rom„=l),  which  provides  an  important 
check  on  the  computation  of  the  Green’s  functions. 


C.  Defect  vibrational  modes 


The  vibrational  modes  of  a  defect  are  characterized  by  a 
peak  in  the  LDOS.  This  peak  is  the  center  of  a  band  of 
vibrational  modes  in  which  the  defect  vibrates  with  signi¬ 
ficantly  greater  amplitude  than  those  atoms  far  removed 
from  it.  The  defect  we  are  considering  consists  of  an  im¬ 
purity  atom  substitutional  for  the  P  atom.  Our  results  are 
for  oxygen  with  a  mass-defect  parameter  of 


Am  m0-mP 
m  m  P 


-0.4834  . 


(15) 


The  springs  which  connect  the  impurity  to  its  nearest 
neighbors  are  described  by  bond-stretching  f\  and  bond- 
bending  f\  interactions.  These  interactions  are  allowed  to 
differ  from  those  of  the  bulk  and  they  are  varied  by  the 
same  fractional  amount.  Thus  the  defect-perturbation 
matrix  depends  on  only  one  parameter,  A f  //,  the  frac¬ 
tional  change  in  spring  constants  for  the  defect  relative  to 
the  bulk  material. 


A f  __  f  i  ~ / 1  fi  —fi 
f  ~  ft  fi 


(16) 


where  f\  and  /2  are  the  bulk  spring  constants.  Negative 
values  of  Af/f  refer  to  a  weakly  bonded  defect  and  posi¬ 
tive  values  refer  to  a  strongly  bonded  defect.  By  defini¬ 
tion  A f/f>  —1.  It  seems  physically  unlikely  that  a  sub¬ 
stitutional  impurity  will  form  bonds  which  are  very  much 
stronger  than  those  of  the  bulk  material.  Therefore,  we 
will  concentrate  our  attention  on  the  range  —  1  <  Af  / 
/<  I- 

For  given  values  of  defect  spring  constants,  we  compute 
the  LDOS  for  each  type  of  vibrational  mode  (Alt  E,  T\, 
and  r2).  The  LDOS  reflects  some  details  of  the  bulk  vi¬ 
brational  modes  and  also  has  peaks  for  each  of  the  defect 
modes.  For  defect-perturbation  values  in  the  range 
—  1  <  Af/f  <  1  we  find  defect  modes  of  A  |  and  T2  sym¬ 
metry  only.  The  A  i  (breathing)  mode  shown  in  Fig.  3  is  a 
resonance  in  the  acoustic  branch.  For  Af /f =0  the  A\ 
LDOS  curve  at  the  top  of  the  figure  is  identical  to  the  A  t 
Green’s  function  for  the  perfect  crystal.  As  Af/f  is  re¬ 
duced  a  resonance  mode  appears  and  moves  to  lower  ener¬ 
gies.  This  resonance  is  rather  broad,  with  a  width  of 
about  5  meV.  The  sum  of  LDOS  for  the  T\,  T\  and  T\ 
modes  is  shown  in  Fig.  4.  These  modes  consist  of  vibra¬ 
tions  of  the  oxygen  atom  itself  along  with  some  response 
of  the  rest  of  the  lattice.  For  Af  /f  ■=■  0  we  obtain  a  strict¬ 
ly  localized  mode  at  56.3  meV.  As  Af/f  is  reduced  this 
r2  mode  falls  in  energy,  entering  the  optical  branch  as  a 
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FIG.  3.  LDOS  of  4 1  modes  for  an  O  atom  on  a  P  site  in 
GaP.  Defect  force  constants  Af //  are  varied  from  Af  / f  =0 
(defect  springs  identical  to  bulk  springs)  to  Af/f  —  —  1  (defect 
springs  to  zero  strength).  Resonance  m  the  acoustic  branch  ap¬ 
pears  as  a  peak  in  the  LDOS.  This  resonance  moves  to  lower 
energies  as  the  defect  force  constant  is  reduced. 

resonance,  becoming  strictly  localized  in  the  acoustic- 
optical  gap,  and  finally  entering  the' acoustic  branch  as  a 
rather  sharp  resonance.  As  Af/f-*—  1,  the  T2  mode  ap¬ 
proaches  a  5  function  at  zero  energy  corresponding  to 
motion  of  the  oxygen  atom  alone  in  the  absence  of  any 
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FIG.  4.  LDOS  of  T:  modes  for  an  O  atom  on  a  P  site  in 
GaP.  For  a  defect  force  constant  of  Af  //  =0  a  strictly  local¬ 
ized  mode  at  56  meV  is  shown  as  a  6  function  in  the  LDOS.  As 
Af/f  is  reduced  this  mode  moves  to  lower  energies,  entering  the 
optical  branch  as  a  resonance,  becoming  strictly  localized  in  the 
acoustic-optical  gap,  and  entering  the  acoustic  branch  as  a  sharp 
resonance. 
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restoring  force.  The  T 2  resonances  in  optical  and  acoustic 
branches  are  relatively  sharp  with  a  width  of  less  than  2 
meV.  The  71-defect  modes  generally  consist  of  at  least 
80%  O  ( 71)  motion.  An  exception  to  this  is  the  strictly 
localized  mode  in  the  energy  range  25—30  meV  which  is 
typically  0.9071  4- 0.0871+0.0271  (presumably  this 
mode  is  comprised  largely  of  bulk  phonons  from  the  very 
top  of  the  LA  branch  which  include  very  little  P-atom 
motion). 

The  results  from  the  LDOS  calculations  are  summa¬ 
rized  in  Fig.  5.  Here  we  plot  the  energies  of  the  defect 
modes  versus  the  defect-spring-constant  parameter.  On 
the  left-hand  side  of  the  figure  a  bulk  DOS  is  shown  for 
reference.  The  solid  lines  in  the  figure  are  the  Green’s- 
function  results.  These  energies  are  defined  as  the  loca¬ 
tion  of  zeros  in  the  real  part  of  the  eigenvalues  of  the  ma¬ 
trix  ( 1  —  G°8L)  [Eq.  (10)].  For  A///  <0  we  see  the  A , 
and  T 2  modes  discussed  above.  For  A///  >0  another  T2 
mode  consisting  of  about  70%  7**  motion  appears  as  a 
resonance  in  the  acoustic  branch.  For  —  1  <  A///  <  1,  no 
defect  modes  of  E  or  T{  symmetry  occur.  Physically,  this 
means  that  E  and  T \  vibrations  of  the  defect  are  strongly 
coupled  to  the  bulk  crystal  so  that  these  types  of  motion 
are  not  localized  around  the  defect.  The  defect  perturba¬ 
tion  for  E  modes  contains  only  bond-bending  interactions 
which  are  relatively  weak  and  produce  defect  modes  only 
for  A///  »1.  The  defect  perturbation  for  T{  modes  is 
identically  zero  as  a  result  of  rotational  invariance  and 
consequently  no  T ,  defect  modes  exist.  The  defect  pertur¬ 
bation  for  A  and  T 2  modes  contains  bond-stretching  and 
mass-defect  terms  which  are  relatively  large  and  produce 
the  modes  shown  here. 


FIG.  5.  Defect  vibrational  inodes  of  oxygen  on  a  P  site  in 
GaP.  Energy  of  the  modes  is  plotted  vs  the  defect  force  con¬ 
stant  A///.  Heavy  solid  lines  show  the  solutions  from  the 
Green 's-function  calculations.  Dashed  and  dotted  lines  give  the 
vibrational  modes  of  an  OGa«  molecule  embedded  in  an  immov¬ 
able  lattice.  E  and  T\  molecular  modes  are  split  by  an  amount 
too  small  to  show  in  the  figure.  Bulk  phonon  DOS  is  shown  on 
the  left-hand  side  of  the  figure. 


The  energy  of  some  Green’s-function  modes  are  quite 
close  to  those  calculated  from  a  simple  molecular  model. 
In  Fig.  5  we  show  the  modes  for  a  OGa4  molecule  imbed¬ 
ded  in  an  immovable  lattice  according  to  the  formulas 
given  in  Appendix  C.  From  the  figure  it  is  apparent  that 
the  A  |  molecular  mode  agrees  very  roughly  with  the  A , 
Green’s-function  mode.  Three  7%  molecular  modes  exist. 
Those  modes  near  25  meV  are  mainly  71-type  (bending 
motion  of  the  Ga  atoms  perpendicular  to  the  O— Ga 
bonds),  those  modes  near  15  meV  are  mainly  T2  type 
(stretching  motion  of  the  Ga  atoms  parallel  to  the  O— Ga 
bonds),  and  the  modes  which  vary  quadratically  with 
A///  are  mainly  71  type  (O-atom  motion).  The  modes 
have  mixed  nature  near  the  crossing  points  of  the 
branches.  Some  of  the  molecular  results  agree  closely 
with  the  Green's-function  modes,  but  other  molecular 
modes  appear  where  there  are  no  Green’s-function  modes. 
In  the  latter  case  the  vibration  of  the  defect  is  strongly 
coupled  to  the  bulk  modes  so  that  a  peak  in  the  LDOS 
does  not  appear.  This  illustrates  the  major  deficiency  of 
the  molecular-model  results — the  number  of  modes  is 
determined  by  the  size  of  the  cluster.  The  Green’s- 
function  method  allows  us  to  couple  these  modes  to  the 
rest  of  the  crystal  in  order  to  see  if  they  remain  somewhat 
localized  around  the  defect.  The  real  importance  of  the 
molecular-model  calculations  is  that  they  provide  a  good 
check  on  the  Green’s  function  results,  since  in  certain  lim¬ 
its  the  results  from  both  methods  agree  quite  well. 

IV.  DISCUSSION 

The  vibrational  states  of  the  neutral  oxygen  defect  O0 
are  seen  in  the  electron-capture  luminescence  of  Dean  and 
Henry.2  In  this  spectrum  the  zero-phonon  line  is  forbid¬ 
den  so  that  the  strong  one-phonon  replicas  must  involve 
T,  phonons.10  As  shown  in  Fig.  2  of  Dean  and  Henry's 
work  vibrational  modes  involving  oxygen  motion  are  seen 
at  24.7  and  28.4  meV  (resonances  in  the  acoustic  branch) 
and  various  phonon  replicas  are  also  seen  in  the  optical 
branch.  Since  oxygen  has  a  smaller  mass  than  phos¬ 
phorous,  it  is  necessary  to  considerably  reduce  the  O— Ga 
force  constants  relative  to  the  bulk  Ga— P  force  constants 
to  produce  such  low-energy  modes.5,6  Comparing  the  ex¬ 
perimental  results  with  the  theoretical  71  modes  shown  in 
Fig.  4  we  see  that  it  is  necessary  to  reduce  the  strength  of 
the  defect  force  constants  to  A///  <  —0.8  to  produce  the 
resonance  in  the  acoustic  branch.  In  Fig.  6  we  show  the 
LDOS  for  71  modes  using  A///  =  —0.85,  and  mass  de¬ 
fects  corresponding  to  l60  and  lsO.  The  theory  predicts 
one  71  resonant  mode  centered  at  21.9  meV  for  l60. 
Dean  and  Henry2  report  a  value  of  —1.6  meV  for  the 
l60— »l80  isotope  shift  of  the  24.7  meV  peak.  As  shown 
in  Fig.  6  we  calculate  a  shift  of  —1.14  meV  for  this  reso¬ 
nance  peak,  which  we  consider  to  be  within  experimental 
uncertainty  of  the  observed  -1.6  meV  considering  the 
limited  resolution  of  the  ,60Ioc  peak  in  the  data.2  Thus  we 
identify  our  computed  21.9-meV  resonance  with  the  ob¬ 
served  24.7-meV  mode.  Our  computed  isotope  shift  is 
83%  of  what  would  be  expected  for  motion  of  the  O  atom 
alone.  Alternatively,  the  LDOS  gives  the  composition  of 
the  71  resonance  mode  to  be  0.0571  +0.8771  +0.877"; 
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FIG.  6.  LDOS  for  T2  modes  of  “O  (solid  line)  and  “O 
(dashed  line)  defects  using  a  defect  force  constant  of 
A/  //  =  —0.85. 


According  to  our  calculations  the  first-shell  Ga  atoms  al¬ 
ways  make  a  significant  contribution  to  the  resonance. 
Motion  of  the  O  atom  alone  is  not  a  normal  mode  of  the 
system  (except  for  A///=  —  1 ). 

The  vibrational  mode  observed2  at  28.4  meV  is  not  well 
described  by  the  present  theory.  However,  we  suggest  that 
this  feature  is  partly  due  to  a  peak  in  the  bulk  T2  DOS 
(LA  phonons)  which  adds  structure  to  the  LDOS.  This 
feature  does  not  appear  in  Fig.  6  because  of  the  simplified 
model  used  for  the  bulk  phonons,  but  better  models8  20  do 
indeed  show  a  strong  peak  near  27  meV.  In  the  optical 
branch  the  experiments  observe  a  number  of  vibrational 
modes.  In  Fig.  6  we  see  some  small  features  in  the  optical 
branch,  but  these  are  associated  only  with  the  bulk  DOS 
and  are  not  localized  near  the  defect.  However,  as  dis¬ 
cussed  above,  the  LO  branch  of  our  theoretical  phonon 
dispersion  curves  is  quite  crude  and  more  sophisticated 
models  may  produce  additional  defect  modes  there. 

Phonons  of  energy  19.5  and  47  meV  have  been  observed 
for  the  0+  defect  in  the  donor-acceptor— pair  (DAP) 
luminescence  of  Dean,  Henry,  and  Frosch.1  Morgan10  ar¬ 
gues  that  these  phonons  are  A  | -type  modes  since  the 
zero-phonon  transition  is  allowed.  We  note  that  phonons 
involving  colinear  motion  of  the  oxygen  and  acceptor 
atoms  are  in  fact  allowed  in  the  DAP  spectra  (these  modes 
are  T2  type  with  respect  to  the  Td  symmetry  of  the  isolat¬ 
ed  O  atom).  This  fact  led  us  originally  to  conclude  that 
the  47-meV  phonon  was  a  T2  mode  involving  resonant  0 
motion.1*  However,  calculations  show  that  the  intensity 
of  a  strictly  localized  7\  phonon  replica  in  a  DAP  spec¬ 
trum  is  very  weak  for  distant  DAP’s.  This  probably  is 
true  for  localized  resonant  modes  also.  We  conclude  that 
the  47-meV  mode  does  not  involve  localized  O  motion. 

In  photoluminescence  excitation  experiments  on  the  0° 
defect,  Monemar  and  Samuelson3  observed  phonons  of  en¬ 
ergy  19  and  48  meV.  These  energies  equal  those  for  the 
0+  defect  within  an  experimental  uncertainty  of  about 
5%.  Alternatively,  Henry  and  Lang21  propose  a  shift  of 
22%  in  the  phonon  energies  of  the  different  charge  states 
in  order  to  explain  observed  nonradiative  capture  cross 
sections.  These  two  results  are  not  necessarily  inconsistent 
since  the  19-  and  48-meV  modes  may  be  relatively  insensi¬ 
tive  to  force-constant  changes,  whereas  the 


capture— cross-section  data  include  the  T2  modes  which 
are  quite  dependent  on  the  defect  force  constants.  A  good 
method  to  determine  the  difference  between  0+  and  O0 
force  constants  would  be  the  observation  of  the  0+  7\ 
mode,  predicted  by  Baraff  et  a/.6  to  lie  at  32.4  meV.  In 
any  case,  considering  the  accuracy  of  the  computations 
present  here,  the  difference  between  0+  and  O0  force  con¬ 
stants  is  relatively  small  and  we  will  henceforth  use  the 
same  force  constants  (A/// =—0.85  as  determined 
above)  for  both  charge  states. 

Turning  now  to  the  calculation,  for  A/// =  —0.85  we 
find  an  A  \  resonance  at  12.2  meV,  which  we  identify  with 
the  observed  19-meV  phonon.  In  the  present  model  the 
computed  energy  of  this  mode  is  somewhat  low,  but  re¬ 
sults  from  more  complicated  models9  yield  an  A !  phonon 
energy  much  closer  to  that  observed.  The  computed  reso¬ 
nance  has  a  width  of  about  5  meV,  which  is  comparable  to 
the  3-meV  lifetime  broadening  estimated  by  Monemar  and 
Samuelson.3  The  breadth  of  this  resonance  indicates  that 
it  is  not  too  strongly  localized  around  the  defect,  so  that 
the  mode  involves  breathing  motion  of  the  first  Ga  shell 
along  with  significant  contributions  from  further  shells. 
In  the  optical  branch  we  find  no  A ,  resonances  for 
A / // a  —0.85.  There  are  no  A ,  modes  here  because  opti¬ 
cal  phonons  involve  mainly  P-atom  motion  (see  Sec. 
Ill  B),  and  in  the  present  calculation  it  is  very  difficult  to 
localize  this  type  of  motion  since  the  P  atoms  are  in  the 
second  shell  around  the  O  atom  but  the  defect  force  con¬ 
stants  extend  only  to  the  first  shell.  Possibly  an  A ,  mode 
in  the  optical  branch  could  be  produced  by  adjusting 
second-nearest-neighbor  O— P  bonds  or  Ga— P  back  bonds. 
However,  from  the  present  model  we  conclude  that  the 
47-meV  mode  does  not  involve  localized  breathing  motion 
of  the  first  Ga  shell. 

In  summary,  for  a  defect  force  constant  of 
A///=—  0.85  our  theory  finds  an  A ,  resonance  at  12.2 
meV  and  a  7;  resonance  at  21.9  meV,  which  we  identify 
with  the  observed  19.5-  (Ref.  1)  and  24.7-meV  (Ref.  2) 
modes,  respectively.  No  other  defect  modes  are  predicted 
by  the  theory  and  we  conclude  that  other  observed  modes 
are  not  localized  at  the  O  impurity  or  in  the  first  Ga  shell. 
As  discussed  above,  the  identification  by  Morgan10  of  the 
observed  47-meV  (Refs.  1  and  3)  phonon  as  an  A  ,  mode 
localized  on  P  atoms  is  consistent  with  our  results.  Also 
we  argue  that  the  observed  28  4-meV  peak  (Ref.  2)  may  be 
partly  due  to  a  peak  in  the  LA-phonon  bulk  DOS.  Simi¬ 
larly,  the  phonons  labeled  A  —C  in  the  capture  lumines¬ 
cence  (Ref.  2)  may  also  be  due  to  peaks  in  the  bulk  DOS 
or  they  may  be  some  type  of  delocalized  vibrational  modes 
of  the  defect. 

V.  CONCLUSIONS 

In  this  paper  we  theoretically  analyze  the  vibrational 
modes  of  substitutional  oxygen  in  GaP  using  nearest- 
neighbor  interactions  only.  Our  theory  predicts  the  pres¬ 
ence  of  basically  one  mode  of  A ,  symmetry  and  one  mode 
of  r,  symmetry  which  are  localized  (strictly  localized  or 
resonant)  around  the  O  impurity.  We  identify  these  local¬ 
ized  modes  with  phonons  observed  in  the  optical  spectra 
of  the  0°  and  O*'  defects,  and  we  argue  that  other  pho- 
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oons  observed  in  these  spectra  are  relatively  delocalized. 
By  the  term  delocalized,  we  mean  a  bulk  phonon  or  a  de¬ 
fect  mode  which  is  not  significantly  localized  at  the  O  im¬ 
purity  or  in  the  first  shell  of  Ga  atoms.  We  should  em¬ 
phasize  that  the  number  of  localized  modes  predicted  by 
the  present  theory  is  largely  determined  by  the  model  used 
for  the  defect  perturbation,  i.e.,  relaxation  of  nearest- 
neighbor  O— Ga  force  constants  only.  More  sophisticated 
defect  models  do  produce  additional  localized  modes,  and 
these  modes  will  be  discussed  in  a  forthcoming  publica¬ 
tion.9  However,  since  the  defect  perturbation  is  basically 
unknown,  we  feel  that  even  the  simple  results  presented 


here  provide  new  insight  concerning  the  nature  of  the  vi¬ 
brational  modes  of  oxygen  in  GaP. 
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APPENDIX  A:  COLLECTIVE  COORDINATES  FOR  A  FIVE-ATOM  TETRAHEDRAL  CLUSTER 

The  collective  coordinates  Q  are  expressed  in  terms  of  the  Cartesian  coordinates  (x,y,z)  of  each  atom.  The  atoms  are 
numbered  0—4  located  along  lattice  directions  (0,0,0),  (1,1,1),  (  —  I,  —  1,1),  (  —  1,1  —  1),  and  (1,  — 1,  — 1),  respectively. 
The  collective  coordinates  form  bases  for.  the  irreducible  representation  of  the  Td  group  as  indicated.  NF  denotes  the 
normalization  factor. 
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APPENDIX  B:  ELEMENTS  OF  THE  DEFECT-PERTURBATION  MATRIX  81  USED 
FOR  THE  GREEN’S-FUNCTION  COMPUTATIONS 


The  defect  consists  of  a  central  O  atom  surrounded  by  four  identical  Ga  atoms.  The  difference  between  the  O— Ga 
force  constants  and  those  of  the  bulk  is  denoted  A/(  for  bond-stretching  and  A/2  for  bond-bending  interactions.  The 
mass  difference  is  denoted  by  Am  =mG— mP.  Bond-bending  interaction  of  strength  fi  are  introduced  between  the  Ga 
atoms  and  \i=  —  A/2/4  by  rotational  invariance.  The  subscripts  a  —c  in  the  7\  matrix  refer  to  the  collective  coordi¬ 
nates  listed  in  Appendix  A.  We  write 


8£(/l,)=— A/,  ,  5A<£)=— ( A/2+3M),  bUT^  )=  —  <A/2  +4M)  ,  8L„( T2 )=  — (A/2  +  j/x>  . 
mo.  mG,  mG,  mG, 

bLbi{T2)=-~( A/,  +  t^>  ,  SZ.«<r2)=—  (|A/,  +  }A/2-<u2Am)  ,  SLa„(T2)  =  bLal>(T2)  =  ~-~~n  , 
mG,  m  p  3m{ 


■G. 


8Lif(  T  2 )  =8Lc4(  r2 )  = 


3mG,mP 


1/2 


A/,  ,  bLac(T2)=bLca(T2)  = 


8 


3  mG,mP 


1/2 


a/2  • 
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APPENDIX  C.  VIBRATION  FREQUENCIES  <o  OF  AN  OGa4  CLUSTER 
EMBEDDED  IN  AN  IMMOVABLE  LATTICE 

Nearest-neighbor  bond-stretching  and  bond-bending  force  constants  /',  and  f\,  respectively,  refer  to  O— Ga  bonds  and 
/i  and  fi  refer  to  those  bonds  connecting  the  Ga  atoms  to  the  immovable  lattice.  Second-nearest-neighbor  bond¬ 
bending  interactions  of  strength  yt  are  introduced  between  the  Ga  atoms  and  =  —  (/':  —f2)/*  by  rotational  invariance. 
The  subscripts  a  -c  in  matrix  Q  refer  to  the  collective  coordinates  listed  in  Appendix  A.  We  have 

|  )= - (/i  +  jf  i  +  7/2 ),  ai1(E)=  (f'i  +  3/i  +  7/1  4-  7/2 ),  <a2(  T\  )  =  (f‘i  +4/i  4-  7/1  +  7/2)  » 

mGa  mGl  mGa 

7*2 )  =  eigenvalues  of  matrix  Q,  C 1 1  =  — - —  (f2  +  yfi  +  7/1  -f  7/2 )  ,  C22  =  — - —  (/ 1  +  7M  +  7/1  +  7/2 )  . 

ma.  mGl 


C]j=  <7/1 +  7/2).  g,2  =c2i  =  —  —  n,  Cn=Cn  — 

mc  an  Gt 


3  mGam0 


/’i.  Ci3=C'3i  — 


3m0,m0 
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The  localized  and  resonant  vibrational  modes  of  a  substitutional  oxygen  defect  in  GaP  have  been 
computed  with  the  use  of  the  Green’s-function  approach.  A  15-parameter  version  of  the 
deformable-dipole  model  was  used  to  describe  the  lattice  dynamics  of  bulk  GaP.  A  defect,  consist¬ 
ing  of  a  mass  perturbation  at  the  substitutional  site  and  a  perturbation  among  the  nearest-  and 
second-nearest-neighbor  short-range  interactions  near  the  oxygen  atom,  was  considered.  The  local 
modes,  labeled  by  symmetry,  are  presented  as  a  function  of  the  interaction  perturbation.  Quantita¬ 
tive  agreement  with  experimental  results  for  0°  and  0+  defects  is  obtained  by  weakening  the  O  (-Ga 
and  -P)  short-range  interactions  to  38%  of  the  bulk  "alues  while  weakening  the  Ga-P  back  bonds. 

With  the  inclusion  of  Coulomb  effects,  this  corresponds  to  an  effective  O-lattice  force  constant  of 
3.5  eV/A2,  which  is  25%  of  the  bulk,  P-lattice  value. 


I.  INTRODUCTION  however,  the  Coulomb  portion  of  the  total  interaction  is 

relatively  small  and  can  be  effectively  contained  in  the 
The  properties  of  the  point  defect  GaP.Op  in  which  ox-  short-range  defect  perturbation  for  sufficiently  localized 
ygen  substituting  for  phosphorous  in  GaP  forms  a  deep  modes. 

donor,  have  been  actively  studied,  often  generating  consid-  The  DBA  and  Green’s-function  theory  are  briefly 
erable  controversy.1-6  In  particular,  numerous  optical  described  in  Sec.  II.  In  that  section  we  also  discuss  the  de¬ 
spectra  have  been  reported.  These  highly  structured  spec-  feet  perturbation  considered  in  this  work,  and  we  derive 

tra  contain  features  due  to  transitions  assisted  by  local  an  exact  expression  for  the  Coulomb  contribution  to  the 

phonons.  By  calculating  the  local  vibrational  modes  of  restoring  force  resulting  from  the  displacement  of  a  single 

the  O  defect,  we  can,  through  comparison  with  the  ob-  atom  in  an  infinite  lattice.  We  find  the  Coulomb  portion 

served  spectra,  directly  identify  mode  symmetries  and  of  the  restoring  force  for  a  P-site  single-atom  displace- 

characterize  the  atomic  motion  near  the  defect.  We  can  ment  to  be  — 17.7%  of  the  bulk  short-range  component, 

also  determine  how  the  short-range  interactions  between  Next,  in  Sec.  Ill,  we  present  our  computed  dispersion 

atoms  near  the  defect  site  are  modified  when  the  defect  is  curves,  density  of  states  (DOS),  and  selected  Green's  func- 

introduced.  In  addition,  assumptions  of  defect-site  tions  for  bulk  GaP.  In  Sec.  Ill  we  also  present  the  results 

geometry  ( Tj  symmetry  and  small  local  lattice  relaxation,  of  our  local-mode  calculations  for  GaP:Op.  (Here,  we 

in  particular)  implicit  in  the  local-mode  calculation  pro-  shall  use  “local  mode”  to  refer  to  both  strictly  localized 

vide  consistency  checks  on  the  interpretation  of  experi-  and  resonant  modes,  as  in  Ref.  7.)  The  interaction  pertur- 

mental  results.  bation  considered  consists  of  varying  certain  O-Ga,  O-P, 

In  this  paper,  we  use  the  Green’s-function  approach  to  and  Ga-P  (back-bond)  force  constants  by  the  respective 

calculate  the  local  phonon  modes  for  oxygen  isotopes  on  a  fractions,  rf,  tj,  and  0. 12 tj,  of  the  corresponding  bulk 

phosphorous  site  in  GaP.  This  work  represents  an  exten-  short-range  values.  Energies  for  local  modes  of  a  given 

sion  of  an  earlier  Green’s-function  calculation7  in  which  symmetry  are  then  presented  as  a  function  of  the  parame- 

bulk  interactions  were  limited  to  nearest-neighbor  (NN)  ter  tj.  The  local  DOS  (LDOS)  and  “eigenvectors”  are 

interactions.  Here,  we  use  the  deformable- bond  approxi-  presented  for  selected  local  modes,  and  a  discussion  of  the 

mation  (DBA)  of  Kune  et  a/.*  as  the  phenomenological  major  features  and  the  choice  of  defect  parameters  is 

description  of  bulk  GaP.  This  15-parameter  model  in-  given. 

eludes  short-range  interactions  out  to  second-nearest  An  extensive  comparison  to  available  photolumines- 
neighbors  (2NN)  and  long-range  Coulomb  effects  arising  cence  data  is  made  in  Sec.  IV.  There,  we  see  that  the 

from  the  electrostatic  interaction  between  induced  di-  electron-capture  luminescence  data  for  the  O0  defect2  is 

poles.*  Our  defect  consists  of  a  mass  perturbation  at  the  well  described  by  our  calculations  for  the  value 

substitutional  site  and  a  perturbation  in  the  short-raige  in-  tj  =  —0.62.  In  particular,  we  identify  the  phonon  replicas 

teractions  among  the  17-atom  cluster  consisting  of  the  O,  at  24.7  and  28.4  meV  with  the  peaks  in  our  T2  LDOS  at 

the  four  neighboring  Ga,  and  the  12  2NN  P  atoms.  We  25.0  and  28.2  meV,  respectively.  Both  peaks  in  the  LDOS 

do  not  explicitly  consider  perturbations  in  the  Coulomb  shift  under  isotopic  substitution  of  l(0  for  160  as  expect- 

interactions  since  their  long  range  would  render  the  ed.  Our  computed  A  t  modes  at  20.4  and  48.2  meV  are  ob- 

Green's-function  approach  intractable.  As  we  shall  see,  served  for  the  charge  states  0+  and  O0. 1,3  In  all  cases  we 
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And  substantial  improvement  in  the  agreement,  both  qual¬ 
itative  and  quantitative,  between  calculated  and  observed 
values  over  that  predicted  by  our  earlier  model.1  Finally, 
the  conclusions  that  can  be  drawn  from  this  work  are 
summarized  in  Sec.  V. 


II.  COMPUTATIONAL  METHOD 
A.  Bulk  phonons 

To  describe  the  lattice  dynamics  of  bulk  GaP  we  have 
used  the  DBA.  This  model  is  a  15-parameter  version  of 
the  deformable-dipole  model  of  Kune  et  al.,  and  has  been 
shown  to  provide  a  good  description  of  several  III-V  com¬ 
pounds.8  The  quantities  of  interest  for  our  calculations 
are  the  force-constant  matrix  <J>  and  the  dynamical  matrix 
Q,  which  is  related  to  <J>  by  the  equation, 
Q=&£~l/2£M.~'/2-  In  this  expression,  M  is  a  diagonal 
matrix  whose  elements  are  equal  to  the  masses  in  the 
problem.  In  the  DBA,  the  force-constant  matrix  has  the 
form 

-£{ti)=(Z+iLt)n-Bia)-'&(Z+!L) .  (l) 

where  <J>tSRI  and  $<el1  correspond,  respectively,  to  the 
short-range  and  electrostatic  terms.  The  short-range 
force-constant  matrix  <£<SRI  contains  .the  non-Coulombic 
interactions  which  in  the  DBA  extend  out  to  2NN’s.  As  a 
consequence  of  the  T}  symmetry  of  GaP,  the  off-diagonal 
Cartesian  blocks  can  be  shown  to  have  the  form' 

ABB 

$lSR,(Ga;P)=$<SRI(P;Ga)=  BAB,  (2a) 

B  B  A 

C|  Z>i  E  | 

<£tSR,(Ga;Ga)=  D,  C,  £,  ,  (2b) 

-E\  -E,  F, 

and 

C2  D2  —Ei 

<|><SR,(P;P)=  Di  C2  -F2  ,  (2c) 

Ei  Ei  Fi 

for  interactions  between  NN’s  [Eq.  (2a)]  and  2NN’s  [Eqs. 
(2b)  and  (2c)],  respectively.  These  Cartesian  blocks  must 
obey  the  transformation  law9 

£af^L,K\L’  ,K')—  .  (3) 

where  £  is  the  point-group  operation  associated  with  the 
space-group  transformation  which  takes  the  site  (/,*)  into 
the  site  (/',*').  In  this  expression,  /,  /',  L,  and  L'  label  the 
N  unit  ceils  in  the  lattice,  *,  K,  and  K'  specify  the  basis 
atom,  and  the  remaining  indices  refer  to  Cartesian  coordi¬ 
nates.  Hence,  from  Eqs.  (2),  (3),  and  the  translation- 
invariance  condition,9 


$</,*;/,*)=-  2'  W 

<j><sR>  for  tj,e  entire  lattice  is  determined.  The  description 
of  the  bulk  short-range  interaction  in  terms  of  these  10  pa¬ 
rameters  (A—  Fz)  is  the  most  general  description  allowed 
by  symmetry,  and  it  physically  corresponds  to  both  cen¬ 
tral  and  noncentral  pair  interactions.8 

In  addition  to  purely  short-range  forces,  the  DBA 
force-constant  matrix  contains  a  contribution  from  classi¬ 
cal  electrostatic  interactions.  Associated  with  each  lattice 
site  is  an  effective  charge  2-  Displacement  of  a  charge 
from  its  equilibrium  site  effectively  results  in  an  electric 
dipole  moment  at  the  site.  Dipoles  may  also  be  induced, 
either  by  the  motion  of  neighboring  atoms  (mechanical 
polarizability  &)  or  by  the  electric  fields  produced  by  oth¬ 
er  dipoles  (electronic  polarizability  a).  The  dipoles,  how¬ 
ever  generated,  produce  electric  fields.  These  fields  exert 
forces  on  the  charges  and  mechanically  induced  dipoles, 
giving  rise  to  the  Coulomb  term  $(el)  of  Eq.  (1).  Five  pa¬ 
rameters  are  used  to  determine  $(d),  giving  a  total  of  15 
for  the  DBA.8  These  15  parameters  have  been  determined 
for  GaP  by  Kune  el  al } 


B.  Green’s  functions  and  group-theoretic  considerations 

The  bulk  Green’s  functions  are  numerically  obtained 
from  the  eigenvalues  and  eigenvectors  of  the  Fourier- 
transformed  dynamical  matrix  in  precisely  the  manner 
described  in  Ref.  7.  From  the  bulk  Green’s  functions  the 
standard  formalism10  is  used  to  obtain  the  local  modes  as 
before.1  The  major  difference  in  the  present  work  lies  in 
the  spatial  extent  of  the  defect  considered.  Here,  the  de¬ 
fect  affects  17  atoms  and,  hence,  51  Cartesian  coordinates. 
Thus,  we  are  required  to  evaluate  a  51x51  block  of  the 
bulk  Green’s-function  matrix. 10 

We  shall  consider  substitutional  defects  which  retain 
the  full  point-group  symmetry  (Td)  of  the  lattice.  Ac¬ 
cordingly,  it  is  natural  and  convenient  to  express  the 
Green’s-function  matrix  in  a  basis  of  “collective  coordi¬ 
nates”  which  transform  according  to  the  irreducible  repre¬ 
sentations  of  Tj.  It  can  be  readily  shown  that  the  51- 
dimensional  reducible  representation  T  of  Td  given  in 
terms  of  Cartesian  coordinates  can  be  reduced  to 

r=3^,©^2®4£®5r,®8r2 .  (5) 


The  collective  coordinates  Q  transforming  according  to 
the  irreducible  representations  (A)—  T2),  are  related  to  the 
original  Cartesian  coordinates  by  an  orthogonal  transfor¬ 
mation.  The  collective  coordinates  have  already  been 
determined  for  the  five-atom  cluster.1  The  additional  col¬ 
lective  coordinates  needed  for  the  12-atom,  2NN  shell  are 
given  in  the  Appendix. 

The  advantage  of  the  collective-coordinate  transforma¬ 
tion  is  that  matrices  reflecting  the  Td  symmetry  are 
block-diagonalized  in  the  form1'9 


6c  xwO(J) 
(jKr)°#ivu<rp  * 


where  (s)  and  (r)  label  the  irreducible  representation.  The 
block  diagonalization  allows  us  to  consider  each  irreduci¬ 
ble  representation  separately,  providing  a  classification 
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scheme  for  the  defect  vibrational  motion  according  to 
symmetry.  For  a  given  irreducible  representation,  p,  and  v 
label  the  dimension,  and  a  and  p  label  the  occurrences 
(e.g.,  *.2,3  and  =  1, ....  8  for  r2).  The  form  of  the 

block  diagonalization  also  restricts  the  number  of  indepen¬ 
dent  matrix  elements  on  the  basis  of  symmetry.  With  the 
use  of  Eq.  (6),  the  number  of  independent  Green’s- 
function  matrix  elements  is  reduced  to  68. 

C.  Defect  perturbation 


tomatically  according  to  rotation  invariance. 

The  dcfect-pcrturbation  matrix  is  now 

where  8M  contains  the  on-site  mass  perturbation.  (Note 
that  we  have  put  the  explicit  mass  dependence  in  64  rath¬ 
er  than  2°,  as  in  Ref.  7.)  We  numerically  construct  64  in 
Cartesian  form,  using  Eqs.  (2)— (4),  (7),  and  (8),  and  then 
apply  the  transformation  to  collective  coordinates.  The 
local  modes  are  given  by  the  frequencies  co  which  solve10 


The  defect  perturbation  considered  in  this  work  consists 
of  the  mass  substitution  m?— *m0  at  the  defect  site  and  a 
perturbation  of  the  short-range  force  constants  among  the 
substituent  and  its  NN’s  and  2NN’s.  We  wish  to  consider 
a  general  force-constant  perturbation  5$  having  the  block 
form  of  Eq.  (2)  and  consistent  with  Td  symmetry.  Let  us 
label  the  O-Ga  perturbation  parameters  &4  and  6 5,  the 
O-P  parameters  6C, . . . ,  5 F,  the  Ga-P  NN  back-bond  pa¬ 
rameters  6 A'  and  65',  the  intra-Ga-shell  2NN  parameters 
6C', . . . ,  BF’,  and  the  intra-P-shell  2NN  parameters 
&C", . . . ,  5 F".  Now,  not  all  these  force-constant  parame¬ 
ters  may  be  chosen  independently.  Applying  the 
infinitesimal-rotation-invariance  condition,9  we  find  that 


5£  =8£'=S£"=0  , 

(7a) 

8/4' -55' =0  , 

(7b) 

6/1  —55  =  — 2[(8C'  — 8D')+6F']  , 

(7c) 

5F=— 2(5C"— SD")  , 

<7d) 

25F''  =  — (8C  —BD)  . 

<7e) 

Redet[l-£2°(<u2)64(o>2)]=0.  (9) 

We  solve  such  a  determinantal  equation  for  each  irreduci¬ 
ble  representation.  A  useful  quantity  for  characterizing 
the  relative  amount  of  motion  of  a  certain  symmetry  type 
and  at  a  particular  frequency  is  the  LDOS.  The  LDOS  of 
the  perturbed  crystal  is  given  by  the  diagonal  elements  of 
—  Im Q/ir,  where  the  perturbed  crystal  Green’s  function  is 
given  by10 

2  =  (l-£054)-1£° .  (10) 

Strong  resonant  modes  show  up  ideally  as  Loren tzi an 
peaks  in  the  LDOS  while  strictly  localized  modes  appear 
as  6  functions.7  For  representations  which  have  multiple 
occurrences  in  the  space  of  the  defect,  the  LDOS  for  each 
occurrence  provides  a  measure  of  the  fraction  of  total  ki¬ 
netic  energy  in  that  particular  type  of  motion,  and  this 
quantity  is  used  in  the  following  sections  to  describe  the 
composition  of  defect  modes. 

D.  Restoring  force  for  one-atom  displacements 


The  linear  combinations  of  the  form  8/4—85  represent 
NN  bond-bending  spring  constants.7  Similarly,  8 C  —6 D, 
8 F,  and  the  corresponding  primed  expressions  may  be 
viewed  as  anisotropic  bond-bending  springs  between  2NN 
pairs,  while  the  parameters  BE  correspond  to  a  more  com¬ 
plicated  pair  interaction.  Essentially,  Eqs.  (7)  impose  con¬ 
straints  on  the  noncentral  pair  interactions  possible  within 
an  isolated  cluster  having  the  geometry  of  our  17-atom  de¬ 
fect.  The  results  up  to  this  point  (given  our  above  as¬ 
sumptions)  are  general,  leaving  nine  free  quantities  to 
specify  8$. 

We  are  now  ready  to  parametrize  the  interaction  pertur¬ 
bation.  To  reduce  the  number  of  free  parameters,  we  as¬ 
sume  the  same  fractional  variation  among  independent  O- 
Ga,  O-P,  and  Ga-P  parameters,  giving 

BA /A  =BB/B  =tjo,g»  >  (8a) 

8C/C2  =BD/D2  —BF/F2  =t)o,p  >  (8b) 

BA' /A  =  i?G..r  .  (8c) 


In  the  discussion  to  follow,  it  will  be  useful  to  know  the 
effective  spring  constant  acting  on  the  O  atom  under  a 
unit  displacement  from  equilibrium.  This  spring  constant 
is  given  t.  (e<r,=k  ,SR)4-A:(el>,  where 

jfc’(SR,=  _4[G4  +&4)+2(C2+8C)+(F2+5F)]  (11) 

denotes  the  modified  short-range  part,  and  /c(el)  is  the  bulk 
P-site  Coulomb  component.  We  wish  to  compute  fcldl.  In 
the  notation  of  Maradudin,9  we  have  for  the  unit  displace¬ 
ment  Ue(l',K')  =  BagBUBtK’, 

Fa{l,K)=-  2  ^hK-,l\K')ba^lVBK„ 

8,1'.*' 

=  -<t>aa(l,K;l,K)  ,  (12) 

where  <PaaU,K;l,K)  =  k{et)  for  <t>=<t>,el).  Now,  using  the 
identity 

Sir—  T7  2e*Pt  — «k-[x(/,K)  —  x(/',k)]|  (13) 

N  r: 
k 


If  we  make  the  analogous  assumptions, 
bC'/Cx=bD'/Dx=bF'/Fx  and  8C,VC2=8D'VD2(  we 
see  from  Eqs.  (7)  that  the  interaction  perturbation  5$  is 
completely  specified  in  terms  of  the  tj’s.  In  this  way,  we 
specify  directly  the  NN  (tjo.oi).  the  adjustable  part  of  the 
2NN  (t7o,p).  and  the  bond-stretching  part, 
8/4 '+285'  =  384',  of  the  back-bond  interactions  (pG,p). 
Other  interactions  within  the  cluster  are  then  modified  au- 


in  Eq.  (12)  gives 
F„(/,k)=--^-  2 


X  exp|  — » k-[x(/,*)  —  x(/',*)]| 


=  I  *t>  , 

k 


(14) 
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where  the  inner  sum  has  been  identified  as  the  Fourier- 
transformed  force-constant  matrix.  Recognizing  that  our 
$(*,*|  k)  transforms  as  T2®T2,  we  can  restrict  the  sum¬ 
mation  over  k  in  Eq.  (14)  to  the  irreducible  wedge  W  of 
the  first  Brillouin  zone  as  in  Ref.  7,  obtaining  finally 

*'*"=77  2  h7  2  £>]  ■  (15) 

N  lew  *-> 

This  sum  was  evaluated  numerically  by  a  method  similar 
to  the  numerical  computation  of  Im£°  in  this  and  in  our 
previous  work.7  The  result  is  k(el,=  — 0. 1771c<SRl  for  the 
bulk  P  site. 

III.  RESULTS 
A.  Bulk  phonons 

In  Fig.  1  we  show  the  bulk  dispersion  curves  and  bulk 
DOS  for  GaP  as  computed  from  the  DBA.  The  DBA  re¬ 
sults  constitute  a  substantial  quantitative  improvement 
over  those  of  the  two-parameter  model  used  previously.7 
The  calculated  LOp  and  TOr  energies,  50.0  and  45.3  meV, 
respectively,  are  split  due  to  Coulomb  effects  and  are  in 
agreement  with  measured  values.8  The  calculated  TA* 
energy  is  12.8  meV  compared  with  the  actual  value,  13.1 
meV.11  The  bulk  DOS  presented  here  is  highly  structured 
compared  with  that  of  the  two-parameter  model.7  As  we 
shall  see,  this  structure  gives  rise  to  local-mode  behavior 
not  seen  in  the  simpler  model.  In  the  units  shown  in  the 
figure,  the  computed  area  under  the  bulk-DOS  curve  is 
6.0004,  compared  to  the  exact  theoretical  value,  6. 

B.  Bulk  Green’s  functions 

We  have  computed  the  bulk  Green’s-function  matrix  2° 
for  both  a  P-ccntered  and  a  Ga-centered  17-atom  cluster. 
The  method  of  computation  is  described  in  Ref.  7.  The 
Green’s  functions  were  calculated  as  functions  of  (fio))2  at 
intervals  of  2.5  meV2  and  then  converted  to  functions  of 
E  =i ico  [see  Eq.  (9)  of  Ref.  7],  This  corresponds  to  a  reso¬ 
lution  in  E  of  about  0.5  meV  near  £  =0  and  0.025  meV  at 


FIG.  1.  Bulk-phonon  dispersion  curves  and  DOS  for  GaP 
calculated  from  the  DBA.  The  calculated  LOp,  TOr,  and  TA* 
energies  are  50,  45.3,  and  12.8  meV,  respectively.  The  peaks  in 
the  bulk  DOS  lie  at  11.3,  27,  44.8,  48.8,  and  52.9  meV.  In  (he 
units  shown,  the  DOS  integrates  to  6.0004,  compared  with  the 
exact  value,  6. 


£  =50  meV.  There  arc  68  bulk  Green’s  functions  each 
for  both  the  P-  and  Ga-centered  clusters — too  many  to 
display  individually.  Instead,  for  each  irreducible  repre¬ 
sentation  (s),  wc  sum  the  diagonal  matrix  elements  and 
present  ^|J,(£)  vs  £  in  Fig.  2,  where 

‘0l,'(£)=--TrIm[2O(,,(£)]  .  (16) 

ir 

The  quantity  &,S>(E)  is  the  partial  DOS  (Ref.  7)  for  the 
representation  (s)  and  is  proportional  to  the  mean-squared 
amplitude  of  modes  associated  with  (s)-type  motion  in  the 
17-atom  cluster  at  energy  £.  In  particular,  the  A2  mode 
consists  of  2NN  motion  only.  From  Fig.  2  we  see  that  the 
P-centered,  A2  partial  DOS  has  a  large  amplitude  in  the 
TO  and  small  amplitude  in  the  TA  phonon  branches, 
while  the  opposite  is  true  for  the  Ga-centered  DOS.  This 
is  because  the  P  atoms  undergo  larger  displacements  in 
optical  modes  than  do  the  heavier  Ga  atoms.  For  the  dis¬ 
cussion  to  follow,  the  important  symmetries  are  A  t  and 
T2.  A  |  modes  consist  of  some  combination  of  the  follow¬ 
ing:  “breathing”  motion  of  the  NN  shell,  2NN  breathing, 
and  a  third,  more  complicated  2NN  motion.  T2  modes 
consist  of  defect-site  motion  along  with  collective  motions 
within  the  NN  and  the  2NN  shells.  Only  T2  modes  can 
involve  defect-site  motion.  (See  the  Appendixes  here  and 
in  Ref.  7  for  a  complete  listing  of  the  collective  coordi- 
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energy  (meV) 

FIG.  2.  Partial  DOS  of  bulk  GaP  for  the  representations. 
A\—T2.  Green’s  functions  for  P-  and  Ga-centered  clusters  are 
shown.  The  A  2  partial  DOS  involves  2NN  motion  only  and  for 
the  P-centered  (Ga-centered)  cluster  is  larger  at  optical  (acoustic) 
frequencies. 
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nates  of  the  17-atom  cluster.) 

Figure  2  gives  the  decomposition  of  the  bulk  DOS  into 
the  various  representations  of  Td  appropriate  for  a  point 
defect.  The  bulk  DOS  can  be  recovered  from  the  various 
partial  densities  by  performing  the  sum 

Tl 

2  d(J,[^<P,,(£)+^oi(£)]  =  17iV(£)  ,  (17) 

m-a, 

where  d{t)  is  the  dimensionality  of  the  representation  (s). 
The  renormalization  factor  of  17  on  the  right-hand  side  of 
Eq.  (17)  occurs  because,  in  summing  over  both  P-  and 
Ga-centered  clusters,  the  two-atom  basis  of  GaP  is  includ¬ 
ed  17  times. 

C.  Local  modes 

Let  us  examine  the  local  phonon  modes  of  GaP  Op  as  a 
function  of  the  interaction  perturbation.  In  Sec.  II C  we 
have  parametrized  the  interaction  perturbation  for 
GaP:Op  in  terms  of  t/0,g.>  Vo.p<  and  t;g,.p.  Now,  we  feel 
physically  motivated  to  require  a  further  constraint  be¬ 
tween  77o,Ga  and  Vo.p-  It  is  natural  to  expect  the  O-P  in¬ 
teraction  to  vanish  when  the  O-Ga  interaction  vanishes, 
and  the  O-P  interaction  to  have  the  bulk  (P-P)  value  when 
the  O-Ga  interaction  has  the  bulk  (P-Ga)  value.  In  this 
way  we  are  led  to  consider  defects  for  which 
Vo.Gt  —  Vo.p^V-  In  contrast,  we  expect  the  relation  be¬ 
tween  the  back  bonds  (r/Gl  P)  and  our  newly  defined  rj  to 
.  depend  on  details  of  the  defect  chemistry  well  beyond  the 
scope  of  this  paper.  In  practice,  we  have  found  that  the 
effect  of  the  parameter  t/g^p  is  to  “turn  on”  certain 
2NN-shell  local  modes  while  not  appreciably  affecting 
those  low-shell  modes  which  have  a  strong  dependence  on 
V-  Therefore,  for  the  purpose  of  presentation,  we  choose 
(see  Sec.  IV)  t/Ga  P=0. 1277  and  present  our  results  as  a 
function  of  the  single  parameter,  17. 

The  local  modes  of  GaP:Op  as  a  function  of  the  per¬ 
turbing  force-constant  parameter  77  are  shown  in  Fig.  3. 
Only  A 1  and  T 2  modes  satisfying  Eq.  (9)  are  shown  here. 
Modes  of  A 2,  £,  and  T ,  symmetry  also  occur,  but  are  om¬ 
itted  for  clarity,  77  ranges  from  —  l,  corresponding  to  re¬ 
moval  of  the  short-range  O-Ga  and  O-P  interactions, 
through  0,  where  all  interactions  have  their  bulk  values,  to 
+  1.  The  effective  O-site  restoring  force,  kiem 
=(  1  +77— 0. 177)/(  1  — 0. 177),  normalized  to  the  bulk  P- 
site  value,  is  indicated  along  the  top  of  the  figure  while 
the  bulk  DOS  adjoins  on  the  left.  Note  that  the  shaded  re¬ 
gion  below  the  value  fc(eff)=0  is  unphysical. 

The  important  features  of  Fig.  3  are  the  steep  T2 
branch  in  the  energy  range  near  12-40  meV,  the  lowest- 
lying  T j  branches  (around  10  meV),  the  low-energy  A  | 
branch  (20—35  meV),  and  the  high-energy  Ax  branch 
(45—55  meV).  The  steep  T2  branch  was  seen  in  our  ear¬ 
lier  model7  and  consists  mainly  of  O-atom  motion.  Also 
seen  in  our  earlier  work  was  the  low-energy  A ,  branch.  In 
the  present  work,  this  A 1  branch,  consisting  of  over  95% 
NN  breathing,  has  been  pulled  up  to  higher  energies 
(77  <0)  in  better  agreement  with  experiment,  as  we  shall 
see  below.  Two  significant  features  appear  in  the  present 
work  not  obtained  previously.  One  is  the  high-energy  A , 


0-site  restoring  force 
0.25 
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per tur Ding  force  constant,  17 


FIG.  3.  A 1  and  T,  local-mode  energies  for  GaPOp  as  a  func¬ 
tion  of  the  short-range  force-constant  perturbation  77  discussed 
in  the  text.  For  ij=  —  0.62  the  A  rmode  energies  are  20.4  and 
48.2  meV  and  the  7'2-mode  energies  ire  7.3,  28.2,  and  43.1  meV. 
The  shaded  region  is  unphysical.  The  DOS  for  bulk  GaP  ad¬ 
joins  on  the  left. 

branch,  which  appears  only  upon  perturbation  of  the 
Ga-P  back-bond  force  constants  and  involves  mostly 
2NN-breathing  motion.  The  other  is  the  low-energy  7\ 
branch.  This  T2  branch  is  primarily  composed  of  Ga- 
shell,  bond-stretching  motion  (£>  6  meV)  and  seems  to 
appear  for  77  <  0  only  when  the  bulk  model  includes  2NN 
interactions. 

In  Sec.  IV  we  shall  see  that  the  O0  defect  is  well 
described  by  the  perturbation  77=— 0.62.  This  corre¬ 
sponds  to  a  kietT)  of  25%  of  the  bulk  value.  The  local¬ 
mode  energies  at  77=  —0.62  are  indicated  in  the  caption  of 
Fig.  3.  Next,  in  Fig.  4,  we  compare  the  A ,  bulk  and  de¬ 
fect  partial  DOS  for  77=— 0.62.  These  are  decomposed 
into  their  Ga-  and  P-shell  parts,  showing  the  relative  ener¬ 
gy  of  motion  associated  with  each  shell  at  a  given  mode 
“energy,”  £  =ha).  From  the  top  curve  in  Fig.  4  we  see 
resonant  modes  at  20.4  and  48.2  meV.  The  20.4-meV 
mode  is  a  rather  broad  resonance  of  width  —  5  meV  and 
involves  mainly  NN  Ga  motion  (96%)  as  mentioned 
above.  The  resonant  mode  at  48.2  meV  is  sharper  ( —  0.5 
meV)  and  consists  of  motion  (97%)  of  the  less  massive  P 
atoms.  It  is  evident  from  the  figure  that  both  resonant 
modes  occur  at  frequencies  appreciably  different  from  the 
bulk  peaks,  and  hence  are  characteristic  of  the  defect.  In 
contrast,  local  modes  confined  near  singular  features  in 
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FIG.  4.  A  i  LDOS  for  GaP:Of  decomposed  into  Ga-  and  P- 
shell  components.  The  corresponding  bulk  partial  DOS's  are  in¬ 
cluded  for  comparison.  The  20.4-meV  resonance  involves  most¬ 
ly  Ga-shell  motion,  while  the  resonance  at  48.2  meV  consists 
mainly  of  P-shell  motion.  All  curves  have  the  same  scale. 


the  bulk  partial  DOS  may  also,  in  general,  occur.  Such 
modes  evidently  couple  strongly  to  bulk  modes  and  are 
thereby  more  “bulklike,”  i.e.,  more  spatially  delocalized  in 
character. 

We  show  the  T2  LDOS  for  GaP:Op  at  77=— 0.62  in 
Fig.  3,  once  again  decomposed  into  contributions  from  O, 
Ga,  and  P  shells.  The  mode  energies  which  solve  Eq.  (9) 
are  indicated  on  the  top  curve.  Twin  resonance  peaks  at 
25.0  and  28.2  meV  are  seen  in  the  O-site  LDOS,  occurring 
just  above  and  just  below  the  27-meV  LA  peak  in  the  bulk 
DOS.  The  28.2-meV  mode  consists  of  86%  0-,  8%  Ga-, 
and  6%  P-shell  motion,  and  the  25.0-meV  mode  consists 
of  70%  O-,  20%  Ga-,  and  10%  P-shell  motion,  within  the 
17-atom  cluster.  We  note  that  at  17=— 0.62  only  the 
mode  at  28.2  meV  is  a  solution  of  Eq.  (9).  The  twin-peak 
structure  appears  to  be  the  result  of  coupling  between 
isolated-defect  and  nearby  bulk  LA-phonon  modes.  The 
27-meV  LA  peak  in  the  bulk  DOS  involves  very  little  P- 
sublattice  motion.  Thus,  a  mode  involving  large  O  motion 
cannot  exist  at  27  meV,  and  when  r\— » —0.62  the  O  mode 
splits  into  two  modes.  Both  O-site  peaks  occur  at  energies 
different  from  bulk  peaks  and  are  well  localized  near  the 
O  atom,  and  upon  the  substitution  of  lsO  for  160,  both 
peak  positions  undergo  a  shift  in  energy.  The  25.0-  and 
28.2-meV  resonances  are,  therefore,  local  rather  than  bulk¬ 
like  in  character.  In  contrast,  the  weak  mode  at  43.1  meV 
is  pinned  to  a  van  Hove  singularity  in  the  bulk  DOS  (see 
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FIG.  5.  r2  LDOS  of  Ga:0P  for  isotopes  ,sO  and  "O  decom¬ 
posed  into  O-atom,  NN-shell,  and  2NN-shell  parts.  The  twin 
resonances  centered  near  25  meV  primarily  involve  O-atom 
motion  and  undergo  isotope  shifts.  Other  resonant  modes  occur 
at  7.3  meV  (Ga  shell)  and  43.1  meV  (bulklike  P  shell).  Al! 
curves  have  the  same  scale. 


Fig.  3)  and  is  evidently  somewhat  more  bulklike.  Finally, 
the  sharp  resonance  at  7.3  meV  primarily  involves  Ga- 
shell  motion  (94%)  along  with  small  P  (5%)  and  0(1% 
contributions. 

TV.  DISCUSSION 

A.  Intensity  of  phonon  replicas 

The  intensity  of  phonon-assisted  optical  transitions  at 
defects  have  been  thoroughly  discussed  by  Rebane.  - 
Selection  rules  for  determining  which  transitions  are  al¬ 
lowed  are  similar  to  those  used  for  molecules,  as  discussed 
by  Herzberg.13  Allowed  transitions  are  those  for  which 
the  dipole  matrix  element 

{'Vm(r,R)Q>mj(R)  |  D,(r)  |  V,(r,R)<l>, ,(/?)> 

is  nonzero.  Here,  'P  is  the  electronic  wave  function,  <J>  is 
the  vibrational  wave  function,  /  and  m  label  electromc 
states,  i  and  j  label  vibrational  states,  r  is  the  set  of  all 
electron  coordinates,  R  is  the  set  of  all  nuclear  coordi 
nates,  and  D,(r)=  —  e  ri  *s  the  electronic  dipole  mo 
ment  operator.  In  place  of  the  nuclear  coordinates  R,  w< 
introduce  new  coordinates  Q;  which  are  the  normal  coor 
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dinates  of  the  vibrational  problem.  These  normal  coordi¬ 
nates  may  be  chosen  to  form  bases  for  irreducible  repre¬ 
sentations  T(  Qit )  of  the  symmetry  group  of  the  imperfect 
crystal,  where  Hx)  refers  to  the  representation  according 
to  which  x  transforms.  We  assume  that  the  same  set  of 
normal  coordinates  may  be  used  for  both  the  initial  and  fi¬ 
nal  electronic  states  (though  the  equilibrium  values  of 
these  coordinates  will  depend  on  the  electronic  state).  The 
total  vibrational  wave  function  can  now  be  written  as 

*M..i=n  *.,<&>•  (i8> 

i 

where  i,  labels  the  occupation  of  the  sth  normal  coordi¬ 
nate.  The  harmonic-oscillator  functions  <j>  are  even  with 
respect  to  their  arguments  for  even  occupation,  and  odd 
for  odd  occupation.  By  inspection,  in  the  vibrational 
ground  state  (i,  =0  for  all  s)  the  wave  function  (  is  to¬ 
tally  symmetric.  For  multiple  one-phonon  occupation,  <t> 
transforms  as  fll'.-il  02«).  where  the  product  runs  over 
all  s  such  that  /,  =  1 . 

As  an  example  of  the  above  discussion,  consider  a  de¬ 
fect  with  Td  symmetry,  in  which  case  the  dipole  moment 
transforms  a s  T2.  An  A  |  -*A  \  electronic  transition  is  for¬ 
bidden,  but  if  one  T2  phonon  is  involved,  the  dipole  ma¬ 
trix  element  transforms  as 

A  |  ®A  i ®  T2®A  i  ®  7*2  ~A  |  ® E <8  7" |  ®  7*2  , 

which  is  allowed.  Similarly,  an  At— *T2  electronic  transi¬ 
tion  is  allowed,  and  participation  of  a  single  A  (  phonon  or 
a  single  T2  phonon  is  also  allowed  since 

\®A  i®  T  2®  T2®A  i  and  Ai®Ai®T2®T2®  T2 

both  contain  A\.  The  intensity  of  phonon  replicas  is 
governed  by  the  Franck-Condon  principle.12-13  Essential¬ 


ly,  this  principle  states  that  the  phonon  replicas  which  will 
be  observed  in  a  transition  are  those  which  describe  the 
change  in  equilibrium  position  of  the  atoms  between  the 
initial  and  final  states.  In  other  words,  if  the  atoms  move 
along  some  coordinate  Q  between  their  initial  and  final 
states,  normal  modes  of  the  type  Q  will  be  strongly  excited 
in  the  transition. 

We  can  use  this  principle  to  demonstrate  that  "antisym¬ 
metric,”  i.e.,  not  totally  symmetric,  phonons  cannot  be 
strongly  excited.  If  some  antisymmetric  vibrational  mode 
is  strongly  excited,  this  necessarily  leads  to  a  change  in  the 
symmetry  of  the  system,  and  in  the  new  symmetry  group 
the  vibrational  mode  is  totally  symmetric.  However,  for 
weak  excitation  of  the  antisymmetric  vibrations,  this 
change  in  symmetry  can  be  neglected.  Thus,  in  electroni¬ 
cally  allowed  transitions,  totally  symmetric  vibrations  will 
predominate,  and  in  electronically  forbidden  transitions, 
the  occupation  of  the  dominant  antisymmetric  modes 
changes  by  one. 

B.  Comparison  of  theory  and  experiment 

In  Table  I  we  compare  our  theoretical  results  with  the 
phonon  energies  observed  in  the  optical  spectra  of  the  0+ 
and  O0  defects.  The  apparent  agreement  between  theory 
and  experiment  is  achieved  partly  by  matching  unknown 
theoretical  parameters  with  experiment,  and  partly  by  the 
predictive  powers  of  the  theory  itself.  If  we  view  the  bulk 
phonon  dispersion  curves  as  input  to  the  theory,  then  the 
only  unknown  is  the  interaction  perturbation  within  the 
17-atom  cluster.  We  have  parametrized  the  interaction 
perturbation  [Eqs.  18)]  in  Sec.  II C  and  have  invoked  the 
further  restriction,  rjQ  Gi  =  t7o.p  =  1/  in  Sec.  IIIC,  so  that 
the  total  number  of  unknown  parameters  is  two.  We 
determine  values  for  these  parameters  by  matching  the  en- 


TABLE  I.  Energies  (meV)  of  the  phonons  associated  with  the  O0  and  0+  defect  in  GaP  l60—  "O 
isotope  shifts  (meV)  are  given  in  parentheses  following  the  energies.  Tentative  identifications  are  indi¬ 
cated  by 


Label 

Expt. 

Theory 

Identification 

Zn-O 

6.0* 

7.3 

defect  ( T 2,  Ga  shell) 

TA 

13.1" 

11.5 

bulk  (TA) 

a,(G  a) 

19.5, c  19* 

20.4 

defect  (A  t,  Ga  shell) 

loc 

24.7(— 1.6)b 

25. 0(  —0.6) 

defect  (T  2,  O  atom) 

loc' 

28.4-0.5)'’ 

28. 2(  —0.7) 

defect  (T;,  O  atom) 

A 

43.0* 

43.1 

defect  IT;,  P  shell  T 

TO 

44. 8b 

44.8 

bulk  (TO) 

e(P) 

45.3 

defect  iE,  P  shell) 

r,(P) 

45.4 

defect  (Ti,  P  shell) 

Cj(P) 

45.8 

defect  (A  2,  P  shell) 

B 

46.1,b  46.4, 46.5' 

48.8 

defect  <r,,  P  shell)? 

a,(P) 

47.5,b  47/  48d 

48.2 

defect  (4 1.  P  shell) 

C 

48.7,b  48. 8‘ 

48.8 

bulk  (LO) 

LOr 

49.8b 

50.0 

bulk  (LOr* 

•Reference  17. 


•Reference  2,  capture-luminescence  spectrum. 
‘Reference  2,  luminescence-excitation  spectrum. 
*  Reference  3 
‘Reference  14. 
r  Reference  1. 
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ergies  of  theoretical  modes  with  those  observed.  The 
theory  then  predicts  a  number  of  other  vibrational  modes 
not  involved  in  the  matching  process  with  which  we  iden¬ 
tify  observed  phonons.  In  this  way  we  determine  the  sym¬ 
metry  and  character  of  almost  all  of  the  phonon  modes  as¬ 
sociated  with  the  0+  and  O0  defects.  Strictly  speaking, 
all  of  the  phonons  observed  in  the  optical  spectra  involve 
some  component  of  defect  motion  (i.e.,  O-atom,  nearest 
Ga-shell,  or  nearest  P-shell  motion).  However,  certain 
phonons  involve  much  more  defect  motion  than  others, 
and  it  is  useful  to  separate  the  modes  into  those  which  are 
associated  largely  with  the  defect  and  those  which  are 
bulklike.  In  Table  I  we  have  classified  the  modes  as  such. 
The  theoretical  defect-mode  energies  are  taken  from  peaks 
in  the  LDOS,  and  the  theoretical  bulk  phonon  energies  are 
taken  from  peaks  in  the  bulk  DOS  or  critical  points  in  the 
dispersion  curves. 

A  set  of  very  interesting  phonon  sidebands  associated 
with  the  neutral  oxygen  defect  O0  occur  in  the  capture 
luminescence  of  Dean  and  Henry.2  Here,  intense  phonon 
resonances  near  25  meV  are  observed  to  involve  signifi¬ 
cant  O  motion  since  they  shift  under  the  isotopic  substitu¬ 
tion,  l60— *l80  and  must  therefore  be  7%  modes.  Since 
the  O  mass  is  less  than  that  of  P,  these  low-energy  reso¬ 
nances  imply  reduced  O-neighbor  force  constants  relative 
to  the  bulk  force  constants.  This  force-constant  reduction 
was  first  estimated  by  Feenstra  and  McGill14  to  be  30%, 
and  later  estimated  by  Baraff,  Kane,  and  Schliiter15  to  be 
20%;  both  estimates  were  based  on  simple  molecular- type 
vibrational  models.  Subsequent  NN  Green’s-function  re¬ 
sults7  yield  an  estimate  of  15%.  In  the  present  model,  if 
we  assume  the  same  fractional  reduction  for  O-Ga  and 
O-P  force  constants,  the  best  value  for  the  purely  short- 
range  force-constant  reduction  is  77=—  0.62.  After  in¬ 
cluding  Coulomb  effects  between  the  O  atom  and  the  rest 
of  the  lattice,  this  corresponds  to  a  total  effective  force 
constant  (k{ttt))  of  3.5  eV/A2,  which  is  25%  of  the  bulk, 
P-site  value.  We  emphasize,  though,  that  kle(f>  includes 
long-range  Coulomb  effects  from  the  entire  lattice  and 
must  not  be  viewed  simply  as  an  equivalent  short-range 
force  constant.  Our  LDOS  for  T2  modes,  summed  over 
the  O  atom,  Ga  shell,  and  P  shell,  is  shown  in  Fig.  5.  The 
theoretical  peak  energies  and  isotope  shifts  listed  in  Table 
I  compare  reasonably  well  with  experiment. 

Having  determined  the  O-neighbor  force  constants,  the 
theory  predicts  an  A  \  mode  at  20.4  meV.  We  identify  this 
mode  with  the  observed  19-meV  phonon  listed  in  Table  I. 
As  discussed  by  Morgan4  and  also  in  Sec.  IV  A,  the  ob¬ 
served  intensities  of  this  19-meV  mode  and  the  prominent 
47-meV  mode  indicate  that  these  modes  have  A  t  symme¬ 
try.  To  produce  a  defect  mode  of  A  |  symmetry  near  47 
meV,  it  is  necessary  to  reduce  the  back-bond  Ga-P  force 
constants  from  their  bulk  values.  We  arbitrarily  choose  to 
reduce  the  Ga-P  force  constants  by  12%  of  the  O- 
neighbor  reduction,  resulting  in  ijai  P= —0.0744  for  the 
Ga-P  interactions  and  yielding  an  A 1  mode  at  48.2  meV. 
Further  reductions  in  the  Ga-P  force  constants  reduce  the 
energy  of  this  mode,  bringing  it  closer  to  the  observed 
value.  However,  considering  the  uncertainties  associated 
with  our  calculated  optical  phonons,  we  feel  that  an  accu¬ 
rate  determination  of  these  back-bond  Ga-P  interactions  is 


outside  the  range  of  our  theory.  Suffice  it  to  say  that  our 
computations  do  yield  an  A ,  mode  in  the  optical  branch 
and  that  the  back-bond  Ga-P  interactions  are  reduced 
from  their  bulk  values. 

A  number  of  phonon  replicas  in  the  optical  branch  are 
seen  in  the  optical  spectra  of  Dean  and  Henry.2  We  iden¬ 
tify  the  symmetry  and  nature  of  these  phonons  as  indicat¬ 
ed  in  Table  I.  Some  of  these  optical-phonon  replicas  are 
weak,  and  considering  the  uncertainty  of  our  optical- 
phonon  computations,  those  identifications  we  feel  should 
be  regarded  as  tentative  are  marked  by  “?.”  The  a^P) 
phonon  has  been  discussed  above.  The  phonon  labeled  B 
we  take  to  be  a  defect  mode  since  its  energy  falls  well 
away  from  peaks  in  the  bulk  DOS.  This  phonon  is  ob¬ 
served  in  capture-luminescence  and  luminescence- 
excitation  spectra,2  and  is  also  seen  in  the  near-neighbor 
(m  =4  and  5)  Zn-O  pair  spectrum14  (labeled  O  in  that 
spectrum).  Considering  the  proximity  of  the  donor- 
acceptor  impurities  in  the  latter  case,  it  seems  that  this 
46.5-meV  mode  corresponds  to  atomic  motion  which  is 
directed  along  the  donor-acceptor  axis,  i.e.,  a  T2  mode 
with  respect  to  the  Tj  symmetry  of  an  isolated  O  impuri¬ 
ty.  The  46.1-meV  B  phonon  in  the  capture-luminescence 
spectrum  most  likely  has  T2  symmetry,4  although  the 
46.4-meV  mode  seen  in  luminescence-excitation  spectra 
could  transform  as  either  /l,  or  Tz,  as  discussed  in  Sec. 
IV  A.  For  now,  we  group  all  three  of  these  modes  togeth¬ 
er  with  the  label  B  in  Table  I,  and  we  tentatively  identify 
them  with  a  computed  7\  resonance  at  48.8  meV  which  is 
strongly  tied  to  a  LO  peak. in  the  bulk  DOS  (see  Figs.  3 
and  5). 

Another  7%  defect  mode  occurs  in  the  theory  at  43.1 
meV,  and  this  mode  we  tentatively  associate  with  the  ob¬ 
served  43.0-meV  phonon  labeled  A.  The  49.8-meV  pho¬ 
non  identified  by  Dean  and  Henry  as  LOr  is  in  agreement 
with  our  computed  50.0-meV  LOr-phonon  energy.  The 
mode  labeled  C  we  associate  with  the  prominent  peak  in 
the  bulk  DOS  occurring  at  48.8  meV,  while  we  associate 
the  observed  44.8-meV  mode  with  the  very  prominent  TO 
peak  in  the  bulk  DOS  occurring  at  44.8  meV.  In  both  of 
these  cases,  the  peaks  in  the  bulk  DOS  do  not  correspond 
to  phonons  at  or  near  critical  points  in  the  Brillouin  zone, 
but  rather  are  associated  with  phonons  in  the  interior  of 
the  zone.  Finally,  with  regard  to  defect  modes  in  the  opti¬ 
cal  branch,  we  note  that  the  present  theory  neglects 
defect-induced  electrostatic  interactions,  and  such  effects 
are  known  to  cause  vibrational  modes  split  off  from  the 
LO  branch.16 

We  now  come  to  the  entry  in  Table  I  labeled  Zn-O. 
This  6.0-meV  mode  is  just  barely  resolved  in  the  NN  Zn- 
O  pair  spectrum  of  Henry,  Dean,  and  Cuthbert17  (a  related 
mode14  is  seen  at  9.2  meV  for  2NN  Zn-O  pairs).  Al¬ 
though  the  computational  results  presented  here  apply  to 
the  isolated  O  defect,  simply  replacing  a  NN  Ga  with  a 
Zn  atom  changes  the  phonon  energies  by  less  than  0.1 
meV  since  Ga  and  Zn  have  almost  the  same  mass.  Fur¬ 
ther  force-constant  variations  of  the  Zn-O,  Zn-P,  or  Zn- 
Ga  interactions  could  affect  the  mode  energies,  but  we  as¬ 
sume  that  these  variations  are  not  significant.  The  pres¬ 
ence  of  the  Zn  atom  will,  however,  drastically  affect  the 
selection  rules  for  observing  phonons.  In  particular,  all 
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modes  which  are  T2  for  isolated  O  are  now  split  to  A  \  ©  E, 
and  strong  phonon  replicas  of  the  allowed  electronic  tran¬ 
sition  transform  as  A  i .  Thus,  the  low-energy  T 2  mode  in 
Fig.  S  should  appear  in  the  Zn-0  spectrum,  and  we  identi¬ 
fy  it  with  the  observed  6.0-meV  mode.  The  theoretical 
“eigenvector”  for  this  mode  is  shown  in  Fig.  6.  We  plot 
the  displacement  eigenvector,  u  =M_I/2|,  satisfying  the 
eigenvalue  equation 

[l-£°(4>2)54<<y2)]£  =  A.(<uJ)£  (19) 

at  the  resonant  frequency,  Au  =  7.3  meV,  for  which 
ReA=0. 18  Now,  the  actual  defect  normal  mode,  given  by 
(1—  G°SL.)£=0,  requires  a  complex  value  of  oi  for  its 
solution.  However,  since  ImA.  is  proportional  to  the  reso¬ 
nance  width,  we  expect,  for  the  very  sharp  7.3-meV  mode, 
that  o)  is  near  the  real  axis  so  that  Eq.  (19)  gives  a  reason¬ 
able  representation  of  the  atomic  motion  in  this  case.  It 
has  been  argued  that  this  particular  resonant  mode  plays 
an  essential  role  in  the  dissociation  of  Zn-0  pairs.14  In 
particular,  the  atomic  motions  which  occur  during  the 
early  stages  of  the  dissociation  may  be  as  shown  in  Fig.  6. 
It  is  seen  that  the  amplitude  of  Zn-atom  motion  is  roughly 
3  times  that  of  oxygen.  Thus,  we  conclude  that  in  the  pair 
dissociation  it  is  the  Zn  atom  which  jumps  to  an  intersti¬ 
tial  site  and  diffuses  away. 

The  vibrational  modes  of  NN  Cd-O  pairs  have  been  ob¬ 
served  in  optical  spectra,1719-20  and  can  be  qualitatively 
understood  on  the  basis  of  our  calculations.  The  Cd-O 
modes  will  be  similar  to  the  Zn-O  modes,  except  that 
those  vibrations  involving  significant  Cd  motion  will  be 
shifted  down  in  energy.  Replicas  of  the  electronically  al¬ 
lowed  A  line  should  have  A\  symmetry  for  the  point 
group  C}„,  i.e.,  the  A ,  modes  and  the  A ,  component  of  the 
T2  modes  for  isolated  oxygen.19  In  emission,  Henry, 
Dean,  Thomas,  and  Hopfield  observe  local  modes  at  7.0, 
47.3,  and  49.7  meV.20  The  7.0-meV  mode  is  similar  to  the 
6.0-meV  Zn-O  mode  discussed  above.  Thus,  this  mode  in¬ 
volves  mainly  Cd  motion  in  agreement  with  the  observed 


FIG.  6.  "Eigenvector”  for  the  7.3-meV  T2  mode.  The  atomic 
displacements  are  indicated  for  the  O  and  its  four  Ga  neighbors. 
The  Ga,4,  and  O  displacements  (in  the  ratio  3.4:1)  both  lie  along 
the  bond  direction.  The  eigenvector  has  Cj,  symmetry  relative 
to  the  O— Ga<4i  bond.  With  Zn  in  place  of  the  Gaui  atom,  this  is 
the  predominant  motion  involved  in  the  dissociation  of  Zn-O 
pairs  (see  text). 


isotope  shift.17  The  optical  resonances  at  47.3  and  49.7 
meV  are  probably  associated  with  the  defect  B  and  a , (P) 
modes,  thus  involving  mainly  P-shell  motion.  From  the 
intensities  of  the  observed  phonon  replicas  we  conclude 
that  in  the  optical  transition  the  major  atomic  motion  is 
that  of  the  Cd  atom  moving  relative  to  the  almost  station¬ 
ary  O  atom.  A  number  of  T3  phonons  are  reported20 
which  produce  replicas  of  the  forbidden  B  line.  TTieoreti- 
cally,  we  expect  these  to  arise  from  the  E  modes  and  the  E 
components  of  the  T ,  and  T2  modes  for  isolated  oxygen. 
As  shown  in  Table  I  we  do  compute  various  E  and  T\ 
modes.  All  these  modes  are  very  sharp  (width  of  —0. 1 
meV)  and  are  tied  to  singularities  in  the  bulk  DOS.  These 
modes  could  account  for  some  of  the  observed  T 3  modes, 
and  the  E  components  of  the  T2  mode  at  43.1  meV  may 
also  contribute.  Our  optical-phonon  energies  really  are 
not  accurate  enough  to  make  a  definitive  identification  of 
these  observed  T 3  phonons. 

Let  us  consider  what  our  calculations  tell  us  about  the 
nature  of  the  oxygen  defect.  In  particular,  two  opposing 
models21,22  presently  exist  for  this  defect,  and  it  is  desir¬ 
able  to  try  to  choose  between  them  based  on  our  results. 
First,  the  defect  we  are  considering  in  our  calculations  is 
only  a  model,  without  the  same  detailed  properties  of  the 
actual  oxygen  defect  although  their  general  properties 
should  be  similar.  The  most  important  result  from  our 
calculations  is  the  O-Iattice  effective  force  constant,  3.5 
eV/A2.  This  number  is  accurate  to  within,  perhaps,  10%. 
and  so  provides  a  test  of  the  opposing  models.  Unfor¬ 
tunately,  it  is  very  difficult  to  accurately  obtain  the  defect 
force  constant  from  the  existing  theoretical  models.  To 
date,  Baraff  et  al. 15  have  attempted  such  a  calculation  for 
the  classical  model6,22  and  have  obtained  an  O-Ga  stretch¬ 
ing  constant  of  4.4  eV/A2,  corresponding  to  an  effective 
force  constant  of  4(4.41/3  =  5.9  eV/A2,  a  result  which  is 
somewhat  too  large.  Nobody  has  yet  attempted  a  similar 
computation  for  the  weak-bonding  model,21  although  the 
computed  bond  strength  should  be  significantly  different 
for  the  two  models  as  they  contain  different  numbers  of 
bonding  and  antibonding  electrons.23  Morgan  has  recently 
argued  that  the  observed  structure  near  25  meV  in  the  O0 
phonon  sidebands  is  due  to  a  dynamic  Jahn-Teller  distor¬ 
tion  of  the  O0  defect.23  Although  such  a  distortion  could 
occur,  our  calculations  demonstrate  that  it  is  not  necessary 
to  invoke  this  mechanism  in  order  to  account  for  the  ob¬ 
served  structure.  With  one  free  parameter  (77 ),  we  have 
obtained  quantitative  agreement  between  experiment  and 
theory  for  those  phonon  energies,  isotope  shifts,  and  inten¬ 
sities. 

One  basis  for  the  weak-bonding  model21  is  the  19-meV 
at(Ga)  phonon,  present  in  donor-acceptor— pair  (DAP' 
spectra,1-3  but  apparently  absent  in  photoluminescence- 
excitation  (PLE)  spectra.2  According  to  our  calculations, 
this  mode  is  rather  broad  as  seen  in  the  Ga-shell  A , 
LDOS  for  O  on  a  P  site,  shown  in  Fig.  4.  Clearly,  this 
resonance  is  distributed  throughout  the  acoustic  branch, 
so  let  us  take  its  "width”  to  be  5  meV.  This  width  has  no 
consequence  in  the  DAP  spectrum  due  to  the  greater- 
than-10-meV  broadening  caused  by  a  distribution  of  pair 
separations.  But  in  the  PLE  spectrum,  the  no-phonon  line 
is  relatively  sharp,  so  that  the  ratio  of  peak  heights  be- 
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tween  an  a  |  tGa)-phonon— assisted  line  and  a  no-phonon 
line  is  significantly  reduced  i rom  the  DAP  spectrum.  Ad¬ 
ditional  broadening  for  the  phonon-assisted  line  is  also 
possible;  the  20.4-meV  Ga-shell  peak  in  the  A  |  LDOS  for 
O  on  a  P  site  shown  in  Fig.  4  represents  mainly  NN  Ga 
motion.  Motion  of  further-lying  Ga  shells  will  be  concen¬ 
trated  at  energies  closer  to  27  meV,  and  the  actual  phonon 
sideband  probably  involves  some  combination  of  these 
motions.  We  feel  that  a  calculation  of  these  broadening 
effects  would  be  useful  in  order  to  quantify  the  apparent 
“anomaly.”21 

V,  CONCLUSION 

The  major  results  of  this  paper  are  concerned  with  the 
existence  of  defect  vibrational  modes,  including  both 
strictly  localized  and  resonant  modes.  It  is  impossible  to 
determine  the  number  of  such  modes  by  group  theory 
alone  since  the  existence  of  a  mode  is  determined  by  a 
quantitative  parameter:  the  degree  of  localization  around 
the  defect.  Thus,  it  is  necessary  to  undertake  a  computa¬ 
tion  of  the  vibrational  modes  of  a  defect.  The  theory  we 
use  for  this  puppose  is  a  phenomenological  one  in  which 
the  force  constants  of  the  perfect  lattice  and  those  of  the 
defect  are  all  chosen  to  match  experiment.  Bulk  phonons 
axe  described  by  a  1 5-parameter  model  which  includes  NN 
and  2NN  plus  electrostatic  interactions.  A  defect  is  then 
introduced  by  perturbing  the  mass  of  an  atom  and  per¬ 
turbing  its  NN  and  2NN  force  constants.  A  Green’s- 
function  technique  is  used  to  evaluate  the  LDOS  near  the 
defect,  and  thereby  determine  the  presence  of  defect 
modes. 

Our  theoretical  results  are  compared  with  experiment 
for  the  GaP:Op  defect.  We  use  a  two-parameter  descrip¬ 
tion  of  the  short-range  interaction  perturbation  near  the  O 
site.  By  a  suitable  choice  of  these  two  defect  force  con¬ 
stants  we  have  obtained  quantitative  agreement  with  ob¬ 
served  phonon  modes  involving  0+  and  O0  defects.  For 
O-neighbor  (i.e.,  O-Ga  and  O-P)  and  Ga-P  (back-bond) 
short-range  defect  force  constants  at  38%  and  93%  of 
their  respective  bulk  values,  our  calculations  predict  two 
A  |  modes  and  two  T2  modes  which  are  reasonably  local¬ 
ized  at  or  near  the  defect,  and  also,  other  modes  which  are 
more  bulldike  in  nature.  We  summarize  the  character  of 
the  four  most  localized  modes  as  follows. 

(1)  O-atom  7*2  mode.  This  mode  involves  mainly  O- 
atom  motion,  and  this  motion  is  localized  near  the  defect 
for  practically  any  value  of  the  defect  force  constants.  To 
match  the  observed  phonon  energies  near  25  meV,  our 
choice  of  defect  force  constants  correspond  to  an  O-iattice 
effective  force  constant  (taking  all  long-range  and  short- 
range  interactions  into  account)  of  3.5  eV/A2,  which  is 
25%  of  the  bulk,  P-lattice  value.  This  force  constant  may 
be  compared  with  that  obtained  by  assuming  that  the  ob¬ 
served  phonon  is  completely  localized  on  the  oxygen  im¬ 
purity,  (25  meV/rt)2m0  =  2.4  eV/A2.  Our  calculations 
demonstrate  that  the  observed  two-peak  structure  in  this 
resonant  mode  is  due  to  coupling  with  LA  bulk  phonons. 

(2)  Ga-shell  Tj  mode.  This  mode  involves  mainly  vec¬ 
torlike  motion  of  the  shell  of  Ga  atoms  surrounding  the  O 
impurity.  The  mode  originates  from  the  TA  peak  in 


bulk-phonon  DOS,  and  is  significantly  localized  around 
the  defect  for  values  of  the  O-neighbor  force  constants 
which  arc  about  80%  or  less  of  the  bulk  values.  Experi¬ 
mentally,  this  mode  is  not  seen  in  spectra  involving  isolat¬ 
ed  O  impurities,  but  is  seen  when  a  Zn  or  Cd  atom  is  lo¬ 
cated  near  the  O  atom,  thereby  lowering  the  symmetry  of 
the  system. 

(3)  Ga-shell  A ,  mode.  This  mode  is  a  breathing  motion 
of  the  Ga  shell.  The  mode  originates  at  the  peak  in  the 
LA  bulk-phonon  DOS  and  is  localized  for  values  of  the 
O-neighbor  force  constants  which  are  about  80%  or  less 
of  the  bulk  values.  This  resonant  mode  is  relatively  broad 
with  a  width  of  about  5  meV  (in  agreement  with  experi¬ 
ment). 

(4)  P-shell  A  i  mode.  This  mode  mainly  involves  breath¬ 
ing  motion  of  the  P  shell  (2NN’s)  surrounding  the  O 
atom.  The  mode  originates  at  the  LO  peak  in  the  bulk 
DOS,  and  appears  when  the  Ga-P  back-bond  interactions 
are  reduced  to  95%  or  less  of  the  bulk  values. 

For  all  of  the  four  defect  modes  discussed  above,  we 
feel  that  our  theoretical  description  is  fairly  close  to  exper¬ 
iment  in  terms  of  energies  and  eigenvectors.  We  have 
demonstrated  that  the  0+  and  O0  defects  are  weakly 
bound  in  GaP,  as  expected,  and  that  the  Ga— P  back 
bonds  appear  to  be  weakened.  Our  predicted  two-peak, 
O-atom  T2  mode  shows  that  the  experimental  results  for 
the  O0  defect  are  consistent  with  a  T^-symmetric 
geometry.  Finally,  our  calculations  identify  the  predom¬ 
inant  atomic  motion  involved  in  the  dissociation  of  Zn-O 
pairs  in  GaP. 
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APPENDIX 

We  have  determined  the  collective  coordinates  Q  for  the 
12-atom,  2NN  shell  seen  by  a  substitutional  defect  in  a 
zinc-blende  lattice.  The  collective  coordinates  form  bases 
for  the  irreducible  representations  of  the  symmetry  group 
of  the  2NN  shell,  namely  Oh.  The  equilibrium  positions 
of  the  12  atoms  in  the  2NN  shell  are  given  (in  units  of 
a0/2)  by 

Fi=  —  r*7  =  ( 0, 1 , 1 )  ,  f2=-F,  =  (  1,0,1)  , 

Fj=  —  r9  =  (  1, 1,0)  ,  r4=  —  F|0=(  —  1, 1,0)  ,  (AD 

r3=  —  F, , =(0,  — 1,1)  ,  F6=  -  Fi2  =  (  1,0,  - 1 )  , 

where  a0  is  the  edge  length  of  the  conventional  unit  cell. 
The  36  collective  coordinates  Q  are  related  by  an  orthogo¬ 
nal  transformation  to  the  36  Cartesian  displacements. 
(x(,yf,z(),  of  atom  i  from  equilibrium,  where 
/  =  1, . . . ,  12.  The  explicit  form  of  this  orthogonal 
transformation  is  given  in  Table  II.  In  this  work,  the 
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TABLE  II.  36  collective  coordinates  Q  for  the  12-atom,  2NN  shell.  The  (unnormalized)  components  associated  with  the  atoms 

i  =  1,  ....  6  are  given  explicitly.  The  components  associated  with  the  remaining  atoms  (  /  =  7 . 12)  are  obtained  by  multiplying 

the  entire  row  of  coefficients  by  the  sign  indicated.  The  irreducible  representation  of  the  Td  group  according  to  which  the  collective 
coordinates  transform  are  also  indicated. 
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